Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated current”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Microfabricated ion trap chip with an integrated microwave antenna

An ion trap chip, which may be used for quantum information processing and the like, includes an integrated microwave antenna. The antenna is formed as a radiator connected by one of its ends to the center trace of a microwave transmission line and connected by its other end to a current return path through a ground trace of the microwave transmission line. The radiator includes several parallel, coplanar radiator traces connected in series. The radiator traces are connected such that they all carry electric current in the same direction, so that collectively, they simulate a single, unidirectionally flowing sheet of current. In embodiments, induced currents in underlying metallization planes are suppressed by parallel slots that extend in a direction perpendicular to the radiator traces.

Nordquist, Christopher↗
DOE OSTI2026

Electrostatic Gating of Ionic Conductance through Heterogeneous van der Waals Nanopores

Nanofluidic ionic transistors typically require gate voltages above 1 V and operate only at submillimolar ionic strengths, limiting their biocompatible applications. We demonstrate ionic transistors consisting of single sub-10 nm nanopores drilled in van der Waals (vdW) heterostructures with internal gate electrodes made of few-layer graphene. These devices deliver up to 10-fold current modulation at gate voltages as low as 0.3 V in 10 mM KCl, and ∼2-fold modulation at near-physiological 100 mM KCl. Baseline conductance with no gate shows surface-charge-dominated transport below 100 mM KCl, consistent with negatively charged hBN walls and ∼5 nm opening of the pores. The surface charge and the electrochemical asymmetry introduced by the three-electrode configuration govern the device’s behavior: negative gate voltage (V G ) enriches ionic concentrations and enhances current, whereas positive V G induces a local depletion zone that suppresses transport. The current modulation by V G is dependent on the polarity of the transmembrane potential and leads to ion current rectification. Molecular dynamics simulations of a nanopore in a hBN–graphene–hBN stack reveal confinement and surface charge-dependent suppression of the relative permittivity of interfacial water. Continuum modeling with radially varying interfacial water permittivity reproduces the asymmetric I–V characteristics and explains how the embedded gate sculpts local potential and ion concentrations. By enabling sub-0.5 V control of ionic transport at up to 100 mM salt concentrations, these devices address a key need in nanofluidics to create low-power ionic circuits and biosensing.

Materials science↗
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗
DOE OSTI2020

Machine learning and ligand binding predictions: A review of data, methods, and obstacles

We report that computational predictions of ligand binding is a difficult problem, with more accurate methods being extremely computationally expensive. The use of machine learning for drug binding predictions could possibly leverage the use of biomedical big data in exchange for time-intensive simulations. This paper reviews current trends in the use of machine learning for drug binding predictions, data sources to develop machine learning algorithms, and potential problems that may lead to overfitting and ungeneralizable models. A few popular datasets that can be used to develop virtual high-throughput screening models are characterized using spatial statistics to quantify potential biases. We can see from evaluating some common benchmarks that good performance correlates with models with high-predicted bias scores and models with low bias scores do not have much predictive power. A better understanding of the limits of available data sources and how to fix them will lead to more generalizable models that will lead to novel drug discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Cluster analysis of occupancy schedules in residential buildings in the United States

The energy performance of residential buildings significantly depends on the building occupants’ behavior, which can be highly variable. When the heating, ventilation and air conditioning (HVAC) system is controlled based on the presence or absence of occupants in a building, occupant behavior is of even further importance to its energy performance. In current practice, building energy simulation tools generally use a single occupancy profile to represent the building’s occupancy schedule, the schedule of which is considered to be the same, regardless of the type of household being modeled. Thus, there is significant potential for improvement to allow for more flexibility and accuracy in calculation of occupancy. The objective of this study is to assess the variations in the typical types of occupancy schedules followed by the U.S. population using cluster analysis. American Time Use Survey data, which statically represents the overall U.S. population’s activities, across 12 years (2006–2017), is used. The ATUS data is segregated into smaller groups based on age and weekday/weekend, then divided into activities that are considered “at home” and “away from home”, which are mapped to the presence or non-presence of occupants in the home. Cluster analysis is then used to identify common types of occupancy schedule patterns for each age group. Three main types of patterns are obtained from cluster analysis for each age group, which together represent approximately 88% of people in the United States. The output of the cluster analysis is further analyzed to evaluate the variation in characteristics, including the number of times leaving home, time of day when leaving the home, and the timespan of absence from the home. The results of this study provide detailed insights on how typical occupants in the United States spend their time in residential spaces which can be used to create occupancy profiles for residential buildings. Finally, these occupancy profiles could be utilized inform an assessment of the energy use impact of occupancy-based controls of energy consuming systems and technologies.

42 ENGINEERING↗
DOE OSTI2021

Coupled physics simulation of fracture in nuclear fuel pellets induced by resistive heating

While it is well known that the UO 2 fuel used in light-water reactors cracks early in its life in the reactor, it is very difficult to collect in-reactor data regarding fuel cracking initiation and propagation to validate fracture models in fuel performance codes. Most fracture data comes from the post-irradiation examination of fuel that has been subjected to long and complex power histories. Replicating the thermal gradients experienced in this fuel to study crack initiation and propagation in a laboratory environment is challenging. This is primarily because the fuel is heated volumetrically by fission, which is difficult to replicate. One way to approximate this volumetric heating out of the reactor is to use electrical resistance (or Joule heating). An electrical conductivity model has been developed and coupled with the existing thermomechanical modeling capability in the BISON nuclear fuel performance code. This will permit simulation of experiments that use resistive heating for the purpose of validating fracture models in BISON. The modeling approach developed here is applied to simulate both historical and current resistive heating experiments on ceramic nuclear fuel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Atomic masses with machine learning for the astrophysical r process

The astrophysical r process plays a vital role in the production of heavy elements. Modeling of the r process is sensitive to masses and further requires knowledge of masses beyond current experimental reach. Therefore, simulations of the r process offer a unique test bed for predicting mass extrapolations. We take a Machine-Learning (ML) approach to model the masses across the entire chart of nuclides. For the first time, we simulate r-process nucleosynthesis with a physics-based ML mass model. We compare simulated abundances to solar data in order to evaluate the model's performance far from stability. The resulting r-process abundances up to thorium and uranium qualitatively match those of the observed solar system abundance pattern, with the characteristic peaks well positioned. We propagate the mass uncertainties obtained from the ML model to r-process abundance yields to estimate an uncertainty band associated with our approach. The size of the uncertainty band is approximately one order of magnitude which aligns with the uncertainty reported using alternative techniques.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Effect of Glacial/Interglacial Recharge Conditions on Flow of Meteoric Water Through Deep Orogenic Faults: Insights Into the Geothermal System at Grimsel Pass, Switzerland

Many meteoric-recharged, fault-hosted geothermal systems in amagmatic orogenic belts have been active through the Pleistocene glacial/interglacial climate fluctuations. The effects of such climate-induced recharge variations on fluid flow patterns and residence times of the thermal waters are complex and may influence how the geothermal and mineralization potential of the systems are evaluated. We report systematic thermal–hydraulic simulations designed to reveal the effects of recharge variations, using a model patterned on the orogenic geothermal system at Grimsel Pass in the Swiss Alps. Fault-bound circulation of meteoric water through the Grimsel fault is driven to depths of ~10 km by the high alpine topography. Simulations suggest that the current single pass flow is typical of interglacial periods, during which 1) meteoric recharge into the fault is high (above tens of centimeters per year), 2) conditions are at or somewhat below the critical Rayleigh number, and 3) the hydraulic connectivity along the fault plane is extensive (an extent of at least 10 km into increasingly higher terrain is required to explain the 10 km penetration depth). The subcritical condition constrains the bulk fault permeability to < 1e-14 m2. In contrast, the limited recharge during the numerous Pleistocene glaciation events likely induced a layered flow system with single pass flow confined to shallow depths while non-Rayleigh convection occurred in the deep fault. The same layering can be observed at low aspect ratios (length/depth) of the fault plane, when the available recharge area limits the flux through the fault.

58 GEOSCIENCES↗
DOE OSTI2022

Reconciling Observed and Predicted Tropical Rainforest OH Concentrations

We present OH observations made in Amazonas, Brazil during the Green Ocean Amazon campaign (GoAmazon2014/5) from February to March of 2014. The average diurnal variation of OH peaked with a midday (10:00–15:00) average of 1.0 × 10 6 (±0.6 × 10 6 ) molecules cm -3 . This was substantially lower than previously reported in other tropical forest photochemical environments (2–5 × 10 6 molecules cm -3 ) while the simulated OH reactivity was lower. The observational data set was used to constrain a box model to examine how well current photochemical reaction mechanisms can simulate observed OH. We used one near-explicit mechanism (MCM v3.3.1) and four condensed mechanisms (i.e., RACM2, MOZART-T1, CB05, CB6r2) to simulate OH. A total of 14 days of analysis shows that all five chemical mechanisms were able to explain the measured OH within instrumental uncertainty of 40% during the campaign in the Amazonian rainforest environment. Future studies are required using more reliable NO x and VOC measurements to further investigate discrepancies in our understanding of the radical chemistry in the tropical rainforest.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Solar cell contacts: quantifying the impact of interfacial layers on selectivity, recombination, charge transfer, and V oc

Interfacial layers (IFLs) are ubiquitous in solar cells, but their precise impact on carrier transfer rates and the relation of these rates to performance metrics and the concepts of selectivity and recombination is lacking. We report the use of a well-defined interdigitated back-contact (IBC) silicon solar cell to determine the precise role of such IFLs. We characterize the action of the common IFL spiro-OMeTAD by making it a third contact to the IBC cell. This architecture creates three solar cells within a single structure that, with numerical simulation, provide the exchange current densities (i.e., charge transfer rates) for electrons ( J 0n ) and holes ( J 0p ) and the quasi Fermi-level splitting in the absorber, which measures the balance of generation and recombination. Further, we describe the relation of V oc to contact recombination, the asymmetry in electron/hole collection rates at a single contact (contact selectivity), and the asymmetry in collection rates of the same carrier at separate contacts (carrier selectivity). Relative to bare gold, neat spiro-OMeTAD reduces J 0n and J 0p (their geometric mean ( J 0n J 0p ) 0.5 decreases by 10 4 ), decreasing contact recombination. Addition of the common dopant Li-TFSI and air increase J 0p / J 0n by 10 6 with little effect on ( J 0n J 0p ) 0.5 , increasing contact hole selectivity. The significant increase in V oc observed by introducing spiro-OMeTAD/Li-TFSI IFLs into the cells studied, however, is due to an increase in the carrier selectivity rather than the contact selectivity or recombination of the spiro-OMeTAD-modified contact. Operando measurements further show voltage-dependent changes in the J 0 s, demonstrating that spiro-OMeTAD contributes to hysteresis. Broader context : Solar cells operate by photogenerating excess charge carriers in an absorber material and, in competition with recombination, asymmetrically extracting them at so-called carrier-selective contacts, one that ideally collects only electrons and the other, only holes. Particularly in emerging photovoltaics such as perovskites, thin layers of organic semiconductors or related materials are introduced between the absorber and contact to improve power conversion efficiency. In terms of interfacial charge transfer, a prevailing view is that such interfacial layers improve performance by helping block the collection of the undesired carrier, considered a form of recombination. We use a novel platform to study the simultaneous impact of spiro-OMeTAD, a common interfacial layer, on the collection of electrons, the collection of holes, and on the recombination of electrons and holes. We quantitatively demonstrate that spiro-OMeTAD layers indeed passivate the contact toward recombination, but that larger improvements in the open-circuit voltage, a key cell metric, can come not only from this but also from increasing the asymmetry of the collection of electrons in the system. Further, operando measurements show transient changes in the properties of spiro-OMeTAD which suggest that it contributes to hysteresis phenomena commonly observed in perovskite and other solar cells.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Understanding and predicting damage and failure at grain boundaries in BCC Ta

Understanding the effect of grain boundaries (GBs) on the deformation and spall behavior is critical to designing materials with tailored failure responses under dynamic loading. This understanding is hampered by the lack of in situ imaging capability with the optimum spatial and temporal resolution during dynamic experiments, as well as by the scarcity of a systematic data set that correlates boundary structure to failure, especially in BCC metals. To fill in this gap in the current understanding, molecular dynamics simulations are performed on a set of 74 bi-crystals in Ta with a [110] symmetric tilt axis. Our results show a correlation between GB misorientation angle and spall strength and also highlight the importance of GB structure itself in determining the spall strength. Specifically, we find a direct correlation between the ability of the GB to plasticity deform through slip/twinning and its spall strength. Additionally, a change in the deformation mechanism from dislocation-meditated to twinning-dominated plasticity is observed as a function of misorientation angles, which results in lowered spall strengths for high-angle GBs.

Chen, J. (ORCID:000000018633229X)↗
DOE OSTI2020

Perspectives on the generation of electron beams from plasma-based accelerators and their near and long term applications

This article first gives the authors' perspectives on how the field of plasma-based acceleration (PBA) developed and how the current experiments, theory, and simulations are motivated by long term applications of PBA to a future linear collider and an x-ray free electron laser. Here, we then focus on some early applications that will likely emerge from PBA research such as electron beam radiotherapy, directional but incoherent x-ray beams for science and technology, near single cycle continuously tunable infrared pulses for spectroscopy, and non-perturbative quantum electrodynamics enabled by PBA electron beams. In our opinion, these near term applications could be developed within the next decade with a concerted effort by the community.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The influence of boundary and edge-plasma modeling in computations of axisymmetric vertical displacement

A number of previously published studies apply extended-MHD computations to model macroscopic dynamics of tokamak vertical displacement events (VDEs). The characteristic plasma-wall contact and resistive-wall diffusion imply sensitivity to boundary conditions in a general sense, but which conditions have significant influence depends on how the plasma is modeled. This work examines these dependencies by comparing results of axisymmetric extended-MHD computations with different sets of boundary conditions, plasma thermal-conduction models, and wall resistivity values. The geometry and plasma parameters of the computations represent a modest-sized tokamak. The forced-displacement scenario and computational setup are similar to those used in previous 3D computations. The results show that for a given wall resistivity, the VDE time-scale is most sensitive to variations in the boundary and thermal-conduction parameters that affect electron thermal transport. The electrical conductivity depends on electron temperature, and the dependence on thermal transport stems from its influence on the electrical circuit that includes the open-field halo current. Lastly, conditions that lead to hotter, broader halo regions slow the evolution. Significant sensitivity to the boundary condition on plasma flow-velocity exists when electron thermal conduction is restricted and electron energy loss is convective, which is expected for conditions at the entrance of the magnetic presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulation of an inductively coupled plasma with a two-dimensional Darwin particle-in-cell code

A two-dimensional particle-in-cell code for the simulation of low-frequency electromagnetic processes in laboratory plasmas has been developed. The code uses the Darwin method omitting the electromagnetic wave propagation. The Darwin method separates the electric field into solenoidal and irrotational parts. The solenoidal electric field is calculated with a new algorithm based on the equation for the electric field vorticity. The system of linear equations in the new algorithm is readily solved using a standard iterative method. The irrotational electric field is the electrostatic field calculated with the direct implicit algorithm. The code is verified by reproducing the two-stream instability, electron electromagnetic waves, and shear Alfvén waves. The code is applied to simulate an inductively coupled plasma with the driving current flowing around the plasma region. In this simulation, a ring of dense plasma forms at the initial stage but then the density becomes maximal in the center and decays monotonically toward the walls. The skin effect is in the transitional mode between local and non-local, and the electron velocity distribution function is non-Maxwellian.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Classification of dynamical Lie algebras generated by spin interactions on undirected graphs

Dynamical Lie algebras (DLAs) are a versatile tool for various topics that span from the expressibility-trainability of variational quantum algorithms (VQAs), to simulation of many body Hamiltonians. Quantum gates and most of the Hamiltonians of interest consist of local interactions; therefore, the analysis of all possible DLAs generated by 1- and 2-local operators is crucial for quantum simulation and VQAs on current hardware. Previously in [R. Wiersema et al ., npj Quantum Inf. 10 , 110 (2024)], we analyzed the DLAs on linear, circular and all-to-all topologies, and obtained results about their dimensions and algebraic structure. Here, in this work, we extend our analysis into any possible hardware topology and provide a classification of all DLAs generated by Pauli strings on any undirected interaction graph. Our results indicate that the DLAs depend solely on whether the connectivity or interaction graph is bipartite or not. In addition, we find that the non-trivial polynomially scaling DLAs appear only on 1D line or circle topologies, and all other DLAs have dimensions scaling exponentially with the system size. Together with the current VQA literature, our results imply that either the majority of VQAs are non-trainable, or we are yet to understand the role of DLAs on the trainability of VQAs.

Algebraic structures↗
DOE OSTI2020

A Bayesian analysis of nuclear deformation properties with Skyrme energy functionals

In spite of numerous scientific and practical applications, there is still no comprehensive theoretical description of the nuclear fission process based solely on protons, neutrons and their interactions. Furthermore, the most advanced simulations of fission are currently carried out within nuclear density functional theory (DFT). In spite of being fully quantum-mechanical and rooted in the theory of nuclear forces, DFT still depends on a dozen or so parameters characterizing the energy functional. Calibrating these parameters on experimental data results in uncertainties that must be quantified for applications. This task is very challenging because of the high computational cost of DFT calculations for fission. In this paper, we use Gaussian processes to build emulators of DFT models in order to quantify and propagate statistical uncertainties of theoretical predictions for a range of nuclear deformations relevant to describing the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

The Atacama Cosmology Telescope: DR6 power spectrum foreground model and validation

We discuss the model of astrophysical emission at millimeter wavelengths used to characterize foregrounds in the multi-frequency power spectra of the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6), expanding on Louis et al. (2025) (2503.14452). We detail several tests to validate the capability of the DR6 parametric foreground model to describe current observations and complex simulations, and show that cosmological parameter constraints are robust against model extensions and variations. We demonstrate consistency of the model with pre-DR6 ACT data and observations from Planck and the South Pole Telescope. We evaluate the implications of using different foreground templates and extending the model with new components and/or free parameters. In all scenarios, the DR6 ΛCDM and ΛCDM+N eff cosmological parameters shift by less than 0.5σ relative to the baseline constraints. Some foreground parameters shift more; we estimate their systematic uncertainties associated with modeling choices. From our constraint on the kinematic Sunyaev-Zel'dovich power, we obtain a conservative limit on the duration of reionization of Δz rei < 4.4, assuming a reionization midpoint consistent with optical depth measurements and a minimal low-redshift contribution, with varying assumptions for this component leading to tighter limits. Finally, we analyze realistic non-Gaussian, correlated microwave sky simulations containing Galactic and extragalactic foreground fields, built independently of the DR6 parametric foreground model. Processing these simulations through the DR6 power spectrum and likelihood pipeline, we recover the input cosmological parameters of the underlying cosmic microwave background field, a new demonstration for small-scale CMB analysis. These tests validate the robustness of the ACT DR6 foreground model and cosmological parameter constraints.

CMBR experiments↗
DOE OSTI2023

Magnetic island formation and rotation braking induced by low-Z impurity penetration in an EAST plasma

Abstract Recent observations of the successive formations of the 4 / 1 , 3 / 1 , and 2 / 1 magnetic islands as well as the subsequent braking of the 2 / 1 mode during a low- Z impurity penetration process in EAST experiments are well reproduced in our 3 D resistive MHD simulations. The enhanced parallel current perturbation induced by impurity radiation predominately contributes to the tearing mode growth, and the 2 / 1 island rotation is mainly damped by the impurity accumulation as results of the influence from high n modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗