SEARCH · Engineering Papers

Results for “criticality simulations”

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 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Coupling localized Noah-MP-Crop model with the WRF model improved dynamic crop growth simulation across Northeast China

Croplands play a critical role in regulating the energy and moisture exchanges between the land surface and atmosphere. However, the interactions between cropland and climate are usually poorly represented due to a lack of detailed representation in crop types and field management. Here, we coupled the Noah-MP-Crop model with the state-of-the-art Weather Research and Forecasting (WRF) model to explore and evaluate the crop growth dynamics in response to climate variations across Northeast China. The default parameters of the crop model were not exactly suitable for the agricultural ecosystems in Northeast China. The detailed cropland distribution, and crop phenology parameters including growing degree days (GDD) and planting (harvesting) date were first created using multi-source remote sensing products and reanalysis data, and was then successfully used to simulate the growth and yield for corn and soybean and associated energy exchanges. We also optimized and calibrated other crop parameters using the time-series of the Moderate Resolution Imaging Spectroradiometer (MODIS) land surface products. The modified crop model substantially improved the simulation of crop growth, plant physiology, and biomass accumulation for both corn and soybean. Coupling the localized dynamic crop model into the WRF led to considerable decreases in the simulated mean-absolute-errors (MAEs) and biases of the leaf area index, evapotranspiration, and gross primary production compared with the MODIS observed values. Compared with the statistical yield from each province, the modified crop model underestimated the corn yield from 11.1% to 48.6%, whereas overestimated the soybean yield from 16.5% to 162.6%.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1993

Numerical Study of Base Pressure Characteristic Curve for a Four-Engine Clustered Nozzle Configuration

Excessive base heating has been a problem for many launch vehicles. For certain designs such as the direct dump of turbine exhaust in the nozzle section and at the nozzle lip of the Space Transportation Systems Engine (STME), the potential burning of the turbine exhaust in the base region has caused tremendous concern. Two conventional approaches have been considered for predicting the base environment: (1) empirical approach, and (2) experimental approach. The empirical approach uses a combination of data correlations and semi-theoretical calculations. It works best for linear problems, simple physics and geometry. However, it is highly suspicious when complex geometry and flow physics are involved, especially when the subject is out of historical database. The experimental approach is often used to establish database for engineering analysis. However, it is qualitative at best for base flow problems. Other criticisms include the inability to simulate forebody boundary layer correctly, the interference effect from tunnel walls, and the inability to scale all pertinent parameters. Furthermore, there is a contention that the information extrapolated from subscale tests with combustion is not conservative. One potential alternative to the conventional methods is computational fluid dynamics (CFD), which has none of the above restrictions and is becoming more feasible due to maturing algorithms and advancing computer technology. It provides more details of the flowfield and is only limited by computer resources. However, it has its share of criticisms as a predictive tool for base environment. One major concern is that CFD has not been extensively tested for base flow problems. It is therefore imperative that CFD be assessed and benchmarked satisfactorily for base flows. In this study, the turbulent base flowfield of a experimental investigation for a four-engine clustered nozzle is numerically benchmarked using a pressure based CFD method. Since the cold air was the medium, accurate prediction of the base pressure distributions at high altitudes is the primary goal. Other factors which may influence the numerical results such as the effects of grid density, turbulence model, differencing scheme, and boundary conditions are also being addressed. Preliminary results of the computed base pressure agreed reasonably well with that of the measurement. Basic base flow features such as the reverse jet, wall jet, recompression shock, and static pressure field in plane of impingement have been captured.

Wang, Ten-See↗

DOE OSTI2020

A Clustering-Based Scenario Generation Framework for Power Market Simulation with Wind Integration

A critical step in stochastic optimization models of power system analysis is to select a set of appropriate scenarios and significant numbers of scenario generation methods exist in the literature. This paper develops a clustering based scenario generation method, which aims to improve the performance of existing scenario generation techniques by grouping a set of correlated wind sites into clusters according to their cross-correlations. Copula based models are utilized to model spatiotemporal correlations and the Gibbs sampling is then used to generate scenarios for day-ahead markets. Our results show that the generated scenarios based on clustered wind sites outperform existing approaches in terms of reliability and sharpness and can reduce the total computational time for scenario generation and reduction significantly. The clustering-based framework can therefore provide a better support for real-world market simulations with high wind penetration.

data visualization↗

DOE OSTI2022

Implementation and evaluation of the unified stomatal optimization approach in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES)

Abstract. Stomata play a central role in regulating the exchange of carbon dioxide and water vapor between ecosystems and the atmosphere. Their function is represented in land surface models (LSMs) by conductance models. The Functionally Assembled Terrestrial Ecosystem Simulator (FATES) is a dynamic vegetation demography model that can simulate both detailed plant demographic and physiological dynamics. To evaluate the effect of stomatal conductance model formulation on forest water and carbon fluxes in FATES, we implemented an optimality-based stomatal conductance model – the Medlyn (MED) model – that simulates the relationship between photosynthesis (A) and stomatal conductance to water vapor (gsw) as an alternative to the FATES default Ball–Woodrow–Berry (BWB) model. To evaluate how the behavior of FATES is affected by stomatal model choice, we conducted a model sensitivity analysis to explore the response of gsw to climate forcing, including atmospheric CO2 concentration, air temperature, radiation, and vapor pressure deficit in the air (VPDa). We found that modeled gsw values varied greatly between the BWB and MED formulations due to the different default stomatal slope parameters (g1). After harmonizing g1 and holding the stomatal intercept parameter (g0) constant for both model formulations, we found that the divergence in modeled gsw was limited to conditions when the VPDa exceeded 1.5 kPa. We then evaluated model simulation results against measurements from a wet evergreen forest in Panama. Results showed that both the MED and BWB model formulations were able to capture the magnitude and diurnal changes of measured gsw and A but underestimated both by about 30 % when the soil was predicted to be very dry. Comparison of modeled soil water content from FATES to a reanalysis product showed that FATES captured soil drying well, but translation of drying soil to modeled physiology reduced the models' ability to match observations. Our study suggests that the parameterization of stomatal conductance models and current model response to drought are the critical areas for improving model simulation of CO2 and water fluxes in tropical forests.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2002

A critical period for functional vestibular development in zebrafish

We have determined a critical period for vestibular development in zebrafish by using a bioreactor designed by NASA to simulate microgravity for cells in culture. A critical period is defined as the briefest period of time during development when stimulus deprivation results in long lasting or permanent sensory deficits. Zebrafish eggs were collected within 3 hours of being laid and fertilized. In experiment 1, eggs were placed in the bioreactor at 3, 24, 30, 36, 48, or 72 hours postfertilization (hPF) and maintained in the bioreactor until 96 hPF. In experiment 2, eggs were placed in the bioreactor immediately after they were collected and maintained in the bioreactor until 24, 36, 48, 60, 66, 72, or 96 hPF. Beginning at 96 hPF, all larvae had their vestibulo-ocular reflexes (VOR) evaluated once each day for 5 days. Only larvae that hatched from eggs that were placed in the bioreactor before 30 hPF in experiment 1 or removed from the bioreactor later than 66 hPF in experiment 2 had VOR deficits that persisted for at least 5 days. These data suggest a critical period for vestibular development in the zebrafish that begins before 30 hPF and ends after 66 hPF. To confirm this, zebrafish eggs were placed in the bioreactor at 24 hPF and removed at 72 hPF. VORs were evaluated in these larvae once each day for 5 days beginning at 96 hPF. These larvae had VOR deficits that persisted for at least 5 days. In addition, larvae that had been maintained in the bioreactor from 24 to 66 hPF or from 30 to 72 hPF, had only temporary VOR deficits. In a final experiment, zebrafish eggs were placed in the bioreactor at 3 hPF and removed at 96 hPF but the bioreactor was turned off from 24 hPF to 72 hPF. These larvae had normal VORs when they were removed from the bioreactor at 96 hPF. Taken as a whole, these data support the idea that there is a critical period for functional maturation of the zebrafish vestibular system. The developmental period identified includes the timeframe during which the vestibular primary afferent neurons are born, innervate their central and peripheral targets, and remodel their central projections. Copyright 2002 Wiley-Liss, Inc.

Non-NASA Center↗

DOE OSTI2022

Advantages of a variable‐resolution global climate model in reproducing the seasonal evolution of East Asian summer monsoon

Abstract The East Asian summer monsoon (EASM) is unique among monsoon systems that it features meridional evolution of the summer monsoon. In this study, we evaluate the performances of a Variable‐Resolution Community Earth System Model (VR‐CESM) regionally refined over eastern China (14 km) in reproducing the seasonal evolution of EASM precipitation over China. Compared with reference datasets, VR‐CESM shows better performance than the corresponding globally uniform coarse‐resolution model CESM (quasi‐uniform 1°), especially over western China where complex local topography exists. The northward monsoon migration is closely related to low‐level southerly flows and vertical moisture advection, which are more reasonably simulated in VR‐CESM. The four critical timings of the EASM (monsoon onset, withdrawal, peak, and duration) are also better captured in VR‐CESM than in CESM. The corresponding spatial Pearson correlation coefficients of the four critical timings with respect to reference datasets are about 0.1 higher in VR‐CESM than those in CESM. Both models are most accurate in simulating monsoon onset and least accurate at simulating the monsoon peak. The overestimated zonal thermal contrast in CESM is responsible for the earlier monsoon onset and excessive precipitation in September over the Yangtze River valley. Finer resolution in VR‐CESM, especially over the Tibetan Plateau (TP), appears to be a main factor in simulating better zonal thermal contrast and seasonal evolution of the EASM.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Final Design of a New Pulsed-Neutron Die-Away Experimental Testbed for Thermal Scattering Law Benchmarks (PNDA) (IER-501 CED-2)

This report presents the final design (CED-2) of IER-501, a testbed for pulsed-neutron die-away (PNDA) experiments at Lawrence Livermore National Laboratory (LLNL). The PNDA experiments will be high-quality benchmarks that serve to optimize and validate thermal neutron scattering laws (TSLs), which are high-priority nuclear data for the Department of Energy’s Nuclear Criticality Safety Program (NCSP). This CED-2 was guided by experimental results gathered from a prototype experimental set-up rather than with detailed simulations as is typically done for criticality experiments. The report presents the PNDA design and its equipment, and it documents the methodology used to postprocess the measured data. It discusses how experimental configurations influenced the decay-constant eigenvalue, or the integral parameter, of the experiment. The studied variations included the distance from the target to the generator, the orientation of the generator relative to the box, the orientation of the target to the generator, the effect of HDPE lining in the shielding box, the settings of the generator, and finally the benefits of using a low-scatter facility. From these results, a final configuration for the target, detectors, shielding box, and neutron generator was chosen. Two factors particularly influenced the chosen configuration: the reproducibility of the experiment’s geometry for benchmark modeling, and the minimization of measured neutrons scattered off the shielding box.

42 ENGINEERING↗

DOE OSTI2022

Including Chi-Nu 235 U PFNS Experimental Data into an ENDF/B-VIII.1 Release Candidate Evaluation

This report documents an evaluation of 235 U prompt fission neutron spectra (PFNS) that is a release candidate for the upcoming U.S. nuclear data library, ENDF/B-VIII.1. This evaluation differs from its predecessor, ENDF/B-VIII.0, mainly by the inclusion of 235 U PFNS measured by the Chi-Nu team of LANL and LLNL. This data set is the first one that covers the 235 U PFNS for continuous incident-neutron energies of 1⁻20 MeV and outgoing-neutron energies from 10 keV⁻10 MeV with high precision. Previous data sets were either measured in a limited energy range or with less precision. Hence, these new Chi-Nu data provide decisive information for the evaluation. The resulting evaluated data correspond well to the new experimental PFNS. The evaluated PFNS also produce average mean energies and 239 Pu/ 235 U PFNS in agreement with associated Chi-Nu data. If one uses the new evaluated data to predict the neutron multiplication factor, k eff , of selected ICSBEP critical assemblies, the differences of simulated values compared to those using ENDF/B-VIII.0 is modest (less than 55 pcm). This difference in k eff can be easily accommodated by changes in the 235 U average prompt fission neutron multiplicity that is currently being re-evaluated. In addition to that, the new PFNS predict on average 235 U LLNL pulsed-sphere neutron-leakage spectra better than ENDF/B-VIIII.0 PFNS.

235U↗

NASA NTRS1999

Monte Carlo Simulation of Nanoparticle Encapsulation in Flames

Two critical challenges facing the application of flames for synthesis of nanopowder materials are: (1) overcoming formation of agglomerates and (2) ensuring that the highly reactive nanopowders that are synthesized in flames can be produced in such a manner that their purity is maintained during subsequent processing. Agglomerates are produced in flames because particle formation occurs in a high temperature and high number density environment. They are undesirable in most advanced applications of powders. For example, agglomerates have a deleterious effect on compaction density, leading to voids when nanopowders are consolidated. Efforts to avoid agglomeration in flames without substantially reducing particle number density and, consequently, production rate, have had limited success. Powder purity must also be maintained during subsequent handling of nanopowders and this poses a significant challenge for any synthesis route because nanopowders, particularly metals and non-oxide ceramic powders, are inherently reactive. Impurities acquired during handling of nanopowders have slowed the advancement of the nanostructured materials industry. One promising approach that has been proposed to address these problems is nano-encapsulation. In this approach, the core particles are encapsulated in a removable material while they are within the flame but before excessive agglomeration has occurred. Condensation can be very rapid so that core particles are trapped within the condensed material and agglomeration is limited. Nano-encapsulation also addresses the handling concerns for post-synthesis processing. Results have shown that when nano-encapsulated powders are exposed to atmosphere the core particles are protected from oxidation and/or hydrolysis. Thus, handling of the powders does not require extreme care. If, for example, at the time of consolidation the encapsulation material is removed by vacuum annealing, the resulting powder remains unagglomerated and free of impurities. In this work, we described a novel aerosol model that has been developed to simulate particle encapsulation in flames. The model will ultimately be coupled to a one-dimensional spherical flame code and compared to results from microgravity flame experiments.

Sun, Z.↗

DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1995

Optimal Trajectories and Control Strategies for the Helicopter in One-Engine-Inoperative Terminal-Area Operations

Engine failure represents a major safety concern to helicopter operations, especially in the critical flight phases of takeoff and landing from/to small, confined areas. As a result, the JAA and FAA both certificate a transport helicopter as either Category-A or Category-B according to the ability to continue its operations following engine failures. A Category-B helicopter must be able to land safely in the event of one or all engine failures. There is no requirement, however, for continued flight capability. In contrast, Category-A certification, which applies to multi-engine transport helicopters with independent engine systems, requires that they continue the flight with one engine inoperative (OEI). These stringent requirements, while permitting its operations from rooftops and oil rigs and flight to areas where no emergency landing sites are available, restrict the payload of a Category-A transport helicopter to a value safe for continued flight as well as for landing with one engine inoperative. The current certification process involves extensive flight tests, which are potentially dangerous, costly, and time consuming. These tests require the pilot to simulate engine failures at increasingly critical conditions, Flight manuals based on these tests tend to provide very conservative recommendations with regard to maximum takeoff weight or required runway length. There are very few theoretical studies on this subject to identify the fundamental parameters and tradeoff factors involved. Furthermore, a capability for real-time generation of OEI optimal trajectories is very desirable for providing timely cockpit display guidance to assist the pilot in reducing his workload and to increase safety in a consistent and reliable manner. A joint research program involving NASA Ames Research Center, the FAA, and the University of Minnesota is being conducted to determine OEI optimal control strategies and the associated optimal,trajectories for continued takeoff (CTO), rejected takeoff (RTO), balked landing (BL), and continued landing (CL) for a twin engine helicopter in both VTOL and STOL terminal-area operations. This proposed paper will present the problem formulation, the optimal control solution methods, and the key results of the trajectory optimization studies for both STOL and VTOL OEI operations. In addition, new results concerning the recently developed methodology, which enable a real-time generation of optimal OEI trajectories, will be presented in the paper. This new real-time capability was developed to support the second piloted simulator investigation on cockpit displays for Category-A operations being scheduled for the NASA Ames Vertical Motion Simulator in June-August of 1995. The first VMS simulation was conducted in 1994 and reported.

Chen, Robert T. N.↗

DOE OSTI2021

Evaluating the performance of WRF urban schemes and PBL schemes over Dallas-Fort Worth during a dry summer and a wet summer

This study evaluated the Weather Research and Forecasting (WRF) model sensitivity to different planetary boundary layer (PBL) schemes (the YSU and MYJ schemes) and urban schemes including the bulk scheme (BULK), single-layer urban canopy model (UCM), multi-layer building environment parameterization (BEP) model, and multi-layer building energy model (BEM). Daily reinitialization simulations were conducted over Dallas-Fort Worth during a dry summer month (July 2011) and a wet summer month (July 2015) with weaker (stronger) daytime (nocturnal) UHI in 2011 than 2015. All urban schemes overestimated the urban daytime 2m temperature in both summers, but BEP and BEM still reproduced the daytime urban cool island in dry summer. All urban schemes reproduced the nocturnal urban heat island, with BEP producing the weakest one due to its unrealistic urban cooling. BULK and UCM overestimated the urban canopy wind speed, while BEP and BEM underestimated it. The urban schemes showed prominent impact on daytime PBL profiles. UCM+MYJ showed a superior performance than other configurations. The relatively large (small) aspect ratio between building height and road width in UCM (BEM) was responsible for the overprediction (underprediction) of urban canopy temperature. The relatively low (high) building height in UCM (BEM) was responsible for the overprediction (underprediction) of urban canopy wind speed. Improving urban schemes and providing realistic urban parameters were critical for improving urban canopy simulation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Application of the Rossi-alpha method to simulations of HEU and organic scintillators

The International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook contains thousands of high-quality benchmark measurements that include detailed descriptions of experiment designs and exhaustive accounting of uncertainties. These benchmarks are used extensively by the criticality safety and nuclear data communities. A subcritical benchmark involving bare, highly enriched uranium (HEU) has yet to be accepted to the handbook. The Measurement of Uranium Subcritical and Critical (MUSiC) was proposed to address this gap. MUSiC will be performed at the National Criticality Experiments Research Center (NCERC) and will leverage bare HEU hemisphere shells that can be assembled in a variety of configurations. The measurement has the following key benefits: it will be the first set of fast, subcritical measurements of HEU to be submitted to ICSBEP; ten configurations will be measured spanning from nearly no multiplication to critical (e.g., it is a measurement of similar systems with different re activities); measurements of an all-HEU system will be of value for nuclear data validation, providing results that assist in determining biases in both measurements and simulations. A variety of detector systems will be used including a Neutron Multiplicity Array Detector (NoMAD) system (similar to the MC-15), four small volume 0.635 cm (#8960;) × 7.59 cm 3He detectors (ideal for measuring prompt neutron decay constants due to their fast recovery speed), and an array of eight 5.08 cm (#8960;) × 5.08 cm EJ-309 organic scintillator detectors. The focus of this work is the simulation of the subcritical MUSiC configurations including the EJ-309 detectors and the application of the Rossi-α method. The EJ-309 organic scintillator detector array is a new detection system that will be deployed at NCERC. Organic scintillators are ideal for measurements of fast, subcritical systems due to their short time response on the order of nanoseconds.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Defining quantum-ready primitives for hybrid HPC-QC supercomputing: a case study in Hamiltonian simulation

As computational demands in scientific applications continue to rise, hybrid high-performance computing (HPC) systems integrating classical and quantum computers (HPC-QC) are emerging as a promising approach to tackling complex computational challenges. One critical area of application is Hamiltonian simulation, a fundamental task in quantum physics and other large-scale scientific domains. This paper investigates strategies for quantum-classical integration to enhance Hamiltonian simulation within hybrid supercomputing environments. By analyzing computational primitives in HPC allocations dedicated to these tasks, we identify key components in Hamiltonian simulation workflows that stand to benefit from quantum acceleration. To this end, we systematically break down the Hamiltonian simulation process into discrete computational phases, highlighting specific primitives that could be effectively offloaded to quantum processors for improved efficiency. Our empirical findings provide insights into system integration, potential offloading techniques, and the challenges of achieving seamless quantum-classical interoperability. We assess the feasibility of quantum-ready primitives within HPC workflows and discuss key barriers such as synchronization, data transfer latency, and algorithmic adaptability. These results contribute to the ongoing development of optimized hybrid solutions, advancing the role of quantum-enhanced computing in scientific research.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Quantum Information for Fusion Energy Sciences (Final Technical Report)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Quantum Information for Fusion Energy Sciences (Final Technical Report)

The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer. In order to efficiently simulate plasma physics, an optimal control method which synthesizes, directly at the pulse level, any quantum gate on qubit and qutrit systems was developed. Using four superconducting transmon quantum processors at Rigetti and LLNL, it was demonstrated that any arbitrary quantum gate on qubits and qutrits could be implemented with high fidelity, leading to a significantly reduced length of a gate sequence. A problem of interest in FES is the nonlinear optical process of laser pulse compression within a plasma. Since quantum physics is linear, simulating nonlinear operations is not naturally feasible on a quantum computer, however it is possible to simulated a quantized version of the nonlinear process. A quantization approach to convert nonlinear wave-wave interaction problems to Hamiltonian simulation problems was developed and demonstrated using two qubits on a Rigetti device. In this experiment, a number of error suppression and mitigation techniques were investigated to determine how best to utilize the finite quantum resources. This study provides an example of how plasma problems may be solved on near-term, noisy quantum computing platforms and identified a promising set of techniques. Building on the insights of these experiments, the investigation turned to linear electron-plasma wave physics. A connection was identified between a local one-dimensional lattice spin model and linear wave phenomena, allowing a plasma physics problem to be efficiently mapped to the quantum computer. In this framework, reflection and transmission of plasma waves at a sharp boundary was studied, as well as the propagation of waves through an inhomogeneous plasma medium. In addition to the suite of error suppression and mitigation techniques developed, this experiment introduced the use of a digital-analog gate scheme designed to efficiently simulate the plasma Hamiltonian. With hardware available at the conclusion of the project, simulation at the scale of 9 qubits and 15 timesteps (60 entangling layers) was achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Nonlinear dynamics and dissipation of a curvilinear vortex driven by a strong time-dependent Meissner current

Here, we report numerical simulations of large-amplitude oscillations of a trapped vortex line under a strong ac magnetic field $\textit{H(t) = H}$ sin$\textit{ωt}$ parallel to the surface. The power dissipated by an oscillating vortex segment driven by the surface ac Meissner currents was calculated by taking into account the nonlinear vortex line tension, vortex mass, and a nonlinear Larkin-Ovchinnikov (LO) viscous drag coefficient $\textit{η(v)}$. We show that the LO decrease of $\textit{η(v)}$ with the vortex velocity $\textit{v}$ can radically change the field dependence of the surface resistance R i ( H ) caused by trapped vortices. At low frequencies, $R_i(H)$ exhibits a conventional increase with $\textit{H}$, but as $\textit{ω}$ increases, the surface resistance becomes a nonmonotonic function of $\textit{H}$ which decreases with $\textit{H}$ at higher fields. The effects of frequency, pin spacing, and the mean-free path $l_i$ on the field dependence of $R_i(H)$ were calculated. It is shown that, as the surface gets dirtier and $l_i$ decreases, the anomalous drop of $R_i(H)$ with $\textit{H}$ shifts to lower fields which can be much smaller than the lower critical magnetic field. Our numerical simulations also show that the LO decrease of $\textit{η(v)}$ with $\textit{v}$ can cause a vortex bending instability at high field amplitudes and frequencies, giving rise to the formation of dynamic kinks along the vortex. Measurements of $R_i(H)$ caused by sparse vortices trapped perpendicular to the surface can offer opportunities to investigate an extreme nonlinear dynamics of vortices driven by strong current densities up to the depairing limit at low temperatures. The behavior of $R_i(H)$, which can be tuned by varying the rf frequency or concentration of nonmagnetic impurities, is not masked by strong heating effects characteristic of dc or pulse transport measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Probing the electronic properties of the electrified silicon/water interface by combining simulations and experiments.

Silicon (Si) is broadly used in electrochemical and photoelectrochemical devices, where the capacitive and Faradaic reactions at the Si/water interfaces are critical for signal transduction or noise generation. However, probing the electrified Si/water interface at the microscopic level remains a challenging task. Here we focus on hydrogenated Si surfaces in contact with water, relevant to transient electronics and photoelectrochemical modulation of biological cells and tissues. We show that by carrying out first principles molecular dynamics simulations of the Si(100)/water interface in the presence of an electric field we can realistically correlate the computed flat-band potential and tunneling current images at the interface with experimentally measured capacitive and Faradaic currents. Specifically, we validate our simulations in the presence of bias by performing pulsed chronoamperometry measurements on Si wafers in solution. Consistent with prior experiments, our measurements and simulations indicate the presence of voltage-dependent capacitive currents at the interface. We also find that Faradaic currents are weakly dependent on the applied bias, which we relate to surface defects present in newly prepared samples.

Ye, Zifan↗