Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large system simulation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Strong Dark Matter Self-interactions Diversify Halo Populations within and surrounding the Milky Way

Abstract We perform a high-resolution cosmological zoom-in simulation of a Milky Way (MW)–like system, which includes a realistic Large Magellanic Cloud analog, using a large differential elastic dark matter self-interaction cross section that reaches ≈100 cm 2 g −1 at relative velocities of ≈10 km s −1 , motivated by the diverse and orbitally dependent central densities of dwarf galaxies within and surrounding the MW. We explore the effects of dark matter self-interactions on satellite, splashback, and isolated halos through their abundance, central densities, maximum circular velocities, orbital parameters, and correlations between these variables. We use an effective constant cross section model to analytically predict the stages of our simulated halos’ gravothermal evolution, demonstrating that deviations from the collisionless R max – V max relation can be used to select deeply core-collapsed halos, where V max is a halo’s maximum circular velocity, and R max is the radius at which it occurs. We predict that a sizable fraction (≈20%) of subhalos with masses down to ≈10 8 M ⊙ is deeply core collapsed in our SIDM model. Core-collapsed systems form ≈10% of the isolated halo population down to the same mass; these isolated, core-collapsed halos would host faint dwarf field galaxies with extremely steep central density profiles. Finally, most halos with masses above ≈10 9 M ⊙ are core-forming in our simulation. Our study thus demonstrates how self-interactions diversify halo populations in an environmentally dependent fashion within and surrounding MW-mass hosts, providing a compelling avenue to address the diverse dark matter distributions of observed dwarf galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

MIRACL Co‐Simulation platform for control and operation of distributed wind in microgrid

Abstract This paper presents a co‐simulation platform (CSP) developed as a part of the “Microgrids, Infrastructure Resilience and Advanced Controls Launchpad (MIRACL)” project, hereafter called MIRACL‐CSP, to allow simulation‐based testing, demonstration, and evaluations of distributed wind under various grid operating conditions. MIRACL‐CSP provides modular interactions among the power distribution system, distributed wind, and utility decision‐making framework. A hierarchical engine for large‐scale infrastructure co‐simulations (HELICS) is used as the core engine of MIRACL‐CSP to establish time and information coordination among the MIRACL‐CSP modules. In this study, MIRACL‐CSP performance is demonstrated using the IEEE 123‐node test distribution grid modelled in GridLAB‐D and the utility decision support application modelled in Python. The functionality of MIRACL‐CSP is demonstrated through various grid operational scenarios in a microgrid and networked microgrid environment. A co‐simulation study is described that demonstrated MIRACL‐CSP capabilities for the microgrid operational scenario in the IEEE 123‐node test system. MIRACL‐CSP is a generic platform that facilitates distributed wind research for different test systems, applications, and valuation of distributed wind.

17 WIND ENERGY↗
DOE OSTI2023

Nuclear–Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics

Simulating the nuclear–electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear–electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear–electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron–proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. Furthermore, these examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Towards a Techno-Economic Analysis of PCM Integrated Hybrid HVAC Systems

Thermal end uses dominate building energy consumption and are a major driver of peak demand. As heating is electrified, peak electrical power required will surge, prompting a need for innovative HVAC system designs and controls.These designs must incorporate novel technologies at the component level and new integration techniques at the system level. One such possibility involves the addition of thermal energy storage (TES) in heating and cooling equipment using a phase change material (PCM) heat exchanger. Here, thermal energy storage via phase change can be used to shift the HVAC system loads to times of lower electricity cost, reduced carbon intensity, and greater energy efficiency. Most of the current utilization of PCM in buildings involves passive components. By actively controlling when heat is stored and released from PCM, we can optimize the building HVAC system to cost-effectively meet consumer needs with the flexibility to draw on renewable energy resources when they are abundant and available. While this combination of technologies is promising in theory, simulation-based evaluation of a prototype can be difficult due to the modeling requirements at the component level and the large number of possible configurations and operating modes at the system level. To conduct this evaluation, we use the Modelica language for modeling and simulation because it enables users to represent the important physics of the problem, interchange and rearrange components in an efficient manner, and implement a range of control configurations. In this work, we considered three case studies: a portable building, a large commercial retail store, and a multifamily residential apartment unit. Each of these employs a different system design, ranging from a single package vertical unit incorporating PCM to a central plant with independent heat pump, evaporative cooling, and thermal energy storage components. This paper describes the technologies in question, presents modeling at the component and system levels, and demonstrates building energy and demand charge cost savings with local time-of-use tariffs in a hot-dry climate.

Helmns, Dre↗
DOE OSTI2021

Modeling Framework for Bulk Electric Grid Impacts from HEMP E1 and E3 Effects (Tasks 3.1 Final Report)

This report presents a framework to evaluate the impact of a high-altitude electromagnetic pulse (HEMP) event on a bulk electric power grid. This report limits itself to modeling the impact of EMP E1 and E3 components. The co-simulation of E1 and E3 is presented in detail, and the focus of the paper is on the framework rather than actual results. This approach is highly conservative as E1 and E3 are not maximized with the same event characteristics and may only slightly overlap. The actual results shown in this report are based on a synthetic grid with synthetic data and a limited exemplary EMP model. The framework presented can be leveraged and used to analyze the impact of other threat scenarios, both manmade and natural disasters. This report d escribes a Monte-Carlo based methodology to probabilistically quantify the transient response of the power grid to a HEMP event. The approach uses multiple fundamental steps to characterize the system response to HEMP events, focused on the E1 and E3 components of the event. 1) Obtain component failure data related to HEMP events testing of components and creating component failure models. Use the component failure model to create component failure conditional probability density function (PDF) that is a function of the HEMP induced terminal voltage. 2) Model HEMP scenarios and calculate the E1 coupled voltage profiles seen by all system components. Model the same HEMP scenarios and calculate the transformer reactive power consumption profiles due to E3. 3) Sample each component failure PDF to determine which grid components will fail, due to the E1 voltage spike, for each scenario. 4) Perform dynamic simulations that incorporate the predicted component failures from E1 and reactive power consumption at each transformer affected by E3. These simulations allow for secondary transients to affect the relays/protection remaining in service which can lead to cascading outages. 5) Identify the locations and amount of load lost for each scenario through grid dynamic simulation. This can be an indication of the immediate grid impacts from a HEMP event. In addition, perform more detailed analysis to determine critical nodes and system trends. 6) To help realize the longer-term impacts, a security constrained alternating current optimal power flow (ACOPF) is run to maximize critical load served. This report describes a modeling framework to assess the systemic grid impacts due to a HEMP event. This stochastic simulation framework generates a large amount of data for each Monte Carlo replication, including HEMP location and characteristics, relay and component failures, E3 GIC profiles, cascading dynamics including voltage and frequency over time, and final system state. This data can then be analyzed to identify trends, e.g., unique system behavior modes or critical components whose failure is more likely to cause serious systemic effects. The proposed analysis process is demonstrated on a representative system. In order to draw realistic conclusions of the impact of a HEMP event on the grid, a significant amount of work remains with respect to modeling the impact on various grid components.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Influences of lidar scanning parameters on wind turbine wake retrievals in complex terrain

Abstract. Scanning lidars enable the collection of spatially distributed measurements of turbine wakes and the estimation of wake properties such as magnitude, extent, and trajectory. Lidar-based characterizations, however, may be subject to distortions due to the observational system. Distortions can arise from the resolution of the measurement points across the wake, the projection of the winds onto the beam, averaging along the beam probe volume, and intervening evolution of the flow over the scan duration. Using a large-eddy simulation and simulated measurements with a virtual lidar model, we assess how scanning lidar systems may influence the properties of the retrieved wake using a case study from the Perdigão campaign. We consider three lidars performing range-height indicator sweeps in complex terrain, based on the deployments of lidars from the Danish Technical University (DTU) and German Aerospace Center (DLR) at the Perdigão site. The unwaked flow, measured by the DTU lidar, is well-captured by the lidar, even without combining data into a multi-lidar retrieval. The two DLR lidars measure a waked transect from different downwind vantage points. In the region of the wake, the observation system reacts to the smaller spatial and temporal variations of the winds, allowing more significant observation distortions to arise. While the measurements largely capture the wake structure and trajectory over its 4–5 D extent, limited spatial resolution of measurement points and volume averaging lead to a quicker loss of the two lobes in the near wake, smearing of the vertical bounds of the wake (< 30 m), wake center displacements up to 10 m, and dampening of the maximum velocity deficit by up to a third. The virtual lidar tool, coupled with simulations, provides a means for assessing measurement capabilities in advance of measurement campaigns.

17 WIND ENERGY↗
DOE OSTI2020

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Driving force and pathway in polyelectrolyte complex coacervation

There is notable discrepancy between experiments and coarse-grained model studies regarding the thermodynamic driving force in polyelectrolyte complex coacervation: experiments find the free energy change to be dominated by entropy, while simulations using coarse-grained models with implicit solvent usually report a large, even dominant energetic contribution in systems with weak to intermediate electrostatic strength. Here, using coarse-grained, implicit-solvent molecular dynamics simulation combined with thermodynamic analysis, we study the potential of mean force (PMF) in the two key stages on the coacervation pathway for symmetric polyelectrolyte mixtures: polycation–polyanion complexation and polyion pair–pair condensation. We show that the temperature dependence in the dielectric constant of water gives rise to a substantial entropic contribution in the electrostatic interaction. By accounting for this electrostatic entropy, which is due to solvent reorganization, we find that under common conditions (monovalent ions, room temperature) for aqueous systems, both stages are strongly entropy-driven with negligible or even unfavorable energetic contributions, consistent with experimental results. Furthermore, for weak to intermediate electrostatic strengths, this electrostatic entropy, rather than the counterion-release entropy, is the primary entropy contribution. From the calculated PMF, we find that the supernatant phase consists predominantly of polyion pairs with vanishingly small concentration of bare polyelectrolytes, and we provide an estimate of the spinodal of the supernatant phase. Finally, we show that prior to contact, two neutral polyion pairs weakly attract each other by mutually induced polarization, providing the initial driving force for the fusion of the pairs.

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

Advancing representations of turbulence in Earth System Models [Slides]

We have created a suite of coastally focused Large Eddy Simulation (LES), including a first ever investigation of the influence of the ocean bottom on turbulence. We have added important forcing not previously considered in LES models. We have evaluated a new more physically realistic and advanced closure against our LES database of simulations. Test cases will constrain constants in the new model and help us transition the mixing model from the deep ocean to the coastal ocean for the first time.

58 GEOSCIENCES↗
DOE OSTI2020

Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method

Recent advances in selected configuration interaction methods have made them competitive with the most accurate techniques available and, hence, creating an increasingly powerful tool for solving quantum Hamiltonians. In this work, we build on recent advances from the adaptive sampling configuration interaction (ASCI) algorithm. We show that a useful paradigm for generating efficient selected CI/exact diagonalization algorithms is driven by fast sorting algorithms, much in the same way iterative diagonalization is based on the paradigm of matrix vector multiplication. We present several new algorithms for all parts of performing a selected CI, which includes new ASCI search, dynamic bit masking, fast orbital rotations, fast diagonal matrix elements, and residue arrays. The ASCI search algorithm can be used in several different modes, which includes an integral driven search and a coefficient driven search. The algorithms presented here are fast and scalable, and we find that because they are built on fast sorting algorithms they are more efficient than all other approaches we considered. After introducing these techniques, we present ASCI results applied to a large range of systems and basis sets to demonstrate the types of simulations that can be practically treated at the full-CI level with modern methods and hardware, presenting double- and triple-ζ benchmark data for the G1 data set. The largest of these calculations is Si$_2$H$_6$ which is a simulation of 34 electrons in 152 orbitals. We also present some preliminary results for fast deterministic perturbation theory simulations that use hash functions to maintain high efficiency for treating large basis sets.

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

Adaptive Machine Learning for Robust Diagnostics and Control of Time-Varying Particle Accelerator Components and Beams

Machine learning (ML) is growing in popularity for various particle accelerator applications including anomaly detection such as faulty beam position monitor or RF fault identification, for non-invasive diagnostics, and for creating surrogate models. ML methods such as neural networks (NN) are useful because they can learn input-output relationships in large complex systems based on large data sets. Once they are trained, methods such as NNs give instant predictions of complex phenomenon, which makes their use as surrogate models especially appealing for speeding up large parameter space searches which otherwise require computationally expensive simulations. However, quickly time varying systems are challenging for ML-based approaches because the actual system dynamics quickly drifts away from the description provided by any fixed data set, degrading the predictive power of any ML method, and limits their applicability for real time feedback control of quickly time-varying accelerator components and beams. In contrast to ML methods, adaptive model-independent feedback algorithms are by design robust to un-modeled changes and disturbances in dynamic systems, but are usually local in nature and susceptible to local extrema. In this work, we propose that the combination of adaptive feedback and machine learning, adaptive machine learning (AML), is a way to combine the global feature learning power of ML methods such as deep neural networks with the robustness of model-independent control. We present an overview of several ML and adaptive control methods, their strengths and limitations, and an overview of AML approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Deep Learning-Based Dynamic Modeling of Three-Phase Voltage Source Inverters

Inverter-based resource (IBR) models are necessary to analyze modern power system stability and create effective control strategies. Modeling IBRs in converter-rich power systems is crucial, yet challenging due to the lack of commercial information on converter topologies and control parameters. This paper proposes novel convolutional neural network (CNN)–based data-driven techniques for modeling IBRs, addressing adaptability and proprietary concerns without requiring internal system physics knowledge. The proposed method is tested using real grid-tied commercial IBR transient data and demonstrates effectiveness and accuracy. Furthermore, the developed modeling approach is integrated and implemented in the open-source power distribution simulation and analysis tool, GridLAB-D, to illustrate the potentiality of dynamic analysis of large-scale power systems with high IBRs.

deep learning, artificial intelligence↗
DOE OSTI2020

Quantum Versus Classical Switching Dynamics of Driven Dissipative Kerr Resonators

Here, we report a first-principles study of the driven dissipative dynamics for Kerr oscillators in the mesoscopic regime. This regime is characterized by large Kerr nonlinearity, realized here using the nonlinear kinetic inductance of a large array of Josephson junctions. The experimentally measured nonlinear resonance lineshapes of the junction array modes show significant deviations from steady-state numerical predictions, and necessitate time-dependent numerical simulations indicative of strong measurement-induced dephasing in the system arising from the large cross-Kerr effect between array modes. Analytical and numerical calculations of switching rate corroborate this by showing the emergence of a slow time scale, which is much longer than the linear decay rate and is set by fluctuation-induced switching times in the bistable regime. Furthermore, our analysis shows that the usual quantum-activated escape treatment is inadequate for prediction of the switching rates at large frequency shifts caused by strong nonlinearities, necessitating a quantum treatment that utilizes the full system Liouvillian. Based on our analysis, we identify a universal crossover parameter that delineates the regimes of validity of semiclassical and quantum descriptions, respectively. Our work shows how dynamical switching effects in strongly nonlinear systems provide a platform to study quantum-to-classical transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Agentic traffic intelligence: Augmented human-in-the-loop scenario generation for microscopic traffic simulation

Traditional microscopic traffic simulation generation often relies on static datasets and manual design, limiting its ability to simulate complex conditions easily. This paper presents a novel framework, Agentic Traffic Intelligence, which combines human approval large language models (LLMs), the Real-Twin tool, and multi-agent systems to perform realistic microscopic traffic simulation scenario generation. The proposed framework incorporates human-in-the-loop (HIL) control, retrieval-augmented generation (RAG), and multi-agent control mechanisms. HIL mechanisms are used to guide multiple LLMs focused on attributes for microscopic simulation generation and to improve the interpretability and transparency of LLM execution for users. RAG enhances context extraction by dynamically integrating external knowledge sources for traffic scenario generation foundations. A multi-agent architecture with supervisory control coordinates the interaction of simulation components, including traffic simulators, control logic, and calibration tools. This enables the synthesis of simulation-ready scenarios that reflect dynamic demand profiles and behavior controls. Furthermore, the framework fuses multisource traffic data with unstructured context and supports iterative refinement through interactive user feedback. Validated through microscopic simulation using Simulation of Urban Mobility, the generated scenarios demonstrate high-fidelity network generation with inflow and turn movement and behavioral calibration, offering a robust and efficient tool for stress-testing and optimizing urban mobility systems.

Hierarchical multi-agent control↗
DOE OSTI2020

Benchmarking and parameter sensitivity of physiological and vegetation dynamics using the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) at Barro Colorado Island, Panama

Plant functional traits determine vegetation responses to environmental variation, but variation in trait values is large, even within a single site. Likewise, uncertainty in how these traits map to Earth system feedbacks is large. We use a vegetation demographic model (VDM), the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), to explore parameter sensitivity of model predictions, and comparison to observations, at a tropical forest site: Barro Colorado Island in Panama. We define a single 12-dimensional distribution of plant trait variation, derived primarily from observations in Panama, and define plant functional types (PFTs) as random draws from this distribution. We compare several model ensembles, where individual ensemble members vary only in the plant traits that define PFTs, and separate ensembles differ from each other based on either model structural assumptions or non-trait, ecosystem-level parameters, which include (a) the number of competing PFTs present in any simulation and (b) parameters that govern disturbance and height-based light competition. While single-PFT simulations are roughly consistent with observations of productivity at Barro Colorado Island, increasing the number of competing PFTs strongly shifts model predictions towards higher productivity and biomass forests. Different ecosystem variables show greater sensitivity than others to the number of competing PFTs, with the predictions that are most dominated by large trees, such as biomass, being the most sensitive. Changing disturbance and height-sorting parameters, i.e., the rules of competitive trait filtering, shifts regimes of dominance or coexistence between early- and late-successional PFTs in the model. Increases to the extent or severity of disturbance, or to the degree of determinism in height-based light competition, all act to shift the community towards early-successional PFTs. In turn, these shifts in competitive outcomes alter predictions of ecosystem states and fluxes, with more early-successional-dominated forests having lower biomass. It is thus crucial to differentiate between plant traits, which are under competitive pressure in VDMs, from those model parameters that are not and to better understand the relationships between these two types of model parameters to quantify sources of uncertainty in VDMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluation of Mesoscale Convective Systems in Climate Simulations: Methodological Development and Results from MPAS-CAM over the U.S.

In this study, a process-oriented approach is developed to evaluate warm-season mesoscale convective system (MCS) precipitation and their favorable large-scale meteorological patterns (FLSMPs) over the U.S. This approach features a novel observation-driven MCS-tracking algorithm using infrared brightness temperature and precipitation feature at 12, 25 and 50 km resolution and metrics to evaluate the model large-scale environment favorable for MCS initiation. The tracking algorithm successfully reproduces the observed MCS statistics from a reference 4-km radar MCS database. To demonstrate the utility of the new methodologies in evaluating MCS in climate simulations with mesoscale resolution, the process-oriented approach is applied to two climate simulations produced by the Variable-Resolution Model for Prediction Across Scales coupled to the Community Atmosphere Model physics, with refined horizontal grid spacing at 50 km and 25 km over North America. With the tracking algorithm applied to simulations and observations at equivalent resolutions, the simulated number of MCS and associated precipitation amount, frequency and intensity are found to be consistently underestimated in the Central U.S., particularly from May to August. The simulated MCS precipitation shows little diurnal variation and lasts too long, while MCS precipitation area is too large and intensity is too weak. The model is able to simulate four types of observed FLSMP associated with frontal systems and low-level jets (LLJ) in spring, but the frequencies are underestimated because of low-level dry bias and weaker LLJ. Precipitation simulated under different FLSMPs peak during daytime, in contrast to the observed nocturnal peak. Implications of these findings for future model development and diagnostics are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Study of Pile-Up Effects in Decay Energy Spectroscopy

Pile-up is an unavoidable complication for cryogenic detectors with relatively large heat capacities and slow rise time, such as systems for decay energy spectroscopy employing large Au absorbers. We have simulated the spectral response of such slow cryogenic detectors using Monte-Carlo algorithms to investigate the effects of pile-up on absolute and relative activity measurements. We focus on the impact of non-distinguishable pile-up that occurs when the rising edges of two waveforms originating from different events overlap in time and are interpreted as a single event. This effect can not be readily identified and corrected in experimental data. We investigated two representative cases of absolute decay counting and plutonium isotopic analysis and find that pile-up can distort the reconstruction of both the absolute and relative activities. This Monte-Carlo methodology quantifies of pile-up effects and provides a systematic methodology for calculating corrective factors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗