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 415 records · Page 23

Charge carrier motion and effect of fixed oxide charge in a microstructured silicon radiation detector

Signal formation in a microstructured semiconductor neutron detector is more complex than in planar diode geometry. Three-dimensional microstructures are laterally smaller than the ionization cloud length, and the electric fields may be weak enough to exhibit plasma time effects. This work is the first detailed treatment of charge carrier motion in these complex semiconductor devices to replicate the time profile and signal magnitude. Simulations were performed using COMSOL Multiphysics to investigate various parameters that affect the propagation of the charge cloud. It was observed that the size of the simulated three-dimensional structure had an impact on the induced current pulse, indicating the importance of simulation geometry optimization to accurately simulate charge cloud expansion. COMSOL Multiphysics was used to replicate accurate charge creation profiles using energy deposition information imported from radiation transport codes. A detailed simulation methodology is presented to benchmark preamplifier event pulses along with complexities in modeling the charge carrier motion along the etched microstructured trenches with Si–SiO 2 boundary conditions, including fixed oxide charge and interface trapping.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Experimental Determination of PV Inverter Response to Grid Phase Shift Events: Preprint

With the continued growth of renewable energy resources which interface to the electric grid via inverters, the understanding of such devices becomes ever more important to the safe and reliable operation of the bulk power system. This work investigates the specific response of a utility-scale PV inverter to grid voltage phase shift-type disturbances which sometimes occur during grid fault events. The role of the PV inverter's phase-locked-loop (PLL) is identified as important to modeling the response. Switching-level simulations of a utility-scale PV inverter with a modeled PLL show a characteristic response when phase shift disturbances require the PLL to track what appear as fast frequency changes. Additionally, in this work a full-scale laboratory testing was carried out using the Opal real time (RT) OP5142 real time simulator and a large grid simulator to create phase shift disturbances with a high degree of repeatability. A photovoltaic (PV) inverter was connected to a grid simulator, and phase shifts were instantaneously implemented on the simulated grid, the results of the currents were then obtained. The experimental results obtained were validated with simulation results obtained from MATLAB/Simulink.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Basic factors for acquiring, correcting, and interpreting probe current-voltage characteristic in moderate-collisional plasma for determining electron energy distribution

Methods for inferring the electron distribution function (EDF) from a probe current–voltage trace are reviewed for moderate-collisional plasma. The explicit accounting of particle collisions, modeling sheath depletion, and electron refilling processes are rarely performed when analyzing probe measurements, yet they are critical for regime validity, EDF-profile accuracy, and space potential quantification precision. Furthermore, the advantages, limitations, and utilization of large-area wall probes in moderate-collisional plasma, capable of highly resolved indications of energetic electrons in plasma, are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

Understanding and predicting divertor heat-load width λq is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for λ q using the extreme-scale edge gyrokinetic code XGC1 [S. Ku et al., Phys. Plasmas 25, 056107 (2018)] in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C. S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013) and R. J. Goldston, Nucl. Fusion 52, 013009 (2012)] and furthermore predicted over six times wider λ q than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for λ q as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide λ q XGC for the full-current Q = 10 ITER plasma. Additionally, the new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider λ q XGC is identified to be the weakly collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stochastic parametric skeletal dosimetry model for humans: Pediatric and adult computational skeleton phantoms for internal bone marrow dosimetry

Currently, computational phantoms that simulate skeletal tissues are used in active red bone marrow (AM) internal dosimetry. Up-to-date reference computational phantoms recommended by the ICRP are based on the analysis of CT-images of cadavers. Such phantoms have significant disadvantages. One disadvantage is that the assessment of uncertainty due to the population variability of skeleton dimensions and microstructure results from the limited availability of autopsy material. Another disadvantage is the simplified modelling of cortical layer and bone microarchitecture. A method of stochastic parametric skeletal dosimetry modelling of the bone structures – SPSD modelling – has been developed as an alternative to the ICRP reference phantoms. In the framework of this approach, skeletal phantom parameters are evaluated based on extensively reviewed results of published measurements of real bones. The SPSD approach allows for the assessment of both population-average values and their variability. SPSD-phantoms of the skeleton are modelled in voxel representation. They consist of smaller phantoms of the bone sites – segments – described by simple geometric shapes with uniform microarchitecture parameters. Such segmentation makes it possible to account for non-homogeneous skeletal microarchitecture and to model the bone structure with the required voxel resolution to elaborate suitable skeletal phantoms. The current study presents the parameters of the SPSD skeletal phantoms for the following age-groups: newborn, 1-year-old, 5-year-old, 10-year-old, 15-year-old (male and female), and adult (male and female). This skeletal phantom can be used for dosimetry as an alternative to available reference phantoms for bone-seeking radionuclides. The above-mentioned age- and sex-specific skeletal phantoms are comprised of 289 unique segments. The characteristics of the SPSD phantoms do not contradict published data and are in good agreement with the measurement results of real bones.

Science & Technology - Other Topics↗

PSTN-054: Updated estimates of the Rubin system throughput and expected LSST image depth

This document presents updated estimates, as of May 2022, of the Rubin system throughput and compares them to throughput requirements from the LSST Science Requirements Document. In addition, it uses these estimates to forecast LSST's median single-visit and co-added image depths for the current (OpSim v2.0) baseline LSST cadence simulation: (23.8, 24.5, 24.0, 23.4, 22.7, 22.0) and (25.6, 26.9, 26.9, 26.4, 25.6, 24.8) in ugrizy, respectively. In addition, it uses these estimates to forecast LSST’s median single-visit and co-added image depths for the current baseline LSST cadence simulation. Estimated system performance relies on actual measurements of the performance of various system hardware components, and on simulations where measurements are still unavailable. The updated performance estimates meet all the relevant requirements from the LSST Science Requirements Document. As our knowledge of the as-buit system continues to improve, updates to these estimates will continue to be provided.

79 ASTRONOMY AND ASTROPHYSICS↗

Kilonova Detectability with Wide-field Instruments

Abstract Kilonovae are ultraviolet, optical, and infrared transients powered by the radioactive decay of heavy elements following a neutron star merger. Joint observations of kilonovae and gravitational waves can offer key constraints on the source of Galactic r -process enrichment, among other astrophysical topics. However, robust constraints on heavy element production require rapid kilonova detection (within ∼1 day of merger) as well as multiwavelength observations across multiple epochs. In this study, we quantify the ability of 13 wide-field-of-view instruments to detect kilonovae, leveraging a large grid of over 900 radiative transfer simulations with 54 viewing angles per simulation. We consider both current and upcoming instruments, collectively spanning the full kilonova spectrum. The Roman Space Telescope has the highest redshift reach of any instrument in the study, observing kilonovae out to z ∼ 1 within the first day post-merger. We demonstrate that BlackGEM, DECam, GOTO, the Vera C. Rubin Observatory’s LSST, ULTRASAT, VISTA, and WINTER can observe some kilonovae out to z ∼ 0.1 (∼475 Mpc), while DDOTI, MeerLICHT, PRIME, Swift/UVOT, and ZTF are confined to more nearby observations. Furthermore, we provide a framework to infer kilonova ejecta properties following nondetections and explore variation in detectability with these ejecta parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

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↗

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↗

A novel optical polarization conversion structure designed by multi-layered thin film structures at the wavelength of 10.6 μm

Polarization conversion performance is one of the most important design contents of optical devices. At present, there are many ways to implement polarization conversion. However, they are mostly based on the form of split ring resonator or equivalent plasma, which are not conducive to the miniaturization of electromagnetic devices, particularly in the field of laser applications. In this paper, a novel optical polarization conversion structure has been designed by multi-layered PS/PTFE thin film structures at the wavelength of 10.6 μm based on the equivalent medium theory. Then it was verified that the proposed structure can effectively convert the linear polarized wave to circular polarized wave by the EM full-wave simulation. Compared with the current laser polarization manipulation devices, the proposed novel polarization converter has advantages of ultrathin thickness of structure, simple structure, good light transmission and no electromagnetic loss. Although the structural parameters in this paper are designed for EM wave the wavelength of 10.6 μm, the proposed design idea of the polarization conversion structure is also suitable for the electromagnetic wave at other wavelengths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗