Search NASA⌕ Search

SEARCH · Search NASA

Results for “Generalized predictive control”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 181 records · Page 10

Data Driven Commercial Building Energy Code Compliance and Technology Inventory for New York City

Building Performance Standards (BPS) are gaining national traction. A BPS will require new processes in the design, construction, and operation of buildings that take the occupants into account and enable predictive analysis to ensure compliance with current and future GHG emissions caps. In New York City, most buildings over 25,000 square feet will be regulated by a BPS starting in 2024, regardless of whether it is new construction permitted under current energy codes or an existing building. This research is one of the first to begin the evaluation of a long-term series of building policies in the context of an open data ecosystem, in cooperation with city agencies. Existing building policies enacted in NYC have ranged from building energy benchmarking and labeling to energy audits to the regulation of GHG emission in buildings. Through the development of a dataset related to building technologies and energy consumption, this project can help to evaluate if meaningful conclusions can be drawn for the data that has been largely self-reported in compliance with city regulations. This project will also provide lessons learned from a deep dive into these types of datasets to provide best practices for municipalities or states seeking to embark on policies like those enacted in NYC. In addition, a Building Automation System (BAS) Stretch Standard of Care (SSOC) for owners, designers, and building operators will enable the measurement and predictive analysis of energy consumption and GHG emissions at the plant, system, or component level, in anticipation of regulated GHG limits on buildings based on energy use. The SSOC is expected to be suitable for use on a national level. The primary feature of an SSOC is a standardized format for a set of BAS points that can be used to control and to gather data from individual plants, systems, or components that are related to building energy consumption. This project examined how measurements compare to prescriptive or simulation-based energy code targets, finding little correlation between predictive 8760-hour energy modeling and actual energy consumption for a small sample (n=27) of buildings constructed after 2015. Other analysis found that, while large multifamily housing (MFH) buildings showed a general trend similar to predicted reductions in energy use from the implementation of model commercial energy codes, this trend was not evident in the office, K-12 school, and hotel use groups in NYC. No upward or downward trends in energy consumption were found when buildings were grouped by size. Energy audit data were analyzed and it appears that there is bias by audit company on measures recommended to clients. Further research should be performed to cross-analyze this with other attributes, such as building size, vintage, and number of stories. Analysis found that for 281 buildings that were permitted and completed after 2015 and had submitted benchmarking data in 2022, between 81% and 96% (by use group) were found to be in compliance with the 2024 to 2029 NYC BPS emission caps, and between 55% and 89% were in compliance with the 2030-2034 caps. This work is beneficial to the public in helping policymakers and building stakeholders better understand the wide-ranging implications of a BPS.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Modeling and Simulation of Inrush Currents in Harmonic Domain

Modeling and simulation capabilities are critical to the stability analysis and evaluation of power distribution systems, with respect to the emphasis on resiliency, microgrids, and distributed energy resources. In this paper, a computational method in the harmonic domain is proposed for the periodic steady-state analysis of the nonlinear inrush current phenomenon. The efficient inrush calculation facilitates the predictions of current amplitudes for the power system operation and control. To demonstrate the accuracy and efficiency, simulation results in the harmonic domain are compared with results from PSCAD in an electromagnetic timescale, as well as the authors’ previous works in the frequency-domain. Impacts of the settings of both offset flux and interested harmonic order are discussed. In addition, within the proposed harmonic-domain method, a general approach that utilizes the discrete Fourier transform to obtain the response of a nonlinear device from a stimulus represented in the frequency-domain is utilized. This method can also be extended to perform the transient analysis in future, using trapezoidal rule for the integration.

Xie, Jing↗

Machine Learning for Automated Weld Quality Monitoring and Control

Resistance Spot Welding (RSW) is a critical process in the automotive industry, valued for its cost-effectiveness, short cycle time, and robustness. However, achieving consistent high-quality joints remains challenging due to the complex interplay of various factors, like materials, processes, and manufacturing uncertainties, etc. Under the collaborative project between Oak Ridge National Laboratory (ORNL) and General Motors (GM), we have developed a robust and expansible machine learning (ML) framework aimed at enhancing quality control in RSW. By harnessing the power of machine learning, we have developed the ability to ensure every aspect of the welding process, from the initial process design stage to the final weld joint quality. The framework operates by analyzing a variety of data streams, including in-line process signals, process parameters, materials, and postprocessed weld joint data. Through this analysis, the models have been trained to detect deviations from optimal quality standards, leveraging their ability to identify signature data patterns and anomalies within in-line signals and construct complex correlations between these signals and weld quality parameters. Meanwhile, the machine learning framework is designed to adapt to a variety of materials, including high strength steels and aluminum alloys, etc. Its flexible architecture facilitates the incorporation of diverse data sources and features, enabling precise modeling and prediction across a broad range of material properties and weld quality variables. The expansible ML frameworks represent a promising transformation in weld quality monitoring and control, empowering industry to achieve high levels of efficiency, consistency, and reliability in manufacturing.

42 ENGINEERING↗

Enhancing the Operational Resilience of Advanced Reactors with Digital Twins by Recurrent Neural Networks

Because of a lack of operation data during abnormal and accident scenarios, along with the existence of uncertainty in the evaluation model for transient and accident analysis, the established abnormal and emergency operating procedures can be biased in characterizing the reactor states and ensuring operational resilience. To improve state awareness and ensure operational flexibility for minimizing effects on the system due to anomaly, digital twin (DT) technology is suggested to support operator's decision-making by effectively extracting and using knowledge of the current and future plant states from the knowledge base. To demonstrate DT's capability for recovering the complete states of reactors and for predicting the future reactor behaviors, this paper develops and assesses both the diagnosis and prognosis DTs in a nearly autonomous management and control system for an Experimental Breeder Reactor-II simulator during different loss-of-flow scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Distinguishing photo-induced oxygen attack on alkyl chain versus conjugated backbone for alkylthienyl-benzodithiophene (BDTT)-based push–pull polymers

Synthetic design has enabled increasing power conversion efficiency advances in organic photovoltaics (OPV). One continuing knowledge gap is detailed understanding of single-material chemical (photo)stability. Many considerations are based on prior OPV-relevant donor homopolymer systems rather than the next-generation push–pull architectures. Generally, energetic offsets between the lowest occupied molecular orbital of the donor and molecular oxygen are assumed to dictate kinetics of photo-induced charge transfer to a super oxide radical. Herein we determine the ambient-induced photo-degradation pathways, presented as proposed site-specific arrow-pushing mechanisms, for five donor polymers all containing the same push unit – alkylthienyl-substituted-benzodithiophene (BDTT) – but with chemically distinct pull units. The donor-only polymer films were subject to controlled photobleaching in air as an accelerated degradation approach, coupled with simultaneous monitoring of absorptance. Sample subsets were periodically removed and analyzed with X-ray photoelectron spectroscopy, to evaluate near-surface chemical composition and oxygen additions to hetero reporter atoms on either the BDTT push unit or distinct pull unit (i.e., sulphur and nitrogen). Here, this methodology allows us to distinguish between different mechanisms of bond cleavage and formation. Overall, we find that neither polymer redox properties nor individual push or pull unit stability are sufficient to predict photo-oxidative degradation of these polymers. Rather, there is a greater dependence on the susceptibility of unique structural groups within a polymeric system. Alkyl chain oxygen addition is generally the first attack site and direct sulfur oxidation on the conjugated backbone occurs after saturation of the alkyl chain initiation sites. This work provides a standard that could be used to evaluate relative photo-oxidative (in)stability for new OPV materials quickly – prior to time-consuming device optimization – and demonstrates an effective methodology for correlating optical degradation with chemical structure alterations via spectroscopic signatures to guide synthetic design.

14 SOLAR ENERGY↗

Mechanistic modeling of in vitro transcription incorporating effects of magnesium pyrophosphate crystallization

The in vitro transcription (IVT) reaction used in the production of messenger RNA vaccines and therapies remains poorly quantitatively understood. Mechanistic modeling of IVT could inform reaction design, scale-up, and control. In this work, we develop a mechanistic model of IVT to include nucleation and growth of magnesium pyrophosphate crystals and subsequent agglomeration of crystals and DNA. To help generalize this model to different constructs, a novel quantitative description is included for the rate of transcription as a function of target sequence length, DNA concentration, and T7 RNA polymerase concentration. The model explains previously unexplained trends in IVT data and quantitatively predicts the effect of adding the pyrophosphatase enzyme to the reaction system. The model is validated on additional literature data showing an ability to predict transcription rates as a function of RNA sequence length.

59 BASIC BIOLOGICAL SCIENCES↗

Transient fuel performance analysis for the preliminary fuel concept of general atomics fast modular reactor

This study investigates the transient fuel performance of General Atomics Fast Modular Reactor (GA-FMR) during accident scenarios, focusing on the behavior of its innovative fuel system that combines high-assay low enriched uranium dioxide (HALEUO2) fuel with SiGA® ceramic matrix composite silicon carbide cladding. The preliminary fuel design’s response was analyzed during reactivity-initiated accidents (RIA) and loss of coolant accidents (LOCA) using BISON fuel performance analysis code, which included both the diffusion enhanced and BISON-FASTGRASS coupled UO 2 models. The RIA analysis demonstrated that effective reactivity control reduced fuel temperature, though with transient fission gas release resulting in additional tensile stress state on the cladding. LOCA simulations revealed differing predictions between the two models: the BISON UO 2 model showed more transient fission gas release but minimal pellet expansion, while the BISON-FASTGRASS UO 2 model predicted less pronounced fission gas release but more fuel swelling and thermal expansion, potentially leading to pellet-cladding mechanical interaction. Here, these findings highlight critical areas for fuel design optimization and identify knowledge gaps requiring further experimental and computational investigation to advance GA-FMR fuel development.

Lee, Soon K. [Argonne National Laboratory (ANL), A↗

Pattern formation on free surfaces via plastic buckling and periodic folding

Spontaneous formation of surface patterns such as folds and creases is widely observed in engineered systems, such as thin films and machined surfaces, as well as in natural systems such as floral petals and leaves. A recognized route for their formation is via surface elastic instabilities. Here we demonstrate an alternate route for their formation —plastic buckling and folding—that occurs when ductile materials are subject to large-strain (>1) deformation. Using in situ imaging of simple-shear deformation of polypropylene, we delineate the mechanics underlying formation of folds and the resulting surface morphology. Here, the fold pattern shows a single wavelength that scales linearly with the size of the plastic zone. An analytical model predicts onset of the plastic buckling instability and characteristics of the folds. Besides obvious applications for controllably generating surface features over large areas, our results suggest a contributing mechanism for origin of surface roughness in solids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Simple Lagrangian Parcel Model for the Initiation of Summer-time Mesoscale Convective Systems over the Central United States

Mesoscale convective systems (MCSs) account for more than 50% of summertime precipitation over the central United States and have a significant impact on local weather and hydrologic cycle. It is hypothesized that the inadequate treatment of MCSs is responsible for the long-standing warm and dry bias over the central United States in coarse-resolution general circulation model (GCM) simulations. In particular, a better understanding of MCS initiation is still lacking. Here a single-column Lagrangian parcel model is first developed to simulate the basic features of a rising parcel. This simple model demonstrates the collective effects of boundary layer moistening and dynamical lifting in triggering convective initiation and reproduces successfully its early afternoon peak with surface equivalent potential temperature as a controlling factor. It also predicts that convection is harder to trigger in the future climate under global warming, consistent with the results from convection-permitting regional climate simulations. Then, a multicolumn model that includes an array of single-column models aligned in the east–west direction and incorporates idealized cold pool interaction mechanisms is developed. The multicolumn model captures readily the cold pool–induced upscale growth feature in MCS genesis from initially scattered convection that is organized into a mesoscale cluster in a few hours. It also highlights the crucial role of lifting effects due to cold pool collision and spreading, subsidence effect, and gust front propagation speed in controlling the final size of mesoscale clusters and cold pool regions. This simple model should be useful for understanding fundamental mechanisms of MCS initiation and providing guidance for improving MCS simulations in GCMs.

54 ENVIRONMENTAL SCIENCES↗

Measurements of the groomed and ungroomed jet angularities in pp collisions at $ \sqrt{s}$ = 5.02 TeV

The jet angularities are a class of jet substructure observables which characterize the angular and momentum distribution of particles within jets. These observables are sensitive to momentum scales ranging from perturbative hard scatterings to nonperturbative fragmentation into final-state hadrons. We report measurements of several groomed and ungroomed jet angularities in pp collisions at $\sqrt{s}$ = 5.02 TeV with the ALICE detector. Jets are reconstructed using charged particle tracks at midrapidity (|η| < 0.9). The anti-$k_T$ algorithm is used with jet resolution parameters R = 0.2 and R = 0.4 for several transverse momentum $p^{ch jet}_{T}$ intervals in the 20–100 GeV/c range. Using the jet grooming algorithm Soft Drop, the sensitivity to softer, wide-angle processes, as well as the underlying event, can be reduced in a way which is well-controlled in theoretical calculations. We report the ungroomed jet angularities, λ α , and groomed jet angularities, λ α,g , to investigate the interplay between perturbative and nonperturbative effects at low jet momenta. Various angular exponent parameters α = 1, 1.5, 2, and 3 are used to systematically vary the sensitivity of the observable to collinear and soft radiation. Results are compared to analytical predictions at next-to-leading-logarithmic accuracy, which provide a generally good description of the data in the perturbative regime but exhibit discrepancies in the nonperturbative regime. Moreover, these measurements serve as a baseline for future ones in heavy-ion collisions by providing new insight into the interplay between perturbative and nonperturbative effects in the angular and momentum substructure of jets. They supply crucial guidance on the selection of jet resolution parameter, jet transverse momentum, and angular scaling variable for jet quenching studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electron interference in atomic ionization by two crossing polarized ultrashort pulses

Formation of geometrically regular interference patterns in the photoelectron momentum distributions (PMDs) corresponding to the photoionization of atoms by two single-color, crossing ultrashort pulses is investigated both analytically and numerically. It is shown that, in contrast to the photoionization by monochromatic pulses, PMDs for the ionization by crossing and co-propagating broadband pulses are essentially different (unless both pulses are linearly polarized), namely, when one pulse is linearly polarized along the propagation direction, ˆk, of the circularly polarized (CP) pulse, then interference maxima (minima) of the ionization probability have the form of three-dimensional single-arm regular spirals which are wound along ˆk. Next, the interference maxima (minima) of the ionization probability by a pair of crossing elliptically polarized pulses have the form of either Newton's rings or two-arm Fermat's spirals, depending on the position of a detection plane. Remarkably, these regular patterns occur only for certain values of the pulse ellipticities, and they become distorted for CP pulses. For both above-mentioned pulse configurations, the features of interference patterns depend on the time delay between pulses, their relative electric field amplitude, and relative carrier-envelope phase. Our predictions, illustrated by the numerical results for the ionization of H and He atoms by two orthogonal pulses, are quite general and we expect them to be valid for the ionization of any randomly oriented atomic or molecular target.

74 ATOMIC AND MOLECULAR PHYSICS↗

Testing quantum gravity using pulsed optomechanical systems

An interesting idea, dating back to Feynman [Report from Chapel Hill Conference, edited by C. M. DeWitt and D. Rickles (1957)], argues that quantum mechanics may break down for large masses if one entertains the possibility that gravity can be “classical,” thereby leading to predictions different from conventional low-energy quantum gravity. Despite the technical difficulty in testing such deviations, a large number of experimental proposals have been put forward due to the high level of fundamental interest. Here, we consider the Schrödinger-Newton (SN) theory and the correlated worldline (CWL) theory, and show that they can be distinguished from conventional quantum mechanics, as well as each other, by performing pulsed optomechanics experiments. For CWL specifically we develop a framework resembling the commonly used “Heisenberg-picture” treatment of coupled oscillators, allowing one to perform simple calculations for such systems without delving into the deeper path-integral formalism. We find that discriminating between the theories will be very difficult until experimental control over low frequency quantum optomechanical systems is pushed much further. Furthermore, the predicted departures of SN and CWL from quantum mechanics occur at the same scale, so both alternative models could in principle be probed by a single experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Predicting core transport in ITER baseline discharges with neon injections

Achieving self-consistent performance predictions for ITER requires integrated modeling of core transport and divertor power exhaust under realistic impurity conditions. We present results from a systematic power-flow and impurity-content study for the ITER 15 MA baseline scenario constrained directly by existing SOLPS-ITER neon-seeded divertor solutions. Using the OMFIT STEP workflow, stationary temperature and density profiles are predicted with TGYRO for $1.5 \unicode{x2A7D} Z_\textrm{eff} \unicode{x2A7D} 2.5$, and the corresponding power crossing the separatrix $P_\textrm{sep}$ is evaluated. We find that $P_\textrm{sep}$ varies by more than a factor of 1.7 across this scan and matches the ${\sim}100$ MW SOLPS-ITER prediction when $Z_\textrm{eff} \simeq 1.6$ or when auxiliary heating is reduced to ${\sim}75\%$ of nominal. Rotation-sensitivity studies show that plausible variations in toroidal flow magnitude modify $P_\textrm{sep}$ by $\lesssim 20\%$, while AURORA modeling confirms that charge-exchange radiation inside the separatrix is dynamically negligible under predicted ITER neutral densities. These results identify a restricted compatibility window, $Z_\textrm{eff} \approx 1.6$ –1.75 and $0.75 \lesssim f_{P_\textrm{aux}} \unicode{x2A7D} 1.0$, in which core transport predictions remain aligned with neon-seeded divertor protection targets. This self-consistent, model-constrained framework provides actionable guidance for impurity control and auxiliary-heating scheduling in early ITER operation and supports future whole-device scenario optimization.

ITER↗

Microkinetic Modeling: A Tool for Rational Catalyst Design

The design of heterogeneous catalysts relies on understanding the fundamental surface kinetics that controls catalyst performance, and microkinetic modeling is a tool that can help the researcher in streamlining the process of catalyst design. Microkinetic modeling is used to identify critical reaction intermediates and rate-determining elementary reactions, thereby providing vital information for designing an improved catalyst. In this review, we summarize general procedures for developing microkinetic models using reaction kinetics parameters obtained from experimental data, theoretical correlations, and quantum chemical calculations. We examine the methods required to ensure the thermodynamic consistency of the microkinetic model. We describe procedures required for parameter adjustments to account for the heterogeneity of the catalyst and the inherent errors in parameter estimation. We discuss the analysis of microkinetic models to determine the rate determining reactions using the degree of rate control and reversibility of each elementary reaction. We introduce incorporation of Brønsted–Evans–Polanyi relations and scaling relations in microkinetic models and the effects of these relations on catalytic performance and formation of volcano curves. We review the analysis of reaction schemes in terms of the maximum rate of elementary reactions, and we outline a procedure to identify kinetically significant transition states and adsorbed intermediates. Here, we explore the application of generalized rate expressions for the prediction of optimal binding energies of important surface intermediates and to estimate the extent of potential rate improvement. We also explore the application of microkinetic modeling in homogeneous catalysis, electro-catalysis, and transient reaction kinetics. We conclude by highlighting the challenges and opportunities in the application of microkinetic modeling for catalyst design.

42 ENGINEERING↗

Benchmark study of a new simplified DFN model for shearing of intersecting fractures and faults

It is challenging to quantitatively predict shearing of intersecting fractures/faults because of dynamic frictional contacts accompanied by possible nonlinear rock deformation. To address such challenges, a new conceptual model—the simplified DFN model—was proposed and validated by Hu et al. 46 to use major paths (MPs) to represent complicated DFNs for calculation of shearing. In this work, we conducted a benchmark study for three examples that involve different levels of complexity of intersecting fractures, and correspondingly different numbers of MPs. The codes and software that were used in the benchmark cover a range of continuum, discontinuum and hybrid numerical methods: NMM (LBNL), FLAC3D (LBNL), GBDEM (KIGAM), FRACOD (DynaFrax), and CASRock (CAS). The general consistency between DFN and MP cases as predicted by all the codes/software demonstrates that major paths can be used to simplify the geometry of DFNs in a wide range of software. Disagreement in results made by some software and potential future improvements are discussed. We show that (1) shearing of one or multiple major fractures can be reduced if there are multiple smaller intersecting fractures in that area, which is a useful basis for understanding and controlling induced seismicity and merits further analysis, and (2) the agreement achieved in the benchmark examples provide confidence that the simplified DFN model is a promising conceptual model that can be used for different types of numerical approaches and software for simplifying the analysis of the shearing of intersecting fractures and faults.

58 GEOSCIENCES↗

Kinetic and X-ray Absorption Spectroscopic Analysis of Catalytic Redox Cycles over Highly Uniform Polymetal Oxo Clusters

Metal–organic framework materials (MOFs) offer an opportunity for investigating catalytic properties of polymetal oxo clusters that are highly well defined and uniform in nature, in contrast to other classes of catalysts that may exhibit a propensity toward active site heterogeneity. We report herein a kinetic and X-ray absorption spectroscopy (XAS) analysis of the two-electron oxidation of CO over divalent metal sites in MIL-100(M = Fe, Cr) (MIL = Materials of Institut Lavoisier) materials carrying μ 3 -oxo bridged trimers, and connect observations about the kinetic relevance of redox steps to density functional theory (DFT) predictions published previously. The high degree of uniformity evident from in situ titration measurements leads to a congruence in mechanistic inferences made from steady-state catalytic, transient stoichiometric, isotopic exchange, and isotopic tracer data that all point to a sequential mechanism comprised of separate oxidation and reduction half-cycle steps conjoined by an active oxygen intermediate. In situ XAS data reinforce mechanistic conclusions derived from kinetic analysis, and suggest that the active oxygen intermediate may be more appropriately characterized as an iron-oxyl (Fe 3+ –O – ) rather than an iron-oxo (Fe 4+ =O 2– ) species. The Cr analogue of MIL-100 exhibits contrasting rate features that can be rationalized using an identical sequence of steps as MIL-100(Fe), but with a highly dissimilar set of kinetic parameters that can also be validated using transient stoichiometric experiments. The larger coverages of active oxygen intermediates on MIL-100(Cr) are consistent with predictions from prior DFT studies that suggest more stable and less reactive active oxygen species for metals with lower d-electron counts, and point to metal identity as a lever for precisely controlling the kinetic relevance of oxidation and reduction half-cycles in catalytic redox sequences over polymetal oxo clusters. Importantly, the results presented point to the utility of developing broadly applicable structure–catalytic property relationships over MOF nodes specifically, and highly uniform catalysts more generally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of equilibrium pressure on plasma response to RMPs in a spherical tokamak

This study presents a comprehensive analysis of the equilibrium pressure on the plasma response to resonant magnetic perturbations (RMPs) in the spherical tokamak (ST) MAST-U, employing both single-fluid and MHD-kinetic hybrid models (implemented via the MARS-F/K codes). As a key finding, the study identifies two different pressure-driven eigenmodes, exhibiting Sturmian property, that affect the Troyon no-wall limits for the onset of the n = 1 and n = 2 ( n is the toroidal mode number) ideal external kink instabilities as well as the corresponding plasma response to the applied RMP. With increasing equilibrium pressure, the plasma response to RMPs is significantly enhanced in the ST plasma, particularly in the high-pressure regime where kinetic effects strongly stabilize the external kink instability. The Troyon no-wall limit divides the plasma response into two regions: well below the limit, the response amplitudes and trends (versus pressure) are similar between the fluid and kinetic models; as the equilibrium pressure approaches the Troyon limit, the kinetic model predicts significant amplification of the RMP field, up to 30 times for cases considered. A relatively weak dependence of the optimal coil phasing on the equilibrium pressure is computed in this ST plasma, similar to the trend obtained for the conventional aspect ratio devices. These findings underscore the importance of incorporating kinetic effects in accurate prediction of the plasma response to RMPs in high-pressure ST tokamak plasmas and provide a theoretical basis for optimizing RMP-based control of the edge-localized modes in future ST devices.

RMP↗

Systematic Construction of Time-Dependent Hamiltonians for Microwave-Driven Josephson Circuits

Time-dependent electromagnetic drives are fundamental for controlling complex quantum systems, including superconducting Josephson circuits. In these devices, accurate time-dependent Hamiltonian models are imperative for predicting their dynamics and designing high-fidelity quantum operations. Existing numerical methods, such as black-box quantization (BBQ) and energy-participation ratio (EPR), excel at modeling the static Hamiltonians of Josephson circuits. However, these techniques do not fully capture the behavior of driven circuits stimulated by external microwave drives, nor do they include a generalized approach to account for the inevitable noise and dissipation that enter through microwave ports. Here, we introduce numerical techniques that leverage classical microwave simulations, efficiently executable in finite-element solvers, to obtain the time-dependent Hamiltonian of microwave-driven superconducting circuits with arbitrary geometries under charge, flux, or mixed electromagnetic modulation. Importantly, our techniques do not rely on a lumped-element description of the superconducting circuit, in contrast to previous approaches to tackling this problem. We demonstrate the versatility of our approach by characterizing the driven properties of realistic circuit devices in complex electromagnetic environments, including coherent dynamics due to charge and flux modulation, as well as drive-induced relaxation and dephasing. Our techniques offer a powerful toolbox for optimizing circuit designs and advancing practical applications in superconducting quantum computing.

Lu, Yao [Yale U.; Yale U. (main); Fermilab] (ORCID↗