SEARCH · Engineering Papers

Results for “Matched”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.

quantum computing, numerical analysis↗

DOE OSTI2022

Approaching PetaVolts per Meter Plasmonics Using Structured Semiconductors

A newly uncovered class of plasmons in the strongly excited limit opens access to unprecedented Petavolts per meter electromagnetic fields with wide-ranging, transformative impact. Unlike conventional plasmons, such plasmons are constituted by non-perturbative, large-amplitude oscillations of the ultradense, delocalized free electron Fermi gas inherent in conductive media. Here structured semiconductors doped to have an appropriate conduction electron density are introduced to tune the properties of the Fermi gas for matched excitation of large-amplitude plasmons using readily available electron beams which enables immediate experimental validation. Specifically, an electrostatic, surface “crunch-in” plasmon is collisionlessly excited by the beam launched inside a tube. Strong excitation due to matching results in relativistic oscillations of the electron gas and unravels unique phenomena. Relativistically induced ballistic electron transport comes about due to relativistic multifold increase in the mean free path and also leads to unconventional heat deposition beyond Ohm’s law. This explains the absence of observed damage or solid-plasma formation in past experiments on conductive samples interacting with electron bunches shorter than 10-13 seconds. Furthermore, relativistic momentum leads to copious tunneling of electron gas across the surface, which then crunches inside the tube. Relativistic effects along with large, localized electron density variations underlying these modes necessitate kinetic approach to theoretical and computational modeling. Kinetic model presented here demonstrates experimental viability of observing tens of gigavolts per meter plasmonic fields excited by matching readily available electron beams to plasmons in semiconductors with 1018cm-3 free electron density, and paves the way for Petavolts per meter plasmonics.

42 ENGINEERING↗

DOE OSTI2021

Analysis of Magnetic Field Emissions and Shield Requirements for Interoperating High-Power EV Wireless Charging System

In this study, the magnetic field emission (MFE) is investigated for matching (identical primary and secondary couplers) and interoperating (different couplers) 11 kW wireless charging systems (WCS). The interoperating charging pads generate higher MFE than the matching charging-pads; hence, using interoperable couplers requires a more effective shielding. In this work, different shield topologies are investigated to suppress the MFE for interoperable conditions. Four 11 kW WCSs with the unipolar (uni) and bipolar (bi) transmitter (Tx) and receiver (Rx) pads has been investigated: 1) uniTx with uniRx, 2) uniTx with biRx, 3) biTx with uniRx, and 4) biTx with bi-Rx. The studied systems are designed for 11 kW WPT3 level with Z2 airgap class (140-150 mm) resembling the standard SAE-J2954. The WCSs are simulated in finite element analysis (FEA) and validated experimentally using a 11 kW biTx-biRx system. The experimental results match with the simulation results with <5% error. The FEA results show that a biTx needs a magnetic shield, and a uniTx needs an aluminum shield to meet the MFE under all interoperability and alignment conditions.

Mohammad, Mostak↗

DOE OSTI2023

A Systematic Investigation of the Effects of Standard-Sample Concentration Mismatch during Fe Isotope Measurement by MC-ICP-MS

Advances in multi-collector inductively coupled plasma-mass spectrometry (MC-ICP-MS) have led to the widespread use of iron (Fe) isotopes to elucidate the (bio)geochemical history of a range of environments. To generate Fe isotope ratio measurements, standard-sample bracketing (SSB) is commonly used to correct for instrumental mass bias inherent to MC-ICP-MS. However, SSB is only accurate when sample and isotope standard Fe concentrations match, in addition to the bulk solution matrix. When the Fe concentrations differ, Fe isotope ratio measurement results may be inaccurate, a phenomenon known as the "self-induced matrix effect." Here, this study systematically characterised the self-induced matrix effect for dry plasma Fe isotope ratio measurements on three MC-ICP-MS instruments and three introduction systems. Our extensive dataset indicates that: (1) the degree of mass bias is consistent regardless of MC-ICP-MS front-end design, (2) the degree of mass bias becomes less reproducible as the concentration difference between the sample and bracketing standard increases, and (3) this applies to both pure Fe solutions and solutions from geological materials. This study reinforces the requirement to match bracketing standard and sample concentrations within 10% and provides a correction method for that fall beyond the recommended concentration range to subsequently allow for proper concentration matching during SSB.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

The diphoton $$q_T$$ spectrum at N$$^3$$LL$$^\prime $$ + NNLO

Abstract We present a $$q_T$$ q T -resummed calculation of diphoton production at order N $$^3$$ 3 LL $$^\prime $$ ′ + NNLO. To reach the primed level of accuracy we have implemented the recently published three-loop $${\mathcal {O}}(\alpha _s^3)$$ O ( α s 3 ) virtual corrections in the $$q\bar{q}$$ q q ¯ channel and the three-loop transverse momentum dependent beam functions and combined them with the existing infrastructure of , a code performing resummation at order N $${}^3$$ 3 LL. While the primed predictions are parametrically not more accurate, one typically observes from lower orders and other processes that they are the dominant effect of the next order. We include in both the $$q\bar{q}$$ q q ¯ and loop-induced gg channel the hard contributions consistently together at order $$\alpha _s^3$$ α s 3 and find that the resummed $$q\bar{q}$$ q q ¯ channel without matching stabilizes indeed. Due to large matching corrections and large contributions and uncertainties from the gg channel, the overall improvements are small though. We furthermore study the effect of hybrid-cone photon isolation and hard-scale choice on our fully matched results to describe the ATLAS $$\text {8}$$ 8 TeV data and find that the hybrid-cone isolation worsens agreement at small $$q_T$$ q T compared to smooth-cone isolation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

cuAlign: Scalable Network Alignment on GPU Accelerators

Given two graphs, the objective of network alignment is to find the best one-to-one mapping of vertices in one graph (??) to vertices in the other (??), such that the number of overlaps is maximized. We say that edges(??, ??) ???and(??', ??') ??? are overlapped if ?? is mapped to ??' and ?? is mapped to??'. Network alignment is an important optimization problem with several applications in bioinformatics, computer vision and ontology matching. Since it is an NP-hard problem, efficient heuristics and scalable implementations are necessary. In this work, we introduce a new framework that combines the concepts of intra-network proximity using vertex embedding,Belief Propagation (BP) and approximate weighted matching, and provides qualitative improvements up to22%over state-of-the-art approaches. We also provide scalable implementations on GPU accelerators, demonstrating up to19×speedup for Belief Propagation and 3× speedup for approximate weighted matching relative to previous multithreaded implementation. A combination of combinatorial and algebraic kernels within the network alignment algorithm poses significant hurdles for parallelization. Load imbalance and irregular DRAM traffic limit achievable performance on GPUs. Our novel approach identifies and exploits unique structural proper-ties of the BP-based algorithm and employs code fusion to reduce data movement between different steps of the algorithm. Using a diverse set of inputs, we demonstrate qualitative improvements of our algorithms, and performance gains of our GPU-accelerated implementation. We believe that our work will enable algorithmic improvements and practical applications of network alignment.

Xiang, Lizhi↗

DOE OSTI2022

Mappymatch [SWR 22-38]

Mappymatch is a pure-python package developed and open sourced by the National Renewable Energy Laboratory. It contains a collection of "Matchers" that enable matching a GPS trace (series of GPS coordinates) to a map. As of 8/14/2023, The current matchers are: LCSSMatcher: A matcher that implements the LCSS algorithm described in this paper. Works best with high resolution GPS traces. OsrmMatcher: A light matcher that pings an OSRM server to request map matching results. See the official documentation for more info. ValhallaMatcher: A matcher to ping a Valhalla server for map matching results. Currently supported map formats are: Open Street Maps NOTE: This software was formerly known as YAMM (Yet Another Map-Matcher)

Reinicke, Nicholas↗

DOE OSTI2024

IsoMatchMS

An R package built to match proteoform isotoping distributions to mass spectrum, quantitatively select matches of high quality, and visualize those matches in a trelliscope display.

Degnan, David↗

DOE OSTI2023

Real-Time Multi-Vehicle Multi-Camera Tracking with Graph-Based Tracklet Features

An essential application in intelligent transportation systems is multi-target multi-camera tracking (MTMCT), where the target’s activity is tracked from different cameras. Although the tracking-by-detection scheme is the primary paradigm in MTMCT, the object association information from the video frames is lost. This is mainly because the multi-camera multi-object matching uses the information from the video frames separately. To solve this problem and leverage this association information, we propose an MTMCT framework, where features are built in the form of a graph and a graph similarity algorithm is used to match multi-camera objects. In this paper, we focus on the real-time scenario, where only the past images are used to match an object. Our method achieves an IDF1 score (the ratio of the number of correctly identified objects to the number of ground truth and average objects) of 0.75 with a rate of 14 frames per second (fps).

Engineering↗

DOE OSTI2021

On-chip non-Hermitian optical parametric amplifiers with a large bandwidth

Recently, our groups have introduced the notion of optical parametric amplification based on non-Hermitian phase matching wherein the incorporation of loss can lead to gain in this nonlinear optical process. Previous simulation results using second-order nonlinear optical coupled-mode theory have demonstrated the potential of this technique as an alternative to the stringent phase-matching condition, which is often difficult to achieve in semiconductor platforms. Here we fortify this notion for the case of third-order nonlinearity by considering parametric amplification in silicon nanowires and illustrate the feasibility of these devices by employing rigorous finite-difference time-domain analysis using realistic materials and geometric parameters. Particularly, we demonstrate that by systematic control of the optical loss of the idler in a four-wave mixing process, we can achieve efficient unidirectional energy conversion from the pump to the signal component even when the typical phase-matching condition is violated. Importantly, our simulations show that a signal gain of ∼ <#comment/> 9 d B for a waveguide length of a few millimeters is possible over a large bandwidth of several hundreds of nanometers ( ∼ <#comment/> 600 n m ). This bandwidth is nearly 2 orders of magnitude larger than what can be achieved in the conventional silicon-photonics-based four-wave mixing process.

Ahmed, Asif (ORCID:0000000152234392)↗

DOE OSTI2022

Quadratic strong coupling in AlN Kerr cavity solitons

Photonic platforms with χ (2) nonlinearity offer new degrees of freedom for Kerr frequency comb development. Here, we demonstrate Kerr soliton generation at 1550 nm with phase-matched quadratic coupling to the 775 nm harmonic band in a single AlN microring and thus the formation of dual-band mode-locked combs. In the strong quadratic coupling regime where the χ (2) phase-matching window overlaps the pump mode, the pump-to-harmonic-comb conversion efficiency is optimized. However, the strong quadratic coupling also drastically modifies the Kerr comb generation dynamics and decreases the probability of soliton generation. By engineering the χ (2) phase-matching wavelength, we are able to achieve a balance between high conversion efficiency and high soliton formation rate under the available pump power and microring quality factors. Our numerical simulations confirm the experimental observations. These findings provide guidance on tailoring single-cavity dual-band coherent comb sources.

Gong, Zheng (ORCID:0000000301704870)↗

DOE OSTI2023

Tunable and efficient ultraviolet generation with periodically poled lithium niobate

On-chip ultraviolet (UV) sources are of great interest for building compact and scalable atomic clocks, quantum computers, and spectrometers. However, few material platforms are suitable for integrated UV light generation and manipulation. Of these materials, thin-film lithium niobate offers unique advantages such as sub-micron modal confinement, strong nonlinearity, and quasi-phase matching. Despite these characteristics, its utilization in the UV has remained elusive because of the substantial sensitivity of standard quasi-phase matching to fabrication imperfections, the photorefractive effect, and relatively large losses in this range. Here, we present efficient (197 ± 5%/W/cm 2 ) second harmonic generation of UV-A light in a periodically poled lithium niobate nanophotonic waveguide. We achieve on-chip UV powers of ~30 µW and linear wavelength tunability using temperature. These results are enabled with large cross section waveguides, which leads to first-order UV quasi-phase-matching with relatively long poling periods (>1.5 µm). By varying the poling period, we have achieved the shortest reported wavelength (355 nm) generated through frequency doubling in thin-film lithium niobate. Furthermore, our results open up new avenues for UV on-chip sources and chip-scale photonics through compact frequency-doubling of common near-IR laser diodes.

42 ENGINEERING↗

DOE OSTI2023

Enhancement of second-order optical nonlinearities and nanoscale periodic domain patterning in ferroelectric boron-substituted aluminum nitride thin films

The discovery and development of CMOS-compatible, nonlinear optical materials is essential to produce integrated photonic devices with advanced functionalities. AlN is a strong candidate for on-chip device demonstration due to its intrinsic second-order optical nonlinearities, large bandgap, and well-established fabrication techniques. However, AlN is not easily phase matched for the largest coefficient d 33 ; the coefficients that could potentially be dispersion phase-matched, d 31 and d 15 , have weak nonlinearities. This work investigates ferroelectric Al 1-x B x N (x = 0 to 0.11) for viability as a large bandgap nonlinear optical material with unique suitability towards ultraviolet light generation using second harmonic generation. The linear and nonlinear optical properties are characterized accounting for material anisotropy. With increasing B concentration, a large enhancement from near negligible values to d 31 = 0.9 ± 0.1 pm/V and d 15 = 1.2 ± 0.1 pm/V is observed. This compares favorably to other large bandgap materials like β-Ba(BO 2 ) 2 , where the largest nonlinear coefficient is d 22 ~ 2.3 pm/V at 800 nm. This is accompanied by a change in the bandgap from 6.1 eV to 5.8 eV as B substitution goes from 0 to 11%. A periodically poled, quasi-phase-matched ferroelectric domain pattern with 400 nm domain size and a wall roughness of <16 nm is demonstrated.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Synthetic method of analogues for emerging infectious disease forecasting

The Method of Analogues (MOA) has gained popularity in the past decade for infectious disease forecasting due to its non-parametric nature. In MOA, the local behavior observed in a time series is matched to the local behaviors of several historical time series. The known values that directly follow the historical time series that best match the observed time series are used to calculate a forecast. This non-parametric approach leverages historical trends to produce forecasts without extensive parameterization, making it highly adaptable. However, MOA is limited in scenarios where historical data is sparse. This limitation was particularly evident during the early stages of the COVID-19 pandemic, where the emerging global epidemic had little-to-no historical data. In this work, we propose a new method inspired by MOA, called the Synthetic Method of Analogues (sMOA). sMOA replaces historical disease data with a library of synthetic data that describe a broad range of possible disease trends. This model circumvents the need to estimate explicit parameter values by instead matching segments of ongoing time series data to a comprehensive library of synthetically generated segments of time series data. We demonstrate that sMOA has competitive performance with state-of-the-art infectious disease forecasting models, out-performing 78% of models from the COVID-19 Forecasting Hub in terms of averaged Mean Absolute Error and 76% of models from the COVID-19 Forecasting Hub in terms of averaged Weighted Interval Score. Additionally, we introduce a novel uncertainty quantification methodology designed for the onset of emerging epidemics. Developing versatile approaches that do not rely on historical data and can maintain high accuracy in the face of novel pandemics is critical for enhancing public health decision-making and strengthening preparedness for future outbreaks.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Effects of soil particles and convective transport on dispersion and aggregation of nanoplastics via small-angle neutron scattering (SANS) and ultra SANS (USANS)

Terrestrial nanoplastics (NPs) pose a serious threat to agricultural food production systems due to the potential harm of soil-born micro- and macroorganisms that promote soil fertility and ability of NPs to adsorb onto and penetrate into vegetables and other crops. Very little is known about the dispersion, fate and transport of NPs in soils. This is because of the challenges of analyzing terrestrial NPs by conventional microscopic techniques due to the low concentrations of NPs and absence of optical transparency in these systems. Herein, we investigate the potential utility of small-angle neutron scattering (SANS) and Ultra SANS (USANS) to probe the agglomeration behavior of NPs prepared from polybutyrate adipate terephthalate, a prominent biodegradable plastic used in agricultural mulching, in the presence of vermiculite, an artificial soil. SANS with the contrast matching technique was used to study the aggregation of NPs co-dispersed with vermiculite in aqueous media. We determined the contrast match point for vermiculite was 66% D 2 O / 33% H 2 O. At this condition, the signal for vermiculite was ~50–100%-fold lower that obtained using neat H 2 O or D 2 O as solvent. According to SANS and USANS, smaller-sized NPs (50 nm) remained dispersed in water and did not undergo size reduction or self-agglomeration, nor formed agglomerates with vermiculite. Larger-sized NPs (300–1000 nm) formed self-agglomerates and agglomerates with vermiculite, demonstrating their significant adhesion with soil. However, employment of convective transport (simulated by ex situ stirring of the slurries prior to SANS and USANS analyses) reduced the self-agglomeration, demonstrating weak NP-NP interactions. Convective transport also led to size reduction of the larger-sized NPs. Therefore, this study demonstrates the potential utility of SANS and USANS with contrast matching technique for investigating behavior of terrestrial NPs in complex soil systems.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"

This data package contains the associated data and scripts for Nagamoto, E., Ombadi, M., Ciulla, F. et al. Widespread drought-driven declines in streamflows and water quality in the Upper Colorado River Basin during 1998-2022. Commun Earth Environ 7, 734 (2026). https://doi.org/10.1038/s43247-026-03890-5. This purpose of this study was to investigate the impact of the 21st century drought on water quantity and quality at catchments throughout the Upper Colorado River Basin (UCRB). We used stream flow, water temperature, specific conductance, air temperature, precipitation, and catchment attribute data for over 200 sites in the UCRB, collected from the National Water Information System using Basin3D (Varadharajan, 2023), GAGESII (Falcone, 2010), and the Google Earth Engine. We identified years of severe drought between 1998 and 2022 using the Standardized Precipitation Evaporation Index (SPEI), then calculated the relative change percentage of the stream flow, water temperature, and specific conductance from drought versus non-drought years. We used the attribute information from GAGESII to investigate what physical traits of catchments are associated streamflow vulnerability (greater relative change) or resilience to drought. We used land cover data from the National Land Cover Database (USGS, 2024) to assess any changes to physical attributes that may not be represented in the static attributes information in GAGESII. To increase data availability, we modeled stream temperature using methods from Willard, 2023. While the study period is water years 1998 to 2022, the raw water quantity and quality data extends to 1950 and the meteorological data extends to 1980. The data and code can be downloaded via the UCRB_drought.zip. Within the zip, the files are organized as follows: - INPUTS: Contains all input data used in UCRB_Drought_Workflow.ipynb - OUTPUTS: Contains all intermediate data created from UCRB_Drought_Workflow.ipynb as well as final products including the calculated Standardized Evapotranspiration Index (SPEI) - climatic_variables: The code used to collect meteorologic data from Google Earth Engine - feature_importance: The code used for the catchment attributes analysis - preprocessing: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - pyeto: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - calculations: Code used in UCRB_Drought_Workflow_Impacts.ipynb - plotting: Code used in UCRB_Drought_Workflow_Impacts.ipynb - README.md - UCRB_Drought_Workflow_Preprocessing.ipynb: The code used to prep raw data for the analysis - UCRB_Drought_Workflow_Impact.ipynb: The code which uses the prepped raw data for analysis, and plots all figures - requirements_ucrb-drought_v2.yml: The requirements file to create a virtual environment and Jupyter Lab kernel to run the code The INPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_RAW" folder contains raw data for streamflow, water temperature, and specific conductance in a ".h5" file. The "NLCD_RAW" folder contains ".csv" files with annual land cover percentages for counties within the UCRB. The "MET_RAW" folder contains a ".csv" file with monthly meteorological data (air temperature and precipitation) for the sites in the UCRB which was obtained from code in the climatic_variables folder. The "GAGESII" folder contains ".csv" files with physical catchment attribute variables for catchments across the country. The "WT_LSTM_data" folder contains ".csv" files with calculated WT (Willard, 2023) and the associated RMSEs. The "Upper_Colorado_River_Basin_Boundary" folder contains geographic data including a shapefile for plotting in the UCRB_Drought_Workflow.ipynb. The "RESERVOIRS_RAW" folder contains ".csv" files for each reservoir in the UCRB with daily reservoir storage. There are also two files in the INPUTS folder that have combined reservoir storage data and reservoir metadata. The OUTPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_data" folder contains a folder "Water_year" with the associated cleaned data, metadata, and data availability information in ".csv" files, a folder "Median_Relchange" with the relative change comparing drought to non-drought years in ".csv" files, and a folder "Peak95_Min5_Relchange" that has ".csv" files for the relative change in peak (95th %) and minimum (5th %) variables. The "NLCD_data" folder contains the difference in land cover from the beginning to end of the study period and the percentage of the county that is within UCRB bounds can be found in Nagamoto et al (2025)). The "MET_data" folder contains separated monthly air temperature and precipitation data and the calculated PET in ".csv" files. The "SPEI_data" folder contains ".csv" files with calculated SPEI values (one restricted to the study period and the other with information from the entire MET data period). The "Paper_Tables" folder contains two ".csv" files containing site information and data availability and information about the GAGESII trait aggregated categories. The base directory includes the file “flmd.csv” for a list and description of all files and the file “dd.csv” for data dictionaries. Scripts for preprocessing, analysis, and figure generation are located in the associated GitHub repository found at [https://github.com/iNAIADS/drought-impacts/tree/develop/UCRB-drought]. UPDATE 1: Title and code file updated to match submitted manuscript 10-15-2025. UPDATE 2: Code and data files updated to match revised manuscript 3-4-2026. UPDATE 3: Code and data files updated to match revised manuscript 6-7-2026. ** NOTE: DD and FLMD have not been updated yet. UPDATE 4: Added associated Manuscript information and DD and FLMD have been updated. To cite this code, please use the following BibTeX: @misc{nagamoto2025drought, author = {Emily Nagamoto and Fabio Ciulla and Mohammad Ombadi and Jared Willard and Rosemary Carroll and Charuleka Varadharajan}, title = {Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"}, year = {2025}, doi = {10.15485/2551894}, publisher = {ESS-DIVE Repository}, url = {https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2551894} }

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗

DOE OSTI2025

Towards Improving luminosity using optics tuning and data-driven methods

The results of Run 24 experiments at Relativistic Heavy Ion Collider (RHIC) for improving luminosity using optics tuning are presented in this study. In the first experiment, MADx matching was used to output magnet strengths corresponding to specific s star movements around Interaction Region 8 (IR8). The corresponding Zero Degree Calorimeter (ZDC) signal was measured in place of luminosity, and Bayesian Optimization aids search of optimal movements. It was found that values retrieved from matching were inaccurate, resulting in negative feedback loops. The second experiment focused on calculating accurate s star movements. The matching method was replaced with a linear sensitivity matrix, directly relating optics to power supply, and its null space was used to fit constraints such as hysteresis effects. At the experiment, beam losses were observed at collimators around boundary of IR8, which were fixed for the third experiment. Dynamic mode decomposition was also introduced to improve quality of turn-by-turn (TBT) data as well as accuracy and consistency of optics measurements at IR8. These improvements will be tested in the experiment of next RHIC run for luminosity optimization.

Accelerator Physics↗