Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation System”

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 973 records · Page 54

Demonstrating the Transient System Impact of Cyber-Physical Events Through Scalable Transmission and Distribution (T&D) Co-Simulation

Modern power systems become more vulnerable to cyber threats due to their growing interconnectivity, interdependence, and complexity. Widespread deployment of distributed energy resources (DERs) further expands the threat landscape to the grid edge, where fewer cybersecurity protections exist. In this article, a systematic cyber-physical events demonstration, enabled by an integrated transmission, distribution, and communication co-simulation framework, is performed. It analyzes cyber risks to power grid under DER-enabled automatic generation control from different angles. Unlike existing works, the simulation captures millisecond-to-minutes frequency and voltage transient dynamics at a cross-region system scale.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Quantum Simulation of an Open System: A Dissipative 1+1D Ising Model

The 1+1D Ising model is an ideal benchmark for quantum algorithms, as it is very well understood theoretically. This is true even when expanding the model to include complex coupling constants. In this work, we implement quantum algorithms designed for the simulation of open or complex coupling quantum field theories on IBM devices with a focus on the measurement of the Lee-Yang edge singularity. This feature corresponds (at large volumes) to a phase transition, and our successful reproduction of the transition represents a non-trivial test for current hardware and its ability to distinguish features of interest in quantum field theories.

Gustafson, Erik↗

Digital quantum computation of fermion-boson interacting systems

We introduce a method for representing the low-energy subspace of a bosonic field theory on the qubit space of digital quantum computers. This discretization leads to an exponentially precise description of the subspace of the continuous theory thanks to the Nyquist-Shannon sampling theorem. The method makes the implementation of quantum algorithms for purely bosonic systems as well as fermion-boson interacting systems feasible. We present algorithmic circuits for computing the time evolution of these systems. The complexity of the algorithms scales polynomially with the system size. The algorithm is a natural extension of the existing quantum algorithms for simulating fermion systems in quantum chemistry and condensed-matter physics to systems involving bosons and fermion-boson interactions and has a broad variety of potential applications in particle physics, condensed matter, etc. Due to the relatively small amount of additional resources required by the inclusion of bosons in our algorithm, the simulation of electron-phonon and similar systems can be placed in the same near-future reach as the simulation of interacting electron systems. We benchmark our algorithm by implementing it for a two-site Holstein polaron problem on an Atos Quantum Learning Machine quantum simulator. The polaron quantum simulations are in excellent agreement with the results obtained by exact diagonalization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploratory Nuclear Reactor Safety Analysis and Visualization via Integrated Topological and Geometric Techniques

A recent trend in the nuclear power engineering field is the implementation of heavily computational and time consuming algorithms and codes for both design and safety analysis. In particular, the new generation of system analysis codes aim to embrace several phenomena such as thermo-hydraulic, structural behavior, and system dynamics, as well as uncertainty quantification and sensitivity analyses. The use of dynamic probabilistic risk assessment (PRA) methodologies allows a systematic approach to uncertainty quantification. Dynamic methodologies in PRA account for possible coupling between triggered or stochastic events through explicit consideration of the time element in system evolution, often through the use of dynamic system models (simulators). They are usually needed when the system has more than one failure mode, control loops, and/or hardware/process/software/human interaction. Dynamic methodologies are also capable of modeling the consequences of epistemic and aleatory uncertainties. The Monte-Carlo (MC) and the Dynamic Event Tree (DET) approaches belong to this new class of dynamic PRA methodologies. The major challenges in using MC and DET methodologies (as well as other dynamic methodologies) are the heavier computational and memory requirements compared to the classical ET analysis. This is due to the fact that each branch generated can contain time evolutions of a large number of variables (about 50,000 data channels are typically present in RELAP) and a large number of scenarios can be generated from a single initiating event (possibly on the order of hundreds or even thousands). Such large amounts of information are usually very difficult to organize in order to identify the main trends in scenario evolutions and the main risk contributors for each initiating event. This report aims to improve Dynamic PRA methodologies by tackling the two challenges mentioned above using: 1) adaptive sampling techniques to reduce computational cost of the analysis and 2) topology-based methodologies to interactively visualize multidimensional data and extract risk-informed insights. Regarding item 1) we employ learning algorithms that aim to infer/predict simulation outcome and decide the coordinate in the input space of the next sample that maximize the amount of information that can be gained from it. Such methodologies can be used to both explore and exploit the input space. The later one is especially used for safety analysis scopes to focus samples along the limit surface, i.e. the boundaries in the input space between system failure and system success. Regarding item 2) we present a software tool that is designed to analyze multi-dimensional data. We model a large-scale nuclear simulation dataset as a high-dimensional scalar function defined over a discrete sample of the domain. First, we provide structural analysis of such a function at multiple scales and provide insight into the relationship between the input parameters and the output. Second, we enable exploratory analysis for users, where we help the users to differentiate features from noise through multi-scale analysis on an interactive platform, based on domain knowledge and data characterization. Our analysis is performed by exploiting the topological and geometric properties of the domain, building statistical models based on its topological segmentations and providing interactive visual interfaces to facilitate such explorations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Structure and Evolution of Mesoscale Convective Systems: Sensitivity to Cloud Microphysics in Convection‐Permitting Simulations Over the United States

Abstract Regional climate simulations over the continental United States were conducted for the 2011 warm season using the Weather Research and Forecasting model at convection‐permitting resolution (4 km) with two commonly used microphysics parameterizations (Thompson and Morrison). Sensitivities of the simulated mesoscale convective system (MCS) properties and feedbacks to large‐scale environments are systematically examined against high‐resolution geostationary satellite and 3‐D mosaic radar observations. MCS precipitation including precipitation amount, diurnal cycle, and distribution of hourly precipitation intensity are reasonably captured by the two simulations despite significant differences in their simulated MCS properties. In general, the Thompson simulation produces better agreement with observations for MCS upper level cloud shield and precipitation area, convective feature horizontal and vertical extents, and partitioning between convective and stratiform precipitation. More importantly, Thompson simulates more stratiform rainfall, which agrees better with observations and results in top‐heavier heating profiles from robust MCSs compared to Morrison. A stronger dynamical feedback to the large‐scale environment is therefore seen in Thompson, wherein an enhanced mesoscale vortex behind the MCS strengthens the synoptic‐scale trough and promotes advection of cool and dry air into the rear of the MCS region. The latter prolongs the MCS lifetimes in the Thompson relative to the Morrison simulations. Hence, different treatment of cloud microphysics not only alters MCS convective‐scale dynamics but also has significant impacts on their macrophysical properties such as lifetime and precipitation. As long‐lived MCSs produced 2–3 times the amount of rainfall compared to short‐lived ones, cloud microphysics parameterizations have profound impact in simulating extreme precipitation and the hydrologic cycle.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of Global Fire Simulations in CMIP6 Earth System Models

Fire is the primary form of terrestrial ecosystem disturbance on a global scale and an important Earth system process. Most Earth system models (ESMs) have incorporated fire modeling, with 19 of them submitting model outputs of fire-related variables to the Coupled Model Intercomparison Project Phase 6 (CMIP6). This study provides the first comprehensive evaluation of CMIP6 historical fire simulations by comparing them with multiple satellite-based products and charcoal-based historical reconstructions. Our results show that most CMIP6 models simulate the present-day global burned area and fire carbon emissions within the range of satellite-based products. They also capture the major features of observed spatial patterns and seasonal cycles, the relationship of fires with precipitation and population density, and the influence of the El Niño–Southern Oscillation (ENSO) on the interannual variability of tropical fires. Regional fire carbon emissions simulated by the CMIP6 models from 1850 to 2010 generally align with the charcoal-based reconstructions, although there are regional mismatches, such as in southern South America and eastern temperate North America prior to the 1910s and in temperate North America, eastern boreal North America, Europe, and boreal Asia since the 1980s. The CMIP6 simulations have addressed three critical issues identified in CMIP5: (1) the simulated global burned area being less than half of that of the observations, (2) the failure to reproduce the high burned area fraction observed in Africa, and (3) the weak fire seasonal variability. Furthermore, the CMIP6 models exhibit improved accuracy in capturing the observed relationship between fires and both climatic and socioeconomic drivers and better align with the historical long-term trends indicated by charcoal-based reconstructions in most regions worldwide. However, the CMIP6 models still fail to reproduce the decline in global burned area and fire carbon emissions observed over the past 2 decades, mainly attributed to an underestimation of anthropogenic fire suppression, and the spring peak in fires in the Northern Hemisphere midlatitudes, mainly due to an underestimation of crop fires. In addition, the model underestimates the fire sensitivity to wet–dry conditions, indicating the need to improve fuel wet-ness estimation. Based on these findings, we present specific guidance for fire scheme development and suggest a postprocessing methodology for using CMIP6 multi-model outputs to generate reliable fire projection products.

Wildfire, Earth system models↗

Application of Artificial Neural Networks to Reliable Nuclear Data for Nonproliferation Modeling and Simulation

Detection and identification of special nuclear materials (SNMs) are an essential part of the US nonproliferation effort. Modern cutting-edge SNM detection methodologies rely more and more on modeling and simulation techniques. Experiments with radiological samples in realistic configurations, is the ultimate tool that establishes the minimum detection limits of SNMs in a host of different geometries. Modern modeling and simulation approaches have the potential to significantly reduce the number of experiments with radioactive sources needed to determine these detection limits and reduce the financial barrier to SNM detection. Unreliable nuclear data is one of the principal causes of uncertainty in modeling and simulating nuclear systems. In particular, nuclear cross sections introduce a significant uncertainty in the nuclear data. The goal of this research is to develop a methodology that will autonomously extract the correct nuclear resonance characteristics of experimental data in a reliable way, a task previously left to expert judgement. Accurate nuclear data will in turn allow contemporary modeling and simulation to become far more reliable, de-escalating the extent of experimental testing. Modeling and simulation techniques reduce the use and distribution of radiological sources, while at the same time increase the reliability of the currently used methods for the detection and identification of SNMs.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

A Perspective on Quantum Computing Applications in Quantum Chemistry Using 25-100 Logical Qubits

The intersection of quantum computing and quantum chemistry represents a promising frontier for achieving quantum utility in domains of both scientific and societal relevance. Owing to the exponential growth of classical resource requirements for simulating quantum systems, quantum chemistry has long been recognized as a natural candidate for quantum computation. This perspective focuses on identifying scientifically meaningful use cases where early fault-tolerant quantum computers, which are considered to be equipped with approximately 25-100 logical qubits, could deliver tangible impact. While recent advances in classical computing have pushed the boundaries of tractable simulations to unprecedented scales, this logical-qubit regime represents the first window where quantum devices can pursue qualitatively distinct strategies, such as polynomial-scaling phase estimation, direct simulation of quantum dynamics, and active-space embedding, that remain challenging for classical solvers, such as multireference charge-transfer and conical-intersection states central to photochemistry and materials design. We highlight near-term opportunities in algorithm and software design, discuss representative chemical problems suited for quantum acceleration, and propose strategic roadmaps and collaborative pathways for advancing practical quantum utility in quantum chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of an Advanced Stimulation / Production Predictive Simulator for Enhanced Geothermal Systems

There are several well-known obstacles to the successful deployment of EGS projects on a commercial scale, of course. EGS projects are expected to be deeper, on the average, than conventional “natural” geothermal reservoirs, and drilling costs are already a formidable barrier to conventional geothermal projects. Unlike conventional resources (which frequently announce their presence with natural manifestations such as geysers, hot springs and fumaroles), EGS prospects are likely to appear fairly undistinguished from the earth surface. And, of course, the probable necessity of fabricating a subterranean fluid circulation network to mine the heat from the rock (instead of simply relying on natural, pre-existing permeable fractures) adds a significant degree of uncertainty to the prospects for success. Accordingly, the basic motivation for the work presented herein was to try to develop a new set of tools that would be more suitable for this purpose. Several years ago, the Department of Energy’s Geothermal Technologies Office recognized this need and funded a cost-shared grant to our company (then SAIC, now Leidos) to partner with Geowatt AG of Zurich, Switzerland and undertake the development of a new reservoir simulator that would be more suitable for EGS forecasting than the existing tools. That project has now been completed and a new numerical geothermal reservoir simulator has been developed. It is named “HeatEx” (for “Heat Extraction”) and is almost completely new, although its methodology owes a great deal to other previous geothermal software development efforts, including Geowatt’s “HEX-S” code, the STAR and SPFRAC simulators developed here at SAIC/Leidos, the MINC approach originally developed at LBNL, and tracer analysis software originally formulated at INEL. Furthermore, the development effort was led by engineers with many years of experience in using reservoir simulation software to make meaningful forecasts for real geothermal projects, not just software designers. It is hoped that, as a result, HeatEx will prove useful during the early stages of the development of EGS technology. The basic objective was to design a tool that could use field data that are likely to become available during the early phases of an EGS project (that is, during initial reconnaissance and fracture stimulation operations) to guide forecasts of the longer-term behavior of the system during production and heat-mining.

15 GEOTHERMAL ENERGY↗

Field test and simulation evaluation of variable refrigerant flow systems performance

Our study aims to compare the performance of a Variable Refrigerant Flow (VRF) system with a Roof Top Unit, (RTU) and a variable-air-volume (VAV) system through field tests and energy simulations. The field test was conducted in as similar conditions as possible between the two systems, such as the installation and operating environments of heating, the ventilation and air conditioning (HVAC) system, including internal heat gain and outdoor conditions, including buildings to compare the performance of the two systems accurately. A VRF system and RTU were installed at the test building located in Oak Ridge, Tennessee, in the USA. The same internal heat gain was generated at the same operating time of the two systems using lighting, electric heaters, and humidifiers inside the building. The HVAC system was alternately operated between cooling and heating operations to acquire energy performance data and to compare energy usage. Furthermore, an hourly building energy simulation model was developed with regard to the VRF system and RTU, and then the model was calibrated using actual measured data. Then, annual energy consumption of the two systems were compared and analyzed using the calibrated model. Moreover, additional analysis was conducted when the controlled discharge air temperature in the RTU was changed. The field test result showed that when energy consumptions of two systems were compared at the same outdoor conditions, using the weather-normalized model, the VRF system exhibited an energy reduction of approximately 17% during cooling operation and of approximately 74% during heating operations. A comparison on the annual energy consumption using simulations showed that the VRF system reduced energy consumption more than that of the RTU by 60%. Furthermore, when the discharge air temperature in the RTU was controlled according to the outdoor air temperature, energy consumption of the RTU was reduced by 6% in cooling operations and by 18% in heating operation. As a result, energy consumption of the VRF system was reduced by more than that of the RTU by 55% approximately.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Engineering Ralstonia eutropha for Production of Isobutanol (IBT) Motor Fuel from Carbon Dioxide, Hydrogen, and Oxygen Project Final Report

This research project is a collaboration between the Sinskey laboratory at MIT and the Worden laboratory at Michigan State University. The goal of the project is to produce Isobutanol (IBT), a branched-chain alcohol that can serve as a drop-in transportation fuel, through the engineered microbial biosynthesis of Carbon Dioxide, Hydrogen, and Oxygen using a novel bioreactor. This final technical report presents the findings of both the biological engineering work at MIT that extended the native branched-chain amino acid pathway of the wild type Ralstonia eutropha H16 to perform this biosynthesis, as well as the unique design, modeling, and construction of a bioreactor for incompatible gasses at Michigan State that enabled the operational testing of the complete system. This 105 page technical report summarizing the three years of research includes 72 figures and 11 tables of findings. Ralstonia eutropha (also known as Cupriavidus necator) is a Gram-negative, facultatively chemolithoautotrophic bacteria. It has been the principle organism used for the study of polyhydroxybutyrate (PHB) polymer biosynthesis. The wild-type Ralstonia eutropha H16 produces PHB as an intracellular carbon storage material while under nutrient stress in the presence of excess carbon. Under this stress, it can accumulate approximately 80 % of its cell dry weight (CDW) as this intracellular polymer. With the restoration of the required nutrients, the cells are then able to catabolize this polymer. If extracted from the cell, this PHB polymer can be processed into biodegradable and biocompatible plastics, however for this research, it is the efficient metabolic pathway channeling the captured carbon that is of interest. R. eutropha is further unique in that it contains two carbon-fixation Calvin–Benson–Bassham cycle operons, two oxygen-tolerant hydrogenases, and several formate dehydrogenases. It has also been much studied for its ability in the presence of oxygen, to fix carbon dioxide into complex cellular molecules using the energy from hydrogen. In this research project, engineered strains of R. eutropha redirected the excess carbon from PHB storage into the production of isobutanol and 3-methyl-1-butanol (branched-chain higher alcohols). These branched-chain higher alcohols can be used directly as substitutes for fossil-based fuels and are seen as alternative biofuels to ethanol and biodiesel. Importantly, these alcohols have approximately 98 % of the energy content of gasoline, 17 % higher than the current gasoline additive ethanol, without impacting corn market production for feed or food. Unlike ethanol, these branched-chain alcohols have low vapor pressure, hygroscopicity, and water solubility, which make them readily compatible with the existing pipelines, gasoline pumps, and engines in our transportation infrastructure. While the use of alternative energies from solar, wind, geothermal, and hydroelectric has spread for stationary power applications, these energy sources cannot be effectively or efficiently employed in current or future transportation systems. With the ongoing concerns of fossil fuel availability and price stability over the long term, alternative biofuels like branched-chain higher alcohols hold promise as a suitable transportation fuel in the future. We showed in our research that various mutant strains of R. eutropha with isobutyraldehyde dehydrogenase activity, in combination with the overexpression of plasmid-borne, native branched-chain amino acid biosynthesis pathway genes and the overexpression of heterologous ketoisovalerate decarboxylase gene, would produce isobutanol and 3-methyl-1-butanol when initiated during nitrogen or phosphorus limitation. Early on, we isolated one mutant R. eutropha strain which produced over 180 mg/L branched-chain alcohols in flask culture while being more tolerant of isobutanol toxicity. After the targeted elimination of genes encoding several potential carbon sinks (ilvE, bkdAB, and aceE), the production titer of the improved to 270 mg/L isobutanol and 40 mg/L 3-methyl-1-butanol. Semicontinuous flask cultivation supplied the cells with sufficient nutrients while minimizing the toxicity caused by isobutanol. Under this cultivation, the R. eutropha mutant grew and produced more than 14 g/L branched-chain alcohols over the duration of 50 days. These results demonstrate that R. eutropha carbon flux can be redirected from PHB to branched-chain alcohols and that engineered R. eutropha can be cultivated over prolonged periods of time for product biosynthesis. While this bioengineering work was being done at the Sinskey laboratory at MIT, the researchers at the Worden laboratory at Michigan State were working on the design and construction of the required specialty bioreactor for incompatible gasses (BIG) that would allow the safe feeding of microbes on Carbon Dioxide, Hydrogen, and Oxygen without explosive results. The early design and assembly work in year 1 incorporated a novel microbubble generator to maximize the bioavailability of gasses within the system comprised of small scale hollow fiber reactors. The early success of the microbubble generator eliminated the need to investigate potentially toxic surfactants within the system. For operational control, the system design incorporated a Opto22-based control network. The researchers also selected the specific hollow fiber material suitable for the bioreactor application. A variety of commercially available hollow fiber membranes were compared with regard to their pore sizes, cell affinity, and potential interference to cell viability assays. The selected membrane with its spongy layer was then tested for diffusivity of O2 and CO2. The instrumented system was then fully assembled for experimentally measuring the heterotrophic growth rate of immobilized R. eutropha cells. The requisite procedures for inoculation, measurement, and cleaning were established enabling the system performance to be validated under controlled laboratory conditions. In year 2, the researchers completed the Opto22 based cross-platform control network, and the system’s communications across the Sartorius fermentation system and Bruker gas chromatograph was established via open platform communications (OPC) protocol. Using the revised system, measurements were taken of the R.eutropha cell growth rate and substrate mass transfer rate in the hollow fiber membrane. Several IBT recovery strategies were explored and resin adsorption was determined to be optimal solution for lab scale operations. The adsorption capacity of the resin column was then measured and IBT desorption using methanol has been demonstrated. With the growing body of experimental data in hand, mathematical models were constructed to demonstrate and map the cellular kinetics, mass transfer of heterotrophic and autotrophic substrates in the hollow fiber, and the adsorption process in the resin column. A structured kinetic model was constructed to describe the competition between cell mass generation and IBT production. The reactor was then scaled up from single fiber to a membrane area of 180 cm2 and then further to 1 ft2. In Year 3 of the research, the IBT mass transfer across the membrane was characterized within the system with experiments to empirically measure the IBT diffusion coefficient in the BIG spongy layer. Using the refined mathematical models, the researchers are now able to explain the experimental observations and predict bioreactor performance under a wide range of experimental conditions. The Big system is able to demonstrate continuous controlled operations with the integrated IBT recovery system. Both heterotrophic and autotrophic production have been shown during continuous operation with heterotrophic and autotrophic stages. Performance of BIG system has been measured during continuous run with alternating heterotrophic growth on fructose and autotrophic product formation on H2, CO2, and O2. Volumetric productivities of IBT at 325 mg/(L day) and of 3M1B at 50 mg/(L day) were achieved, which were comparable to that achieved under heterotrophic conditions. Using the mathematical model, researchers are able to predict system performance for scaled-up BIG system. The apparent diffusion coefficient of IBT in the spongy layer of XM-50 hollow fiber membranes has been measured at various lumen liquid flow rates. The experiment is simulated in COMSOL to validate the results. The constructed COMSOL model is able to simulate BIG system performance in both batch and continuous mode. Mathematical simulations of the system performance have been run to identify the most crucial operational conditions, identifying the rate-limiting factors in autotrophic production of IBT, and quantitating the rate of IBT catabolism. Investigations of the productivity of the production system have suggested and the modeling of the system has revealed a particular sensitivity to the catabolism of the produced IBT by the engineered R. eutropha. Experiments have been designed and executed to quantify the IBT catabolism of R. eutropha, which open up possibilities for further system improvements through future, targeted bioengineering of the strain. Finally, the researchers at Michigan State performed an economic analysis of the system, based on the collective results, and their findings are presented in full within the report.

08 HYDROGEN↗

Incorporating Residential Smart Electric Vehicle Charging in Home Energy Management Systems

Electric vehicles (EVs) are expected to drastically increase residential electricity consumption and could provide a significant source of flexible demand. Aggregating smart EV charge controllers with other smart home devices through a home energy management system can lead to more optimal outcomes that benefit homeowners, utilities, and grid operators. Control strategies should consider occupant convenience by accounting for the need for fully charged EVs near the EV departure time. In this paper, we develop an EV charging framework that accounts for occupant convenience using OCHRE, a residential energy model, and foresee, a home energy management system. We simulate a community with high EV penetration and show that integrated, smart EV charging reduces peak demand and smooths night-time energy consumption. Simulation results show that the proposed control strategy nearly eliminates peak period EV charging and reduces the daily peak demand from EVs by 23%.

ADVANCED PROPULSION SYSTEMS,ENERGY CONSERVATION, C↗

Massively Parallel Simulation of Large Molecular Systems with Long-Range Interactions Final Report CRADA No. TC-0297-92-B

The purpose of this CRADA project is to develop and implement a large-scale molecular dynamics model for floppy molecules with long-ranged interaction on massively parallel supercomputers and accelerate the use of parallel supercomputers within the biomolecular and polymer industries. Within three years of the project, the partners will: create a new molecular dynamics computer simulation code (LAMMPS) from existing codes at SAndia National Laboratory-Albuquerque (SNL-A) and Lawrence Livermore National Laboratory (LLNL) and new algorithms and code developed during the CRADA project; design the resulting code to optimally use the emerging generation of massively parallel supercomputers; and apply the code to the solution of biomolecular problem suggested by Bristol-Myers Squibb (BMS) and polymer problems suggested by DuPont.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spatiotemporal quenches for efficient critical ground state preparation in the two-dimensional transverse field Ising model

Quantum simulators have the potential to shed light on the study of quantum many-body systems and materials, offering unique insights into various quantum phenomena. Although adiabatic evolution has been conventionally employed for state preparation, it faces challenges when the system evolves too quickly or the coherence time is limited. In such cases, shortcuts to adiabaticity, such as spatiotemporal quenches, provide a promising alternative. This paper numerically investigates the application of spatiotemporal quenches in the two-dimensional transverse field Ising model with ferromagnetic interactions, focusing on the emergence of the ground state and its correlation properties at criticality when the gap vanishes. We demonstrate the effectiveness of these quenches in rapidly preparing ground states in critical systems. Our simulations reveal the existence of an optimal quench front velocity at the emergent speed of light, leading to minimal excitation energy density and correlation lengths of the order of finite system sizes we can simulate. These findings emphasize the potential of spatiotemporal quenches for efficient ground state preparation in quantum systems, with implications for the exploration of strongly correlated phases and programmable quantum computing.

2-dimensional systems↗

iTOUGH2: A multiphysics simulation-optimization framework for analyzing subsurface systems

We report iTOUGH2 is a simulation-optimization framework for the TOUGH suite of nonisothermal multiphase flow models and related simulators of geophysical, geochemical, and geomechanical processes. After appropriate parameterization of subsurface structures and their properties, iTOUGH2 runs simulations for multiple parameter sets and analyzes the resulting output for parameter estimation through automatic model calibration, local and global sensitivity analyses, data-worth analyses, and uncertainty propagation analyses. Development of iTOUGH2 is driven by scientific challenges and user needs, with new capabilities continually added to both the forward simulator and the optimization framework. This paper provides a summary description of methods and features implemented in iTOUGH2, and discusses the usefulness and limitations of an integrated simulation-optimization workflow in support of the characterization and analysis of complex multiphysics subsurface systems.

58 GEOSCIENCES↗