SEARCH · Engineering Papers

Results for “Simplified Simulator”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Impact of fusion reactor neutronics modeling for transmutation and thermal feedback *

Abstract Fusion neutronics calculations provide important metrics pertinent to fusion device operations, such as tritium breeding ratios (TBRs) and data on heat deposition, material activation, and damage. Because of the high computational burden required to generate a high-fidelity Monte Carlo simulation of a 3D fusion device, various assumptions are made to reduce computational time by simplifying the reactor model or the calculation iteration. This paper explores the impact of fusion neutronics metrics such as the TBR and decay heat of structural materials based on assumptions of material composition in the fusion reactor and temperature modeling of materials. Results show that for compact tokamaks with high power and long operational cycles, the transmutation of structural materials is significant enough to cause a substantial change in the flux spectrum and decrease the TBR by 1.68% after 2 years of full power operation. Additionally, assuming a constant temperature and material density can impact the TBR calculations up to 3%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Nonlinear quantum behavior of ultrashort-pulse optical parametric oscillators

The quantum features of ultrashort-pulse optical parametric oscillators (OPOs) are investigated theoretically in the nonlinear regime near and above threshold. Viewing the pulsed OPO as a multimode open quantum system, we rigorously derive a general input-output model that features nonlinear coupling among many cavity (i.e., system) signal modes and a broadband single-pass (i.e., reservoir) pump field. Under appropriate assumptions, our model produces a Lindblad master equation with multimode nonlinear Lindblad operators describing two-photon dissipation and a multimode four-wave-mixing Hamiltonian describing a broadband, dispersive optical cascade, which we show is required to preserve causality. To simplify the multimode complexity of the model, we employ a supermode decomposition to perform numerical simulations in the regime where the pulsed supermodes experience strong single-photon nonlinearity. We find that the quantum nonlinear dynamics induces pump depletion as well as corrections to the below-threshold squeezing spectrum predicted by linearized models. We also observe the formation of non-Gaussian states with Wigner-function negativity and show that the multimode interactions with the pump, both dissipative and dispersive, can act as effective decoherence channels. Finally, we briefly discuss some experimental considerations for potentially observing such quantum nonlinear phenomena with ultrashort-pulse OPOs on nonlinear nanophotonic platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Microgrid Modeling for Stability Analysis

Here this document is a summary of a report prepared by the IEEE PES Task Force (TF) on Microgrid (MG) Dynamic Modeling, IEEE Power and Energy Society, Tech. Rep. PES-TR106, 2023. In this paper, the major issues and challenges in microgrid modeling for stability analysis are discussed, and a review of state-of-the-art modeling approaches and trends is presented. In the context of the IEEE 1547 standard, the document covers issues associated with component models for MG dynamic studies and simulations, including generator and grid modeling, full and average converter models, unbalanced and balanced system conditions, dynamic and static loads, protection requirements, and detailed and simplified controls considering communications delays, packet losses, and security issues. Considering the future integration of grids and MGs to form broad integrated networks, a discussion is presented of the use of phasor vis-à-vis electromagnetic transient simulation tools for MG dynamic stability studies, as well as modeling scale-up issues and MG equivalent models. Specifically white-, grey-, and black-box models, are presented. This TF paper and companion report constitute a modeling guide for R&D groups working on developments and standards of MGs with a focus on stability issues.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

A High Resolution Simulation of a Single Shock-Accelerated Particle

We report particle drag models, which capture macroviscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a two-dimensional axisymmetric Navier–Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted using the hydrocode FLAG. FLAGs capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated—a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAGs Navier–Stokes and heat diffusion solutions and to provide a rationale for recent experimental particle drag measurements.

42 ENGINEERING↗

DOE OSTI2020

A Numerical Simulation of a Single Shock-Accelerated Particle

Particle drag models, which capture macro viscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a 2D axisymmetric Navier-Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted in the hydrocode FLAG. FLAG’s capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated – a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAG’s Navier-Stokes and heat diffusion solutions, and to provide a rationale for recent experimental particle drag measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Using Pilot Jobs and CernVM File System for Simplified Use of Containers and Software Distribution

High Energy Physics (HEP) experiments entail an abundance of computing resources, i.e. sites, to run simulations and analyses by processing data. This requirement is fulfilled by local batch farms, grid sites, private/commercial clouds, and supercomputing centers via High Throughput Computing (HTC). The growing needs of such experiments and resources being prone to trends of heterogeneity make it difficult for physicists to handle these resources directly. Additionally, HEP collaborations heavily rely on data and software releases, typically in the order of tens of gigabytes, while conducting simulations and analyses. Hence, aspects of scalability, reliability, and maintenance become crucial with regards to the distribution of the necessary data and software stack. The GlideinWMS [4] framework helps with the resource management problem by using pilot jobs, aka Glideins, to provision reliable elastic virtual clusters. Glideins are submitted to unreliable heterogeneous resources which are validated and customized by the Glideins to make the worker nodes available for end-user job execution. On the other hand, the CernVM File System (CernVM-FS or CVMFS) [1] helps with data distribution. It is a write-once, read-everywhere filesystem used to deploy scientific software to thousands of nodes on a worldwide distributed computing infrastructure. CVMFS is based on the Hyper Text Transfer Protocol and has been widely used within the particle physics community for (1) distributing experiment software and data such as calibrations, and (2) facilitating containerization by efficiently hosting container images along with providing containerization software, especially Singularity [3] GlideinWMS relies on CVMFS installed locally on the computing resources to satisfy the experiments' software needs. This requires system administrators' effort to install and maintain CVMFS at the sites and limits the use of sites, especially HPC resources, that do not have CVMFS installed. This poster presents a solution, taking advantage of Glideins to provide CVMFS at most sites without the need for a local installation. Doing so expands the pool of resources available for HEP experiments and reduces the effort of system administrators for current resources. Additionally, the proposed solution allows GlideinWMS to also start Singularity [3], a containerization software that can run unprivileged, on sites where neither CVMFS nor Singularity are available, including HPC sites. The benefits provided by this solution are: (1) lower overhead for site administrators in that they have less software to install, (2) an expanded pool of resources that run user jobs with easy access to software and data provided by CVMFS, thus making life easier for the scientists, and (3) improved flexibility to use HPC resources by enabling GlideinWMS pilot jobs to support HPC sites.

Urs, Namratha↗

DOE OSTI2023

Realistic fuel spray modeling for gasoline direct injection engine applications

Fuel spray modeling plays a critical role during modern gasoline direct injection (GDI) engine development due to fuel injection’s dominant impact on engine performance and emissions as well as the complex physical processes involved. In engineering three-dimensional (3D) computational fluid dynamics (CFD) simulations, the liquid-phase fuel atomization, evaporation, and mixing are usually modeled with the discrete droplet model (DDM) adopting a Lagrangian approach for multiphase CFD simulations. To this end, general practices heavily depend on the reduced order characterization of the injector nozzle flow. However, such simplified injector modeling may lead to insufficient representations of the complex spray dynamics. To tackle this problem, this study proposes a novel workflow to numerically evaluate GDI sub-cooled and flash-boiling sprays under engine-relevant conditions using a side-mounted GDI injector together with real gasoline fuel properties. The workflow introduces a one-way coupling (OWC) method leveraging high-fidelity nozzle flow simulations to provide realistic boundary conditions to the Lagrangian injector model. The proposed workflow was first verified in a constant volume chamber (CVC) environment and then implemented in a practical GDI engine setup to study spray morphology, fuel-air mixing, and wall-wetting propensity. In addition, detailed comparison was performed between the OWC method and the conventional rate of injection (ROI) routine. In conclusion, quantitative analysis of spray characteristics was conducted to highlight possible source of discrepancies of the conventional ROI method.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

thevenin: An Equivalent Circuit Modeling Package [SWR-24-132]

This package is a wrapper for the well-known Thevenin equivalent circuit model for simulating battery behavior with a focus on ease of use. The model is comprised of a single series resistor followed by any number of parallel RC pairs. The package includes an intuitive "experiment" interface, which simplifies the programming of constant or dynamic current, voltage, and/or power-driven loads, making it ideal for a wide range of applications, from basic tests to complex simulations. It also supports state-of-charge (SOC) and temperature-dependent properties, allowing the model to be more easily calibrated against real cells.

Randall, Corey↗

DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs: Preprint

Subsurface data analysis, reservoir modeling, and machine learning (ML) techniques have been applied to the Brady Hot Springs (BHS) geothermal field in Nevada, USA to further characterize the subsurface and assist with optimizing reservoir management. Hundreds of reservoir simulations have been conducted in TETRAD-G and CMG STARS to explore different injection and production fluid flow rates and allocations and to develop a training data set for ML. This process included simulating the historical injection and production since 1979 and prediction of future performance through 2040. ML networks were created and trained using TensorFlow based on multilayer perceptron (MLP), long short-term memory (LSTM), and convolutional neural network (CNN) architectures. These networks took as input selected flow rates, injection temperatures, and historical field operation data and produced estimates of future production temperatures. This approach was first successfully tested on a simplified single fracture doublet system, followed by the application to the BHS reservoir. Using an initial BHS dataset with 37 simulated scenarios, the trained and validated network predicted the production temperature for 6 production wells with the mean absolute percentage error of less than 8%. In a complementary analysis effort, the principal component analysis applied to 13 BHS geological parameters revealed that vertical fracture permeability shows the strongest correlation with fault density and fault intersection density. A new BHS reservoir model was developed considering the fault intersection density as proxy for permeability. This new reservoir model helps to explore under-exploited zones in the reservoir. A data gathering plan to obtain additional subsurface data was developed; it includes temperature surveying for three idle injection wells, at which the reservoir simulations indicate high bottom-hole temperatures. The collected data assist with calibrating the reservoir model and may lead to converting these wells to producers to access under-exploited zones in the reservoir. Data gathering activities are planned for the first quarter of 2021.

40 EE - Geothermal Technologies Office (EE-4G)↗

DOE OSTI2020

Atomistic insights into the synergistic effects of tensile strain and hydroxyl group on increasing the thermal conductivity of monohydric alcohols as latent heat storage materials

The low thermal conductivity of linear monohydric alcohols significantly limits their performance when used as latent heat storage materials. It was found that imposing certain mechanical strains on nanomaterials with highly ordered structures can lead to appreciable thermal conductivity enhancement. To explore the potential of enhancing the intrinsic thermal conductivity of monohydric alcohols, a simplified ideal crystal model was employed for molecular arrangement in this work. The variation of atomistic heat transfer as a function of tensile strains was exploited by non-equilibrium molecular dynamics simulations. Furthermore, the influence of stretching on the vibration density of states and the corresponding molecular morphology were revealed. The results suggested that the thermal conductivity of monohydric alcohols with such simplified ideal crystalline structures can increase by about 170 % at the strain of 0.1. The distance between neighboring interfaces was found to determine the energy of inter-molecular interactions at the crystal grain boundaries. The hydrogen bonds energy has not reached the critical value when increasing tensile strains from 0 to 0.1. In contrast, the strongest van der Walls interactions for alkanes were obtained at the strain of 0.06 where the interfacial distances of alkanes have achieved their equilibrium distances. In addition, the interfacial heat transfer coefficient of ideal crystal monohydric alcohols at hydroxyl interfaces was found to be threefold than that between methyl interfaces, in consistent with the strength of inter-molecular interactions between functional groups at the interfaces.

36 MATERIALS SCIENCE↗

DOE OSTI2020

A Simplified Building Modeling Approach for Identifying Whole Building Retrofits

Building Energy Modeling (BEM) is an effective strategy for optimizing new and retrofit building design, evaluating its energy savings potential and rating its energy performance. However, because of the complexity and expense of modeling, only a small percentage of buildings are simulated. Over 80% of buildings are <25,000 ft2, and the cost associated with BEM can be a big deterrent for reaching this subsector, which usually ends up complying with energy codes and qualifying for incentives using prescriptive measures or generalized design guides. Similarly, efficiency considerations for small building retrofits are often identified through deemed measures or Technical Reference Manuals, which though a scalable and convenient approach, can limit innovation in design and optimization of measures. Simplification of the modeling process, where appropriate, can lead to increased energy savings and more informed decision making. This paper will discuss the development of a ruleset for a simplified approach based on the Performance Rating Method (PRM) contained in ASHRAE Standard 90.1. The approach aims to simplify and lower the cost of the energy modeling process in a controlled and documented manner. An anticipated outcome is an increased use of BEM and its effectiveness in the design, retrofit and operation of commercial buildings. In addition, simplification of the modeling process will result in fewer errors from misinterpretations of both program requirements and program intent in simulation software. The paper will discuss the need and advantages of a simplified modeling approach for whole building retrofits and the technical process for achieving the same.

Goel, Supriya↗

DOE OSTI2020

Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling

Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Decoding the formation of hammerhead ion populations observed by Parker Solar Probe

Context. In situ observations by the Parker Solar Probe (PSP) have revealed new properties of the proton velocity distributions (VDs), including hammerhead features that suggest a non-isotropic broadening of the beams. Aims. The present work proposes a very plausible explanation for the formation of hammerhead proton populations through the action of a proton firehose-like instability triggered by the proton beam. Methods. We investigated a self-generated firehose-like instability driven by the relative drift of ion populations using a simplified moment-based quasi-linear (QL) theory. While simpler and faster than advanced numerical simulations, this toy model provided rapid insights and concisely highlighted the role of plasma micro-instabilities in relaxing the observed anisotropies of particle VDs in the solar wind and space plasmas. Results. The QL theory proposed here shows that the resulting transverse waves are right-hand polarized and have two consequences on the protons: (i) They reduce the relative drift between the beam and the core, but above all, (ii) they induce a strong perpendicular temperature anisotropy specific to the observed hammerhead ion beam. Moreover, the long-run QL results suggest that these hammerhead distributions are rather transitory states that are still subject to relaxation mechanisms, in which instabilities such as the one discussed here are very likely involved. This foundational work motivates future detailed studies using advanced methods.

Shaaban, Shaaban M. (ORCID:000000030465598X)↗

DOE OSTI2021

Peculiar velocity estimation from kinetic SZ effect using deep neural networks

ABSTRACT The Sunyaev–Zel’dolvich (SZ) effect is expected to be instrumental in measuring velocities of distant clusters in near future telescope surveys. We simplify the calculation of peculiar velocities of galaxy clusters using deep learning frameworks trained on numerical simulations to avoid the independent estimation of the optical depth. Images of distorted photon backgrounds are generated for idealized observations using one of the largest cosmological hydrodynamical simulations, the Magneticum simulations. The model is tested to determine its ability of estimating peculiar velocities from future kinetic SZ observations under different noise conditions. The deep learning algorithm displays robustness in estimating peculiar velocities from kinetic SZ effect by an improvement in accuracy of about 17 per cent compared to the analytical approach.

Wang, Yuyu (ORCID:0000000202458547)↗

DOE OSTI2025

Evaluating the Impact of Tritium Permeation Membrane Performance and Direct Internal Recycling on Fusion Fuel Cycle Efficiency Using TMAP8

An efficient fuel cycle is vital to sustainable and cost-effective energy generation in fusion systems. Since tritium is not widely available, fusion systems must breed their own tritium for sustainable fusion deuterium-tritium reactions. An inefficient fuel cycle increases the tritium inventory needed for operations, which increases costs, constraints on tritium management systems, and safety concerns. A fuel cycle model is a powerful tool for understanding tritium inventories and flow rates across all systems in the fuel cycle. By simplifying the technical details into time-dependent tritium flow rates and inventories, the model can simulate the entire fuel cycle with high computational efficiency, even for technologies that are still under development. It can therefore quantify the impact of new tritium management technologies on fuel cycle efficiency. To evaluate the impact of key components on reducing tritium inventory, we are using and expanding an existing fuel cycle models based on latest advancements in fuel cycle research. The new model integrates Tritium Permeation Membrane (TPM) and Direct Internal Recycling (DIR) to enhance tritium transport from blanket breeders and plasma exhaust. These fuel cycle models are implemented in TMAP8 (Tritium Migration Analysis Program, version 8), a MOOSE-based open-source application designed to provide cutting-edge capabilities for tritium transport and fuel cycle modeling. The study aims to demonstrate the extensibility of existing fuel cycle modeling capability in TMAP8 and to offer a proof-of-principle design for future fusion plant systems. The presentation will cover the performance of fuel cycle modeling capabilities available in TMAP8, highlight advancements in fuel cycle research, and present a sensitivity analysis of these models. The results underline potential approaches and technology solutions to lower tritium inventory requirements, highlighting their role in shaping the future of fusion energy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Leveraging generative artificial intelligence to bridge domain gaps in wind turbine research

A central challenge in wind turbine health monitoring is the scarcity of real-world data due to limited instrumentation, leading researchers to rely on simulation models that often suffer from reduced fidelity. However, even within simulation environments, discrepancies arise because of modeling assumptions, and configuration fidelities, creating domain gaps that limit the transferability of learned representations. Here, to investigate domain translation under controlled conditions, this project explores the use of generative artificial intelligence, specifically cycle-consistent generative adversarial networks (CGANs), to bridge the gap between OpenFAST simulation models representing 1.5 MW and 5 MW wind turbines. A physics-informed CGAN architecture is introduced, where a simplified turbine tower dynamics model is incorporated into the training loss to ensure physically consistent outputs. Quantitative results showed moderate to high agreement in frequency-domain features. Incorporating the physics-informed loss function improved the R 2 values by 30%, reduced the RMSE from 1.39 to 1.1 m/s 2 , and reduced training time by 82%. Furthermore, under increased turbulence intensity (IEC Category A), the RMSE remained stable at approximately 1.1 m/s 2 . While the present study is entirely simulation-based, it establishes a pipeline for evaluating physics-informed generative domain translation, which may serve as a foundation for future simulation-to-reality validation studies.

17 WIND ENERGY↗

DOE OSTI2020

Deep learning interfacial momentum closures in coarse-mesh CFD two-phase flow simulation using validation data

Multiphase flow phenomena have been widely observed in the industrial applications while it remains a challenging yet unsolved problems. Three-dimensional computational fluid dynamics (CFD) approaches resolve the flow fields on a finer special and temporal scales which can complement the dedicated experimental study. However, closures have to be introduced to reflect the underlying physics in multiphase flow. Among them, the interfacial forces, including drag, lift, turbulent dispersion and wall lubrication forces, play in important role on the bubble’s distribution and migration in liquid-vapor two-phase flow. Development of those closures traditionally rely on the experimental data and analytical derivation with simplified assumptions which usually cannot deliver a universal solution across wide range of flow conditions. In this paper, a data-driven approach, named as Feature Similarity Measurement (FSM), is developed and applied to improve the simulation capability of two-phase flow with coarse-mesh CFD approach. Interfacial momentum transfer in adiabatic bubbly flow serves as the focus of the present study. Both a mature and a simplified set of interfacial closures are taken as the low fidelity data. Experimental data and fine mesh CFD simulations results are adopted as high-fidelity data. Qualitative and quantitative analysis are performed in this paper which reveals that FSM can substantially improve the prediction of coarse mesh CFD model regardless of the choice of interfacial closures and it provides scalability and consistency across discontinuous flow regimes. Furthermore, it demonstrates that data-driven method can aid the multiphase flow modeling by exploring the connections between local physical features and simulation errors.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Dissipative ground state preparation in ab initio electronic structure theory

Dissipative engineering is a powerful tool for quantum state preparation, and has drawn significant attention in quantum algorithms and quantum many-body physics in recent years. In this work, we introduce a novel approach using the Lindblad dynamics to efficiently prepare the ground state for general ab initio electronic structure problems on quantum computers, without variational parameters. These problems often involve Hamiltonians that lack geometric locality or sparsity structures, which we address by proposing two generic types of jump operators for the Lindblad dynamics. Type-I jump operators break the particle number symmetry and should be simulated in the Fock space. Type-II jump operators preserves the particle number symmetry and can be simulated more efficiently in the full configuration interaction space. For both types of jump operators, we prove that in a simplified Hartree-Fock framework, the spectral gap of our Lindbladian is lower bounded by a universal constant. For physical observables such as energy and reduced density matrices, the convergence rate of our Lindblad dynamics with Type-I jump operators remains universal, while the convergence rate with Type-II jump operators only depends on coarse grained information such as the number of orbitals and the number of electrons. To validate our approach, we employ a Monte Carlo trajectory-based algorithm for simulating the Lindblad dynamics for full ab initio Hamiltonians, demonstrating its effectiveness on molecular systems amenable to exact wavefunction treatment.

Quantum chemistry↗