SEARCH · Engineering Papers
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.
There are versions of QUIC for 64-bit Windows, 64-bit Linux, and 64-bit Intel Mac OS X. As compared to earlier 32-bit versions of QUIC, the 64-bit versions allow the QUIC transport and dispersion codes to access greater than 2 GB of RAM and thus larger problems can be run.
A new technology called Shear Assisted Processing and Extrusion (ShAPE™), pioneered at the Pacific Northwest National Laboratory (PNNL), has the potential to revolutionize the US extrusion industry. Maturation of the ShAPE™ process would create a new cross cutting US manufacturing technology that advances three key the objectives within AMO; namely, nanomaterials processing, process intensification, and materials for extreme and harsh environments. PNNL and SCM Metals Products, Inc. (SCM), a division of Kymera International, will collaborate to mature the ShAPE™ extrusion process beyond TRL 3 to increase the likelihood of its adoption by industry. PNNL’s custom (one-of-a-kind in the world) ShAPE™ machine will be used to extrude round rods of high strength aluminum alloy feedstock powder provided by SCM. As proof-of-concept, PNNL has previously extruded 5 mm diameter rods with SCM aluminum alloy powder AL-12.4TM. This project aims to scale up the process to a diameter of 25 mm while retaining the excellent mechanical properties achieved for the 5 mm diameter rods. The degree of process intensification with ShAPE™ will also be estimated. The primary objectives of this project are as follows. 1) Fabricate 25 mm diameter nanostructured extrusions in bulk, directly from AL-12.4TM aluminum alloy powder without the need for canning, degassing, compaction, and other intermediate steps typical of powder metallurgy extrusion. 2) Achieve a microstructure with an ultrafine grained (UFG <1 µm) aluminum matrix having refined and homogenously distributed nanoscale second phases. 3) Show that the ram force can be significantly reduced compared to conventional powder metallurgy extrusion of the same material and estimate the extent of process intensification that may be possible with ShAPE™. 4) Demonstrate improved high-temperature mechanical properties compared to conventional extrusion of the same material. The following tasks will be undertaken to achieve these objectives. PNNL will develop the die set and process parameters necessary to fabricate Al-12.4TM extrusions with a 25 mm diameter. Microstructural analysis will be performed by PNNL while mechanical testing will be conducted by SCM. With information learned from the process development and characterization, PNNL will perform an estimate of the degree of process intensification that may be possible with ShAPE™.
The degradation of antenna performance during hypersonic re-entry is a well known phenomenon that can lead to complete radio blackout. Recent additions to the Empire code establish it as a tool for the study and analysis of the problem. Coupling to the Sandia Parallel Aerodynamics and Reentry Code (SPARC) enables the electromagnetic analysis of realistic re-entry plasma profiles. The geometric flexibility afforded by both Empire and SPARC allow the consideration of arbitrary vehicle and antenna configurations. We have used this tool to study antenna performance during re-entry when the boundary layer becomes turbulent. A concise description of line-of-sight transmissions, which employs advanced statistical methods, was developed. New insights into the low altitude reflectometer readings of RAM-C2 are offered. Techniques for the reconstruction of the re-entry plasma profile from reflectometer data were explored.
Radioactive Material (RAM) Packagings Used to transport radioactive material • Weapons components • Medical isotopes • Spent nuclear fuel • Etc. • Packages must adhere to 10 CFR 71 which drive to NUREGs, ASME codes, etc. • Transport Index (A2 values) determine if a package is Type A, Type AF, Type B, Type B(U)F
Recent advances in probabilistic programming languages (PPLs) have provided the capability for exact inference: computing a closed-form probability distribution for a given probabilistic program. In particular, the new language Roulette uses a language oriented programming (LOP) approach, wherein analysts build new programming languages on top of a set of primitives provided by Roulette, which then translates these structures into a weighted model counting problem which can be solved by automated reasoning tools. However, because Roulette provides few convenience features, developing these new languages is challenging even for expert users. We developed a standard library of common probability functions for Roulette with the goal of improved usability. This included approximation of continuous probability density functions using discrete probability mass functions. We demonstrated this approach by modeling a cosmic ray striking a RAM controller. We found that Roulette provides a powerful interface for highly expressive probabilistic programs to be generated. In collaboration with the NNSA Advanced Simulation and Computing program, which resulted in development of a tool called Circulette, we were able to model complex circuits expressed in Verilog using probabilistic programs with an expressivity not previously possible. Our research question that motivated the development of a Roulette standard library was to determine whether non-experts could use a PPL to model relevant problems regarding radiation effects on microelectronics. This standard library improved the expressivity of Roulette by implementing common probability density functions, mathematical operators on distributions, and support for empirical distributions. While Roulette is a powerful modeling language, the untyped, LOP approach makes error messages difficult to understand and requires expert aid. We recommend further research on Roulette, especially with its error messages, to enable improved usability. At the same time, this project demonstrated that for users familiar with Roulette and the LOP approach, Roulette provides powerful new capabilities that can be integrated with other Sandia modeling capabilities.
In Fermilab's PIP-II machine protection system, beam loss signals from various detectors are digitized at 125 MS/s. Noise from both high-frequency sources and low-frequency 60 Hz AC power equipment can con-taminate the data. To suppress noise across these ranges especially 60 Hz and its harmonics, which overlap with beam loss signal frequencies advanced digital processing beyond standard filtering is re-quired. Several real-time functional blocks were simu-lated and tested on an FPGA: (1) a dual time-constant discharging integrator filter, (2) a de-ripple baseline extraction and storage block, and (3) a fast-recovery discharging integrator. The nonlinear IIR integrator filter removes high-frequency noise and feeds into the baseline extractor. Upon detecting abrupt beam loss, it switches to a longer time constant to prevent baseline distortion. The de-ripple block calculates a valid base-line by averaging over multiple 60 Hz periods, storing results in a 4096-word FPGA RAM. This baseline is subtracted from raw data before integration by the fast-recovery block, which resets quickly after use. All blocks achieved expected performance.
Ensuring powder mixture homogeneity is paramount for modern powder processing techniques. A promising recent development in the field of powder processing, resonant acoustic mixing (RAM), proposes to use low-frequency vertical vibration to deliver well-mixed powders without relying on mechanical agitators. However, an undesirable outcome of vibrations applied to granular materials comprised of non-uniform particle sizes is the tendency towards size segregation. In this work, we experimentally ex
Electromagnetic ion cyclotron (EMIC) waves are widely believed to play an important role in influencing the radiation belt and ring current dynamics. Most studies have investigated the effects or characteristics of EMIC waves by assuming their left-handed polarization. However, recent studies have found that the reversal of polarization, which occurs at higher latitudes along the wave propagation path, can change the wave-induced pitch angle diffusion coefficients. Whether such a polarization reversal can influence the global ring current dynamics remains unknown. In this study, we investigate the ring current dynamics and proton precipitation loss in association with polarization-reversed EMIC waves by using the ring current–atmosphere interactions model (RAM). The results indicate that the polarization reversal of H-band EMIC waves can truly decrease the scattering rates of protons of 10 to 50 keV or >100 keV in comparison with the scenario in which the EMIC waves are considered purely left-handed polarized. Additionally, the global ring current intensity and proton precipitation may be slightly affected by the polarization reversal, especially during prestorm time and the recovery phase, but the effects are not large during the main phase. This is probably because the H-band EMIC waves contribute to the proton scattering loss primarily at E < 10 keV, an energy range that is not strongly affected by the polarization reversal.
A pin style bottom anode employs conductive steel rods that serve as the pathway for the high electrical power through rammed refractory at the bottom of a DC Electric Arc Furnace (EAF). Anode wear during operation is important to monitor, as anode replacement is expensive and impacts EAF productivity. Liquid steel penetration into the un-sintered refractory layer can result from rapid electrical power ramp-up, dips in furnace temperature, or operating the anode for too long between EAF campaigns. In extreme cases, the liquid steel may penetrate the bottom of the furnace when anode wear progresses too close to the bottom shell, which is extremely dangerous and must be avoided. The current state of the art for monitoring bottom anode wear employs thermocouples imbedded in the anode pins at points in the anode. However, this approach is not sensitive enough to detect localized damage to the anode, especially when cracking occurs. Here, the present work utilizes fiber optic sensors to monitor the health of the anode, by creating a real-time spatially distributed temperature map of the anode. Unlike the traditional thermocouples, these sensors can be mounted at significantly greater depths, provide distributed temperature measurements, and can withstand temperatures of up to 900°C. Additionally, they are able to perform temperature measurements with a spatial resolution of 1.3 mm at a 5 Hz acquisition rate, providing unprecedented high-density real time monitoring of anode health and increasing the efficiency, and safety of EAF operation.
The collection and analysis of volatile memory is a vibrant area of research in the cybersecurity community. The ever-evolving and growing threat landscape is trending towards fileless malware, which avoids traditional detection but can be found by examining a system’s random access memory (RAM). Additionally, volatile memory analysis offers great insight into other malicious vectors. It contains fragments of encrypted files’ contents, as well as lists of running processes, imported modules, and network connections, all of which are difficult or impossible to extract from the file system. For these compelling reasons, recent research efforts have focused on the collection of memory snapshots and methods to analyze them for the presence of malware. However, to the best of our knowledge, no current reviews or surveys exist that systematize the research on both memory acquisition and analysis. We fill that gap with this novel survey by exploring the state-of-the-art tools and techniques for volatile memory acquisition and analysis for malware identification. For memory acquisition methods, we explore the trade-offs many techniques make between snapshot quality, performance overhead, and security. For memory analysis, we examined the traditional forensic methods used, including signature-based methods, dynamic methods performed in a sandbox environment, as well as machine learning-based approaches. We summarize the currently available tools, and suggest areas for more research.
We present a novel machine learning (ML)-based method to accelerate conservative-to-primitive inversion, focusing on hybrid piecewise polytropic and tabulated equations of state. Traditional root-finding techniques are computationally expensive, particularly for large-scale relativistic hydrodynamics simulations. To address this, we employ feedforward neural networks (NNC2PS and NNC2PL), trained in PyTorch (2.0+) and optimized for GPU inference using NVIDIA TensorRT (8.4.1), achieving significant speedups with minimal accuracy loss. The NNC2PS model achieves 𝐿 1 and 𝐿 ∞ errors of 4.54 × 10 −7 and 3.44 × 10−6, respectively, while the NNC2PL model exhibits even lower error values. TensorRT optimization with mixed-precision deployment substantially accelerates performance compared to traditional root-finding methods. Specifically, the mixed-precision TensorRT engine for NNC2PS achieves inference speeds approximately 400 times faster than a traditional single-threaded CPU implementation for a dataset size of 1,000,000 points. Ideal parallelization across an entire compute node in the Delta supercomputer (dual AMD 64-core 2.45 GHz Milan processors and 8 NVIDIA A100 GPUs with 40 GB HBM2 RAM and NVLink) predicts a 25-fold speedup for TensorRT over an optimally parallelized numerical method when processing 8 million data points. Moreover, the ML method exhibits sub-linear scaling with increasing dataset sizes. We release the scientific software developed, enabling further validation and extension of our findings. By exploiting the underlying symmetries within the equation of state, these findings highlight the potential of ML, combined with GPU optimization and model quantization, to accelerate conservative-to-primitive inversion in relativistic hydrodynamics simulations.
Although many galaxies in the Virgo cluster are known to have lost significant amounts of H i gas, only about a dozen features are known where the H i extends significantly outside its parent galaxy. Previous numerical simulations have predicted that H i removed by ram pressure stripping should have column densities far in excess of the sensitivity limits of observational surveys. We construct a simple model to try and quantify how many streams we might expect to detect. This accounts for the expected random orientation of the streams in position and velocity space as well as the expected stream length and mass of stripped H i. Using archival data from the Arecibo Galaxy Environment Survey, we search for any streams that might previously have been missed in earlier analyses. We report the confident detection of 10 streams as well as 16 other less-certain detections. We show that these well match our analytic predictions for which galaxies should be actively losing gas; however, the mass of the streams is typically far below the amount of missing H i in their parent galaxies, implying that a phase change and/or dispersal renders the gas undetectable. By estimating the orbital timescales, we estimate that dissolution rates of 1–10 M {sub ⊙} yr{sup −1} are able to explain both the presence of a few long, massive streams and the greater number of shorter, less-massive features.
We present the first in-depth dynamical analysis of the archetypal wide-angle tailed (WAT) cluster A562. We have combined Gemini observations with archival data from the literature to form a sample of 76 cluster members and derived a mean redshift of 0.1088 ± 0.0004 and a velocity dispersion of 919 ± 116 km s{sup −1}. This relatively large velocity dispersion suggests either a very massive cluster (M {sub dyn} > 6.9 × 10{sup 14} M {sub ☉}) and/or a merger system. The merger model is supported by a non-Gaussian galaxy velocity distribution, an elongated spatial distribution of likely cluster members, and an elongated X-ray emitting gas. This scenario would generate the bulk flow motion of the intracluster medium that can exert enough ram pressure to bend the radio jets. Thus, our observations support the model in which a recent off-axis merger event produced the cluster wide conditions needed to shape the WAT in A562.
In this work, we explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)--mass systems. We consider a generic self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, w, sets the velocity scale above which momentum transfer due to DM self-interactions becomes inefficient. We perform high-resolution zoom-in simulations of an MW-mass halo for values of w that span scenarios in which self-interactions either between the host and its subhalos or only within subhalos efficiently transfer momentum, and we study the effects of self-interactions on the host halo and on the abundance, radial distribution, orbital dynamics, and density profiles of subhalos in each case. The abundance and properties of surviving subhalos are consistent with being determined primarily by subhalo--host halo interactions. In particular, subhalos on radial orbits in models with larger values of the cross section at the host halo velocity scale are more susceptible to tidal disruption owing to mass loss from ram pressure stripping caused by self-interactions with the host. This mechanism suppresses the abundance of surviving subhalos relative to collisionless DM simulations, with stronger suppression for larger values of w. Thus, probes of subhalo abundance around MW-mass hosts can be used to place upper limits on the self-interaction cross section at velocity scales of ~200 km s-1, and combining these measurements with the orbital properties and internal dynamics of subhalos may break degeneracies among velocity-dependent SIDM models.
Warped disk galaxies are classified into two morphologies: S and U types. Conventional theories routinely attribute both types to galactic tidal interaction and/or gas accretion, but reproducing U types in simulations is extremely challenging. Here we investigate whether both types are governed by the same mechanisms using the most extensive sample of ~8000 nearby (0.02 < z < 0.06) massive (M*/M ⊙ > 10 9 ) edge-on disks from the Sloan Digital Sky Survey. We find that U types show on average bluer optical colors and a higher specific star formation rate (sSFR) than S types, with more strongly warped U types having a higher sSFR. We also find that while the S-type warp properties correlate with the tidal force by the nearest neighbor regardless of the environment, there is no such correlation for U types in groups/clusters, suggesting a nontidal mechanism could be at play for U types, such as ram pressure stripping (RPS). Indeed, U types are more common in groups/clusters than in fields and they have stellar mass, gas fraction, sSFR enhancement, and phase-space distribution closely analogous to RPS-induced jellyfish galaxies in clusters. We furthermore show that the stellar disks of most RPS galaxies in the IllustrisTNG simulation are warped in a U shape and bent in the opposite direction of stripped gas tails, satisfying theoretical expectations for stellar warps embedded in jellyfishes. We therefore suggest that despite the majority of U types that live in fields being still less explained, RPS can be an alternative origin for those in groups/clusters.
The circumgalactic medium (CGM) is often assumed to exist in or near hydrostatic equilibrium, with the regulation of accretion and the effects of feedback treated as perturbations to a stable balance between gravity and thermal pressure. We investigate global hydrostatic equilibrium in the CGM using four highly resolved L* galaxies from the Figuring Out Gas & Galaxies in Enzo (FOGGIE) project. The FOGGIE simulations were specifically targeted at fine spatial and mass resolution in the CGM (Δx ≲ 1 kpc h –1 and M ≃ 200M ⊙ ). We develop a new analysis framework that calculates the forces provided by thermal pressure gradients, turbulent pressure gradients, ram pressure gradients of large-scale radial bulk flows, centrifugal rotation, and gravity acting on the gas in the CGM. Thermal and turbulent pressure gradients vary strongly on scales of ≲5 kpc throughout the CGM. Thermal pressure gradients provide the main supporting force only beyond ~0.25R 200 , or ~50 kpc at z = 0. Within ~0.25R 200 , turbulent pressure gradients and rotational support provide stronger forces than thermal pressure. More generally, we find that global equilibrium models are neither appropriate nor predictive for the small scales probed by absorption line observations of the CGM. Local conditions generally cannot be derived by assuming a global equilibrium, but an emergent global equilibrium balancing radially inward and outward forces is obtained when averaging over the nonequilibrium local conditions on large scales in space and time. Approximate hydrostatic equilibrium holds only at large distances from galaxies, even when averaging out small-scale variations.
The engulfment of substellar bodies (SBs), such as brown dwarfs and planets, by giant stars is a possible explanation for rapidly rotating giants, lithium-rich giants, and the presence of SBs in close orbits around subdwarfs and white dwarfs. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an engulfed SB. We model the SB as a rigid body with a reflective surface because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram-pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory inside the star. We find that engulfment can increase the luminosity of a 1 M ⊙ star by up to a few orders of magnitude. The time for the star to return to its original luminosity is up to a few thousand years when the star has evolved to ≈10 R ⊙ and up to a few decades at the tip of the red giant branch (RGB). No SBs can eject the envelope of a 1 M ⊙ star before it evolves to ≈10 R ⊙ if the orbit of the SB is the only energy source contributing to the ejection. In contrast, SBs as small as ≈10 MJup can eject the envelope at the tip of the RGB. The numerical framework we introduce here can be used to study planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB.
The environments where galaxies reside crucially shape their star formation histories. We investigate a large sample of 1626 cluster galaxies located within 105 galaxy clusters spanning a large range in redshift (0.26 < z < 1.13). The galaxy clusters are massive (M 500 ≳ 2 × 10 14 M ⊙ ) and uniformly selected from the SPT and ACT Sunyaev–Zel'dovich surveys. With spectra in hand for thousands of cluster members, we use the galaxies' position in projected phase space as a proxy for their infall times, which provides a more robust measurement of environment than quantities such as projected clustercentric radius. We find clear evidence for a gradual age increase of the galaxy's mean stellar populations (~0.71 ± 0.4 Gyr based on a 4000 Å break, D n 4000) with the time spent in the cluster environment. This environmental quenching effect is found regardless of galaxy luminosity (faint or bright) and redshift (low or high-z), although the exact stellar age of galaxies depends on both parameters at fixed environmental effects. Such a systematic increase of D n 4000 with infall proxy would suggest that galaxies that were accreted into hosts earlier were quenched earlier due to longer exposure to environmental effects such as ram pressure stripping and starvation. Compared to the typical dynamical timescales of 1–3 Gyr of cluster galaxies, the relatively small age increase (~0.71 ± 0.4 Gyr) found in our sample galaxies seems to suggest that a slow environmental process such as starvation is the dominant quenching pathway. Our results provide new insights into environmental quenching effects spanning a large range in cosmic time (~5.2 Gyr, z = 0.26–1.13) and demonstrate the power of using a kinematically derived infall time proxy.