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.
This presentation provides an overview of REopt Lite, and how it can be used to evaluate the economic and resilience benefits of hybrid systems.
We present the numerical formulation, verification, and performance portable implementation of an immersed interface method for microstructure scale electrochemical modeling of batteries. The innovation in this approach is the resolution of chemical species and electrostatic potential discontinuities at active interfaces without the use of interface conforming unstructured grids. A unified formulation on Cartesian grids for all domains (electrodes and electrolyte) is used with interfacial flux conditions applied using volume fraction or “color” function gradients. We have developed one dimensional and two dimensional test cases with analytic solutions for electrochemical modeling using which we verified the consistency and accuracy of our scheme. Our solver is also validated against solutions from a macroscale model and an unstructured multi-subdomain solver for a full lithium ion cell. We then demonstrated the utility of our solver on an image-based complex battery electrode microstructure at high charging rate. Our technique also exhibits good scalability on distributed memory architectures using central processing units (CPU), with problem sizes up to 1.8 billion degrees of freedom and with 5400 ranks. Initial performance studies of our open-source performance portable solver showed about 70 times speed up using a graphics processing unit (GPU) compared to single compute core for a problem with 4 million cells.
This study investigates the role of the background meridional moisture gradient (MMG) on the propagation of the Madden–Julian oscillation (MJO) across the Maritime Continent (MC) region. It is found that the interannual variability of the seasonal mean MMG over the southern MC area is associated with the meridional expansion and contraction of the moist area in the vicinity of the MC. Sea surface temperature anomalies associated with relatively high and low seasonal mean MMG exhibit patterns that resemble those of El Niño–Southern Oscillation. Therefore, by contrasting the years with anomalously low and high MMG, we show that MJO propagation through the MC is enhanced (suppressed) in years with higher (lower) seasonal mean MMG, although the effect is less robust when MMG anomalies are weak. Column-integrated moisture budget analysis further shows that sufficiently large MMG anomalies affect MJO activity by modulating the meridional advection of the mean moisture via MJO wind anomalies. Our results suggest that the background moisture distribution has a strong control over the propagation characteristics of the MJO in the MC region.
Abstract To‐date, only two UIC Class VI permits have been issued by the US Environmental Protection Agency. We illustrate a four‐phase workflow to first identify regional storage resources and then down‐select sites to yield permit‐ready locations that can accept and store large volumes of CO 2 . Specific permit requirements should guide objectives and define deliverables of respective workflow phases. In the first phase we used available regional data and screened structure and injection zones to locate resources that match CO 2 volumes planned to be captured. Available data were also used to assess presence and depth of usable groundwater, the key resource being protected via permitting. We then used advanced, closed‐form, analytical solutions (EASiTool) to estimate CO 2 injectivity into each hydrologically connected injection compartment. In the second phase we acquired and conditioned additional wireline logs and leased available seismic datasets. We interpreted the depositional systems from wireline well‐log character and mapped sandbody geometry to interpolate injection and confining‐zone distribution. Using available data, we mapped faults and locations of freshwater and overpressure (or other capacity‐limiting geologic parameters) in more detail. In the third phase, we used the augmented geologic data to develop a static model for the selected area, extracted the areas of highest interest, and generated and ran dynamic (flow) models. In a fourth phase, we reduced major uncertainties identified in earlier phases. Our case study indicates that to complete preparation of a permit application requires (1) improved lithologic characterization information (thicknesses and horizontal and vertical connectivity) and (2) better definition of poorly defined local faults. © 2023 The Authors. Greenhouse Gases: Science and Technology published by Society of Chemical Industry and John Wiley & Sons Ltd.
Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.
A high performance infrared (HPIR) system was developed and demonstrated for the infrared absorption analysis of {sup 235}U and {sup 238}U isotopes in uranium hexafluoride gas samples. Sweeping the quantum cascade laser light source over the spectral range and sampling via a high-rate analog-to-digital converter provided 0.0005 cm-1 spectral resolution, which allowed for high-precision measurements of the isotopic peak shift. A data analysis method was developed using principal component analysis to predict the isotope weight % content of {sup 235}U. The HPIR precision, accuracy, and error were evaluated for a wide range of isotopic ratio samples (0.287 – 93.7 weight % {sup 235}U), and the results were compared to the International Target Values (ITVS) set forth by the International Atomic Energy Agency (IAEA) for non-destructive and destructive analytical techniques. The method meets or surpasses the IAEA ITVs for non-destructive analysis of samples with isotopic content of depleted to highly enriched. The results also demonstrated the capability of the HPIR system to correctly predict the {sup 235}U weight % content of a mislabeled sample whose isotopic distribution was validated by mass spectroscopic measurements. The HPIR measurement is nondestructive and, thus, allows for confirmatory analyses of the exact sample at a designated IAEA lab if higher-resolution or a certified analysis is needed. (author)
Abstract In this article, we examine the validity range of the Effective Field Theory (EFT) description of high-energy Drell–Yan processes at the LHC. To this purpose, we consider explicit mediators that contribute to these processes in the s - and t -channels, comparing their effects in Drell–Yan distributions with the ones obtained by matching onto the corresponding EFT. We determine the conditions for the EFT results to accurately describe these scenarios. In particular, we explore the impact of including dimension-eight $$(d=8)$$ ( d = 8 ) operators in the faster convergence of the EFT series, at the analytical and numerical level, considering contributions to the cross section up to the square of $$d=8$$ d = 8 EFT operator insertions. Moreover, we discuss the possible implications of clipping LHC data and illustrate results for a specific New-Physics scenario motivated by low-energy flavor data.
Mega-ampere class pulsed power machines drive intense currents into small volumes to study high energy and density environments. Power lost during these events is a difficult and paramount problem to solve. For example, facilities such as Sandia National Laboratories’ Z machine experience meaningful power loss, which can be linked to non-linear ohmic heating at high currents (i.e., 26 MA on Z) leading to thermal desorption of contaminants and subsequent shunt plasma formation. Characterizing and understanding this type of thermal desorption is key to design optimizations necessary to minimize current loss, which will be even more important for next generation pulsed power. This type of characterization requires the ability to identify and determine concentration of analytes with nanosecond resolution given the pulse width of Z is on the order of 100 ns. This report summarizes progress on a small exploratory project focused on investigating options to meet this challenge using mass spectrometry. The main focus of these efforts utilized an Energy and Velocity Analyzer for Distributions of Electric Rockets intending to determine how quickly transient data could be resolved. This probe combines an electrostatic analyzer with a Wien velocity filter (ExB) to obtain ion energy and velocity distributions. Primary results from this exploratory project indicate significant additional work is needed to demonstrate a nanosecond time scale mass spectrometer for this application and also highlight that alternative detection methods such as laser-based diagnostics should be considered to meet the need for ultra-fast detection.
The aim of this research program has been to develop model for qualitative and quantitative understanding of complex processes in small molecules, such as vibronic coupling, predissociation, photodissociation, and isomerization.
Consequence-driven Cybersecurity Analysis for Extreme Fast Charging Electric Vehicle Infrastructure Electric vehicle (EV) development and associated charging infrastructure are expected to advance rapidly. Thirty percent of all global vehicle sales may be EVs and hybrid EVs by 2025, and they will rely on increasingly sophisticated strategies for grid integration. Next-generation EV charging infrastructure is expected to include interconnected renewable resources, such as photovoltaic (PV) arrays and battery storage systems, along with grid-edge devices. Although distributed energy resources (DERs) are useful in several ways, such as peak shaving at high demand times and backup supply for added resilience, the integration of vehicle charging and DERs could create more avenues for cyberattack. Physical and/or remote access to EV charging station components, including charge ports, power electronics, controllers, and local generation (e.g., PV and energy storage) could be paths to cause power fluctuations, leading to altered operations at the charging station, escalated privileges to administrative systems, exfiltration of financial information (including personally identifiable information), and reduced grid stability. One compromised EV supply equipment component can open the door to a variety of exploitable vulnerabilities. Cloud computing and mobile application control have the potential to expand the threat surface to non-repudiation and firmware integrity challenges. Vendor clouds have access to hundreds of chargers, and if compromised, can scale the attack surface exponentially. The high power and voltage levels of xFC infrastructure (e.g., 400 kW at 1000- V DC) increase the hazards and ability to impact the grid and vehicles more than lower-power charging systems. Legacy communications systems and protocols could also put EV infrastructure at risk of cyberattacks requiring a robust patch management process. Communications networks link EVs and chargers to several stakeholders - including charging station operators, grid operators, vendors/manufacturers, and aggregators - who have both physical and network access to share information for control, monitoring, and analytics. Information in these networks that is vulnerable to compromise includes the state of charge, charging duration, payment information, electricity price, and load control. Analyzing and prioritizing these interconnections risks could help address cybersecurity related to data leakage and manipulation.
Carbon molecular sieve (CMS) membranes have impressive separation properties; however, both chemical and morphology structures need to be understood better. Here we characterize CMS with the simplest polyimide (PI) PMDA/pPDA (PMDA=pyromellitic dianhydride, pPDA=p-phenylenediamine), using FTIR, solid-state 15 N-NMR and 13 C-NMR, XPS, XRD, and Raman spectra to study chemical structure. Here we also compare gas separation properties for this CMS to a CMS derived from a more conventional PI precursor. The detailed characterization shows the presence of aromatic pyridinic, pyrrolic rings as well as graphitic, pyridonic components and a few other groups in both CMS types derived from the very different precursors. The CMS morphologies, while related to precursor and pyrolysis temperature details, show similarities consistent with a physical picture comprising distributed molecular sieving plate-like structures. These results assist in understanding diverse CMS membrane separation performance.
Non-uniform secondary current distribution at rotating disk electrodes (RDE) is a common problem when using resistive electrolyte media or large applied currents. In a recent publication, we have shown that auxiliary electrodes such as the ring of a rotating ring disk electrode (RRDE) can help suppress current non-uniformities at the disk enabling reliable electroanalytical measurements. However, this previous work considered linear kinetics where current distribution non-uniformities were moderate. In the present contribution, we consider current distribution non-uniformities encountered under Tafel kinetics. We show, for the case of 2,5-dihydroxy-1,4-benzoquinone (DHBQ) reduction, that optimally chosen ring conditions serve to provide effective shielding at the disk edge rendering the overall disk current distribution to be uniform. Numerical modeling and scaling analysis (using the Wagner number) are presented to aid a user in determining the optimal ring current density for achieving uniform disk current distribution under Tafel kinetics. This approach is especially useful when studying soluble-soluble redox transitions for which, unlike deposit distribution in electrodeposition, the current distribution non-uniformity is not visually apparent.
Mass transport limitations in zeolite catalysts pose major hurdles for their optimal performance in diverse chemical reactions. Most approaches to reduce these restrictions focus on the synthesis of either hierarchical or nanosized zeolites. Here we demonstrate that the existence of a siliceous, catalytically inactive exterior rim on ZSM-5 particles dramatically reduces the diffusion limitations, which leads to an enhanced catalyst lifetime for the methanol-to-hydrocarbon reaction. In this study, our findings reveal that binary inorganic and organic structure-directing agents enable a one-pot synthesis of Si-zoned ZSM-5 catalysts with diffusion properties that are characteristic of particles with a much smaller size. Operando ultraviolet–visible light diffuse reflectance spectroscopy reveals a marked reduction in external coking among Si-zoned samples. Molecular dynamics simulations to assess the diffusion of methanol and benzene in siliceous pores and in those with Brønsted acids reveal substantially reduced transport limitations in zoned regions, consistent with the improved catalyst activity of Si-zoned zeolites relative to that of ZSM-5 with a homogeneous acid-site distribution.
For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.
An accurate full-dimensional quantum state-to-state simulation of the six-atom title reaction based on first-principles theory is reported. Counterintuitive effects are found: Increasing the energy in the reactant’s CD 3 umbrella vibration reduces the energy in the corresponding product vibration. An in-depth analysis reveals the crucial role of the effective dynamical transition state: Its geometry is controlled by the vibrational states of the reactants and subsequently controls the quantum state distribution of the products. This finding enables generalizing the concept of transition state control of chemical reactions to the quantum state–specific level.
Multifidelity modeling has been steadily gaining attention as a tool to address the problem of exorbitant model evaluation costs that makes the estimation of failure probabilities a significant computational challenge for complex real-world problems, particularly when failure is a rare event. To implement multifidelity modeling, estimators that efficiently combine information from multiple models/sources are necessary. In past works, the variance reduction techniques of control variates (CV) and importance sampling (IS) have been leveraged for this task. In this paper, we present the CVIS framework—a creative take on a coupled CV and IS estimator for bifidelity reliability analysis. The framework addresses some of the practical challenges of the CV method by using an estimator for the control variate mean and sidestepping the need to estimate the covariance between the original estimator and the control variate through a clever choice for the tuning constant. Furthermore, the task of selecting an efficient IS distribution is also considered, with a view towards maximally leveraging the bifidelity structure and maintaining expressivity. Additionally, a diagnostic is provided that indicates both the efficiency of the algorithm as well as the relative predictive quality of the models utilized. Finally, the behavior and performance of the framework is explored through analytical and numerical examples.
Presentation provides an overview of how facility energy and water efficiency audits can provide resilience benefits, how NREL's REopt Lite web tool optimizes the economic and resilience benefits of distributed energy resources, and how FEMP's Technical Resilience Navigator assesses risk to a site's critical functions and prioritizes solutions.
The magnetic pre-sheath (MPS) width, L MPS , is a critical parameter to define the sheath potential, which controls the ion trajectory of low-Z species (D, T, He, and C), as well as the prompt re-deposition of high-Z species. In this work, to determine L MPS , we fabricated micro-trenches (30×30×4 µm) via focused ion beam (FIB) milling on a silicon surface and exposed them to L-mode deuterium plasmas in DIII-D via the Divertor Material Evaluation System (DiMES) removable sample exposure probe. The areal distribution of impurity depositions, mainly consisting of carbon, was measured by energy-dispersive X-ray spectroscopy (EDS) to reveal the deuterium ion shadowing effect on the trench floors. The carbon deposition profiles showed that the erosion was maximized for the azimuthal direction of φ = -40° (referenced to the toroidal magnetic field direction) as well as the polar angle of θ = 80°. A Monte Carlo equation-of-motion model, based on a collisionless MPS, was used to calculate the azimuthal and polar deuterium ion angle distributions (IADs) for a range of L MPS = k × ρ i , where ρ i is the ion gyro radius and k = 0.5-4. Then, gross erosion profiles were calculated by a Monte Carlo micro-patterning and roughness (MPR) code for ion sputtering using as input the calculated azimuthal and polar IADs for each value of k . Good agreement with the experimental C deposition profiles was obtained for the case k = 2.5-3.5. This result is consistent with a previous kinetic modeling prediction of k ~ 3, as well as previous analytical investigations that predicted the L MPS to be several ion gyro radii. A validation of theoretical sheath models supports its applicability to ITER and pilot plant divertors to successfully predict plasma-materials interactions.