Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sim”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Ultrahigh-energy tau neutrino cross sections with GRAND and POEMMA

Next-generation neutrino experiments will push the limits in our understanding of astroparticle physics in the neutrino sector to energies orders of magnitude higher than the current state-of-the-art high-energy neutrino experiment, IceCube. These experiments will use neutrinos to tell us about the most extreme environments in the Universe while simultaneously leveraging these extreme environments as probes of neutrino properties at the highest energies accessible in the foreseeable future: $E \sim 10^9$ GeV. At these energies neutrinos are readily absorbed in Earth. Assuming an isotropic distribution, by looking at how the flux varies as a function of the angle through Earth, we show here that it is possible to extract the $ν_τ$-N cross section with precision at the ~20% level (1σ assuming Wilks' theorem) given $N_{events} \sim 100$ events.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Gravitational Contact Interactions and the Physical Equivalence of Weyl Transformations in Effective Field Theory

Theories of scalars and gravity, with non-minimal interactions, $\sim (M_P^2 +F(\phi) )R +L(\phi)$, have graviton exchange induced contact terms. These terms arise in single particle reducible diagrams with vertices $\propto q^2$ that cancel the Feynman propagator denominator $1/q^2$ and are familiar in various other physical contexts. In gravity these lead to additional terms in the action such as $\sim F(\phi) T_\mu^\mu(\phi)/M_P^2$ and $F(\phi)\partial^2 F(\phi)/M_P^2$. The contact terms are equivalent to induced operators obtained by a Weyl transformation that removes the non-minimal interactions, leaving a minimal Einstein-Hilbert gravitational action. This demonstrates explicitly the equivalence of different representations of the action under Weyl transformations, both classically and quantum mechanically. To avoid such "hidden contact terms" one is compelled to go to the minimal Einstein-Hilbert representation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Utilizing cosmic-ray positron and electron observations to probe the averaged properties of Milky Way pulsars

Pulsars have long been studied in the electromagnetic spectrum. Their environments are rich in high-energy cosmic-ray electrons and positrons likely enriching the interstellar medium with such particles. In this work we use recent cosmic-ray observations from the AMS-02, CALET and DAMPE collaborations to study the averaged properties of the local Milky Way pulsar population. We perform simulations of the local Milky Way pulsar population, for interstellar medium assumptions in agreement with a range of cosmic-ray nuclei measurements. Each such simulation contains $\sim 10^{4}$ pulsars of unique age, location, initial spin-down power and cosmic-ray electron/positron spectra. We produce more than $7 \times 10^3$ such Milky Way pulsar simulations. Further, we account for and study i) the pulsars' birth rates and the stochastic nature of their birth, ii) their initial spin-down power distribution, iii) their time evolution in terms of their braking index and characteristic spin-down timescale, iv) the fraction of spin-down power going to cosmic-ray electrons and positrons and v) their propagation through the interstellar medium and the Heliosphere. We find that pulsars of ages $\sim 10^{5}-10^{7}$ yr, have a braking index that on average has to be 3 or larger. Given that electromagnetic spectrum observations of young pulsars find braking indices lower than 3, our work provides strong hints that pulsars' braking index increases on average as they age, allowing them to retain some of their rotational energy. Moreover, we find that pulsars have relatively uniform properties as sources of cosmic-ray electrons and positrons in terms of the spectra they produce and likely release $O(10\%)$ of their rotational energy to cosmic rays in the ISM. Finally, we find at $\simeq$12 GeV positrons a spectral feature that suggests a new subpopulation of positron sources contributing at these energies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Low-mass dark matter search with CDMSlite

The SuperCDMS experiment is designed to directly detect WIMPs (Weakly Interacting Massive Particles) that may constitute the dark matter in our galaxy. During its operation at the Soudan Underground Laboratory, germanium detectors were run in the CDMSlite (Cryogenic Dark Matter Search low ionization threshold experiment) mode to gather data sets with sensitivity specifically for WIMPs with masses ${<}10$ GeV/$c^2$. In this mode, a large detector-bias voltage is applied to amplify the phonon signals produced by drifting charges. This paper presents studies of the experimental noise and its effect on the achievable energy threshold, which is demonstrated to be as low as 56 eV$_{\text{ee}}$ (electron equivalent energy). The detector biasing configuration is described in detail, with analysis corrections for voltage variations to the level of a few percent. Detailed studies of the electric-field geometry, and the resulting successful development of a fiducial parameter, eliminate poorly measured events, yielding an energy resolution ranging from ${\sim}$9 eV$_{\text{ee}}$ at 0 keV to 101 eV$_{\text{ee}}$ at ${\sim}$10 keV$_{\text{ee}}$. New results are derived for astrophysical uncertainties relevant to the WIMP-search limits, specifically examining how they are affected by variations in the most probable WIMP velocity and the galactic escape velocity. These variations become more important for WIMP masses below 10 GeV/$c^2$. Finally, new limits on spin-dependent low-mass WIMP-nucleon interactions are derived, with new parameter space excluded for WIMP masses ${\lesssim}$3 GeV/$c^2$.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Preheating after multifield inflation with nonminimal couplings. II. Resonance structure

This is the second in a series of papers on preheating in inflationary models comprised of multiple scalar fields coupled nonminimally to gravity. We work in the rigid-spacetime approximation and consider field trajectories within the single-field attractor, which is a generic feature of these models. We construct the Floquet charts to find regions of parameter space in which particle production is efficient for both the adiabatic and isocurvature modes, and analyze the resonance structure using analytic and semianalytic techniques. Particle production in the adiabatic direction is characterized by the existence of an asymptotic scaling solution at large values of the nonminimal couplings, ${{\xi}}_{I}{\gg}1$, in which the dominant instability band arises in the long-wavelength limit, for comoving wave numbers $k{\rightarrow}0$. However, the large-${{\xi}}_{I}$ regime is not reached until ${{\xi}}_{I}{\ge}\mathcal{O}(100)$. In the intermediate regime, with ${{\xi}}_{I}{\sim}\mathcal{O}(1-10)$, the resonance structure depends strongly on wave number and couplings. The resonance structure for isocurvature perturbations is distinct and more complicated than its adiabatic counterpart. An intermediate regime, for ${{\xi}}_{I}{\sim}\mathcal{O}(1-10)$, is again evident. For large values of ${{\xi}}_{I}$, the Floquet chart consists of densely spaced, nearly parallel instability bands, suggesting a very efficient preheating behavior. The increased efficiency arises from features of the nontrivial field-space manifold in the Einstein frame, which itself arises from the fields’ nonminimal couplings in the Jordan frame, and has no analog in models with minimal couplings. Quantitatively, the approach to the large-${{\xi}}_{I}$ asymptotic solution for isocurvature modes is slower than in the case of the adiabatic modes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Diffusion-dominated mixing in moderate convergence implosions

High-$Z$ material mixed into the fuel degrades inertial fusion implosions and can prevent ignition. Mix is often assumed to be dominated by hydrodynamic instabilities, but in this paper we report Omega data, using shells with ${\sim}150\mathrm{nm}$ deuterated layers to gain unprecedented resolution, which give strong evidence that the dominant mix mechanism is diffusion for these moderate temperature (${\lesssim}6$ keV) and convergence (${\sim}12$) implosions. Finally, small-scale instability-driven or turbulent mix is negligible.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Highly Resolved Measurements of a Developing Strong Collisional Plasma Shock

In this paper, the structure of a strong collisional shock front forming in a plasma is directly probed for the first time in laser-driven gas-jet experiments. Thomson scattering of a 526.5 nm probe beam was used to diagnose temperature and ion velocity distribution in a strong shock ($M{\sim}11$) propagating through a low-density (${\rho}{\sim}0.01\text{ }\text{ }\mathrm{mg}/\mathrm{cc}$) plasma composed of hydrogen. A forward-streaming population of ions traveling in excess of the shock velocity was observed to heat and slow down on an unmoving, unshocked population of cold protons, until ultimately the populations merge and begin to thermalize. Finally, instabilities are observed during the merging, indicating a uniquely plasma-phase process in shock front formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Kinetic Theory of Spin Diffusion and Superdiffusion in X X Z Spin Chains

We address the nature of spin transport in the integrable XXZ spin chain, focusing on the isotropic Heisenberg limit. We calculate the diffusion constant using a kinetic picture based on generalized hydrodynamics combined with Gaussian fluctuations: we find that it diverges, and show that a self-consistent treatment of this divergence gives superdiffusion, with an effective time-dependent diffusion constant that scales as $D(t) \sim t^{1/3}$. This exponent had previously been observed in large-scale numerical simulations, but had not been theoretically explained. We briefly discuss XXZ models with easy-axis anisotropy $\Delta > 1$. Our method gives closed-form expressions for the diffusion constant $D$ in the infinite-temperature limit for all $\Delta > 1$. We find that $D$ saturates at large anisotropy, and diverges as the Heisenberg limit is approached, as $D \sim (\Delta - 1)^{-1/2}$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Quantification of Dopant Profiles in SiGe HBT Devices

We report on the use of atom probe tomography (APT), scanning transmission electron microscopy (STEM), and secondary ion mass spectroscopy (SIMS) to characterize doping profiles in the base region of SiGe HBT devices. We compare SIMS profiles obtained from large regions (400 um 2 ) of the device wafer to profiles obtained from individual devices of different emitter window widths (0.25 and 0.18 um 2 ) using APT. From this comparison we show how APT can provide a deeper insight into evaluating the fabrication process and its effects on electrical models of device performance and enabling the building of higher performance systems. We also demonstrate that APT can be used to characterize defects within the intrinsic regions of a device.

Jones, Eric↗
DOE OSTI2018

Internal Resistive Barriers Related to Zinc Diffusion During the Growth of Inverted Metamorphic Multijunction Solar Cells

Majority carrier barriers at heterointerfaces are a common source of non-linear resistance that hinders concentrator solar cell performance. The source of a particular barrier is often unclear in a multijunction device with numerous heterointerfaces. In this work, we demonstrate Zn-dopant diffusion during inverted metamorphic multijunction (IMM) device growth to be one key cause of internal barrier formation. Using an inverted GaAs/GaAs tandem solar cell with a high temperature annealing layer grown in between each subcell, we simulate the annealing conditions of a multijunction growth in a simplified structure. Through analysis of the device by secondary ion mass spectrometry (SIMS) and electrochemical capacitance-voltage profiling, we show that annealing causes Zn to diffuse out of the top cell Ga 0.5 In 0.5 P back surface field (BSF) and accumulate in the GaAs base. Through equilibrium band modeling, we show that the resultant doping profile forms an energetic barrier to hole flow in the valence band, which correlates with fill factor losses in the current-voltage curves measured under concentration. When we, instead, employ a C-doped Al 0.2 Ga 0.8 As BSF layer in the top cell, we do not observe evidence of a heterojunction barrier. We attribute this difference to the reduced diffusivity of carbon, confirmed by SIMS, as well as more favorable valence band offsets between GaAs and Al 0.2 Ga 0.8 As. Finally, we compare 5-junction IMM cells with Al 0.2 Ga 0.8 As:C and Ga 0.5 In 0.5 P:Zn BSF layers in the GaAs third junction, respectively, and show a significantly improved concentrator device performance when Al 0.2 Ga 0.8 As:C is employed. Lastly, we demonstrate the importance of designing annealing tolerance into multijunction structures that are subjected to extended annealing during growth.

14 SOLAR ENERGY↗
DOE OSTI2019

PPT-GPU: Scalable GPU Performance Modeling

Performance modeling is a challenging problem due to the complexities of hardware architectures. In this paper, we present PPT-GPU, a scalable and accurate simulation framework that enables GPU code developers and architects to predict the performance of applications in a fast, and accurate manner on different GPU architectures. PPT-GPU is part of the open source project, Performance Prediction Toolkit (PPT) developed at the Los Alamos National Laboratory. We extend the old GPU model in PPT that predict the runtimes of computational physics codes to offer better prediction accuracy, for which, we add models for different memory hierarchies found in GPUs and latencies for different instructions. To further show the utility of PPT-GPU, we compare our model against real GPU device(s) and the widely used cycle-accurate simulator, GPGPU-Sim using different workloads from RODINIA and Parboil benchmarks. The results indicate that the predicted performance of PPT-GPU is within a 10% error compared to the real device(s). In addition, PPT-GPU is highly scalable, where it is up to 450x faster than GPGPU-Sim with more accurate results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Comparison of PID Shunting in Polycrystalline and Single-Crystal Silicon Modules via Multi-Scale, Multi-Technique Characterization

We used the methods we reported last year to investigate potential-induced degradation (PID). We have now applied these methods to single-crystalline silicon modules that have degraded during field deployment, as well as in minimodules stressed in the laboratory. We will compare these results to the polycrystalline results presented last year. Small cores have been removed from the modules and subjected to analysis. We use a combination of photoluminescence and dark lock-in thermography imaging, laser marking, electron-beam induced current measurements, and subsequent focused ion-beam marking to allow analysis of individual defects via time-of-flight secondary-ion mass spectrometry (TOF-SIMS) to investigate the root-cause mechanism for PID shunting. We see a direct correlation between recombination active shunts and sodium content. The sodium content in shunted areas peaks at the SiN/Si interface and is consistently observed at a concentration of 0.1%-1% in shunted areas. TOF-SIMS data taken on degraded and non-degraded single-crystalline sample areas show a similar trend as in the polycrystalline samples: more sodium is seen in the degraded areas.

14 SOLAR ENERGY↗
DOE OSTI2024

An Evaluation of the Effect of Network Cost Optimization for Leadership Class Supercomputers

Dragonfly-based networks are an extensively deployed network topology in large-scale high-performance computing due to their cost-effectiveness and efficiency. The US will soon have three Exascale supercomputers for leadership class workloads deployed using dragonfly networks. Compared to indirect networks of similar scale, the dragonfly network has considerably reduced cable lengths, cable counts, and switch counts, resulting in significant network cost savings for a given system size, however, these cost reductions result in reduced global minimal paths and more challenging routing. Additionally, large scale dragonfly networks often require a taper at the global link level, resulting in less bisection bandwidth than is achievable in other traditional non-blocking topologies of equivalent scale. While dragonfly networks have been extensively studied, they have yet to be fully evaluated in an extreme scale (i.e., exascale) system that targets capability workloads. In this paper, we present the results of the first large scale evaluation of a dragonfly network on an exascale system (Frontier) and compare its behavior to a similar scale fat-tree network on a previous generation TOP500 system (Summit). This evaluation aims to determine the effect of network cost optimizations by measuring a tapered topology’s impact on capability workloads. Our evaluation is based on a collection of synthetic microbenchmarks, mini-apps, and full scale applications. It compares the scaling efficiencies of each benchmark between the dragonfly-based Frontier and the fat-tree-based Summit systems. Our results show that a dragonfly network is $\sim \mathbf{3 0 \%}$ more cost efficient than a fat-tree topology, which amortizes to $\sim 3 \%$ of an exascale system cost. Furthermore, while tapered dragonfly networks impose significant tradeoffs, the impacts are not as broad as initially thought and are mostly seen in applications with global communication patterns, particularly all-to-all (e.g., FFT-based algorithms), but also local communication patterns (e.g., nearest-neighbor algorithms) that are sensitive to network performance variability.

Khan, Awais↗
DOE OSTI2023

Characterization of MKIDs for CMB Observation at 220 GHz With the South Pole Telescope

We present an updated design of the 220 GHz microwave kinetic inductance detector (MKID) pixel for SPT-3G+, the next-generation camera for the South Pole Telescope. We show results of the dark testing of a 63-pixel array with mean inductor quality factor $Q_{i} = 4.8 \times 10^{5}$, aluminum inductor transition temperature $T_{c} = 1.19$K, and kinetic inductance fraction $\alpha _{k} = 0.32$. We optically characterize both the microstrip-coupled and CPW-coupled resonators, and find both have a spectral response close to prediction with an optical efficiency of $\eta \sim 70\%$. However, we find slightly lower optical response on the lower edge of the band than predicted, with neighboring dark detectors showing more response in this region, though at level consistent with less than 5% frequency shift relative to the optical detectors. The detectors show polarized response consistent with expectations, with a cross-polar response of $\sim 10\%$for both detector orientations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Secondary Ion Mass Spectrometry Methodology for Isotopic Ratio Measurement of Micro‐Grains in Thin Sections: True Grain Size Estimation and Deconvolution of Inter‐Grain Size Gradients and Intra‐Grain Radial Gradients

In a thin section, grains that were approximately spherical in situ appear circular in cross section, and the distribution of apparent diameters frequently assumed to be their size distribution. Scanning ion imaging by secondary ion mass spectrometry (SIMS) is capable of providing precise (< 1‰) stable isotope ratio measurements of such grains, but, importantly, also registers their rate of evolution in apparent size as they are ablated by the primary beam. By assessing rates of radius change with depth, the described methodology enables the ‘true’ size of grains to be estimated, as well as the distance of the sectioned surface from the original grain centre. Transects in three dimensions are made possible, and this capability enables better identification (and thus separation) of both inter‐grain chemical signatures as a function of grain size, and intra‐grain radial trends. In this example, we highlight the specific application to pyrite (FeS 2 ) minerals, which are frequently analysed by SIMS to determine their inter‐grain and intra‐grain geochemical variations, particularly in their sulfur stable isotopic ratios (δ 34 S). Benefits of the new methodology over the Faraday cup ‘spot mode’ are described. Data correction algorithms and precision considerations are discussed.

Jones, Clive↗
DOE OSTI2019

Assessment of Long Term Creep Using Strain Rate Matching From the Stepped Isostress Method

Creep testing is an ongoing need, particularly with the development of new candidate alloy systems for advanced energy systems. The conventional creep test (CT) is regarded as a proven method to gather creep data however, the test is impractical due to being real-time: lasting up to 105 hours to characterize the service of long-lived turbomachinery components. Accelerated methods to gather the long-term creep properties of materials are needed to reduce the time to qualification of new materials. The time-temperature-stress-superposition principle (TTSSP) and the derivative time-temperature superposition principle (TTSP), time-stress superposition principle (TSSP), stepped isothermal method (SIM), and stepped isostress method (SSM) are accelerated creep tests (ACT) commonly used to predict the long-term creep behaviors of polymers and composites. The TTSP and TSSP tests require multiple specimen tested at various temperatures/stresses whereas the SIM and SSM tests employ a single specimen where temperature/stress are periodically step increased until rupture. The stepped creep deformation curve can then be time and strain shifted to produce a master creep curve. While these ACTs are useful tools to predict long-term creep, the drawback is the lack of mathematical laws to determine the virtual start time and time shift factors, especially for different materials. In this paper, a new self-calibration approach is developed and compared to existing SSM data for Kevlar 49. This new approach focuses on matching the creep strain rates between stress steps and fitting the data to a master curve using a modified theta projection model. This is performed using a MATLAB code consisting of five subroutines. The first subroutine takes the stress, time, and creep strain from SSM/SIM tests, and segregates the data intro arrays corresponding to each stress level. The second subroutine finds the constants for the modified theta projection model for each stress level. The third subroutine performs a time shift adjustment using creep strain rate matching. The fourth subroutine calculates the accelerated time of rupture. The last subroutine generates accelerated creep versus time plots. Kevlar 49 SSM data is gathered from literature and run through the MATLAB code. The master curves generated from the MATLAB are compared to the conventional creep curve of Kevlar 49 as well as the master curve gathered from literature in order to validate the feasibility of this new approach. The goal of this project is to vet if the self-calibration approach can produce results similar to the reference calibration approach.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Analysis of a Wave-Powered, Reverse-Osmosis System and Its Economic Availability in the United States: Preprint

A wave energy converter (WEC) system has the potential to convert the wave energy resource directly into the high-pressure flow that is needed by the desalination system to permeate saltwater through the reverse-osmosis membrane to generate clean water. In this study, a wave-to-water numerical model was developed to investigate the potential use of a wave-powered desalination system (WPDS) for water production in the United States. The model was developed by coupling a time-domain radiation-and-diffraction-method-based numerical tool (WEC-Sim) for predicting the hydrodynamic performance of WECs with a solution-diffusion model that was used to simulate the reverse-osmosis process. To evaluate the feasibility of the WPDS, the wave-to-water numerical model was applied to simulate a desalination system that used an oscillating surge WEC device to pump seawater through the system. The annual water production was estimated based on the wave resource at a reference site on the coast of northern California to investigate the potential cost of water in that area, where the cost of water and electricity is high compared to other regions. In the scenario evaluated, for a 100-unit utility-scale electricity-producing array, the estimated levelized cost of energy for these WECs is about 3-6 times the U.S.'s current, unsubsidized electricity rates. However, with clean water as an end product and by directly producing pressurized water with WECs, rather than electricity as an intermediary, it is presently only 12 percent greater than typical water cost in California. This study suggests that a WEC array that produces water may be a viable, near-term solution to the nation's water supply, and the niche application of the WPDS may also provide developers with new opportunities to further develop technologies that benefit both the electric and drinking water markets.

16 TIDAL AND WAVE POWER↗
DOE OSTI2018

Numerical Analysis on Hydraulic Power Take-Off for Wave Energy Converter and Power Smoothing Methods: Preprint

One of the primary challenges for wave energy converter (WEC) systems is the fluctuating nature of wave resources, which require the WEC components to be designed to handle loads (i.e., torques, forces, and powers) that are many times greater than the average load. This approach requires a much greater power take-off (PTO) capacity than the average power output and indicates a higher cost for the PTO. Moreover, additional design requirements, such as battery storage, are needed, particularly for practical electrical grid connection, and can be a problem for sensitive equipment (e.g., radar, computing devices, and sensors). Therefore, it is essential to investigate potential methodologies to reduce the overall power fluctuation while trying to optimize the power output from WECs. In this study, a detailed hydraulic PTO model was developed and coupled with a time-domain hydrodynamics model (WEC-Sim) to evaluate the PTO efficiency for WECs and the trade-off between power output and fluctuation using different power smoothing methods, including energy storage, pressure relief mechanism, and a powerbased setpoint control method. The study also revealed that the maximum power fluctuation for WECs can be significantly reduced by one order of magnitude when these power smoothing methods are applied.

16 TIDAL AND WAVE POWER↗