Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “charging analysis”

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 55 records · Page 3

JuTrack: A Julia package for auto-differentiable accelerator modeling and particle tracking

Efficient accelerator modeling and particle tracking are key for the design and configuration of modern particle accelerators. In this work, we present JuTrack, a nested accelerator modeling package developed in the Julia programming language and enhanced with compiler-level automatic differentiation (AD). With the aid of AD, JuTrack enables rapid derivative calculations in accelerator modeling, facilitating sensitivity analyses and optimization tasks. Here we demonstrate the effectiveness of AD-derived derivatives through several practical applications, including sensitivity analysis of space-charge-induced emittance growth, nonlinear beam dynamics analysis for a synchrotron light source, and lattice parameter tuning of the future Electron-Ion Collider (EIC). Through the incorporation of automatic differentiation, this package opens up new possibilities for accelerator physicists in beam physics studies and accelerator design optimization.

43 PARTICLE ACCELERATORS↗

Surface Wettability Analysis from Adsorption Energy and Surface Electrical Charge

Surface wettability is determined by the attraction of a liquid phase to a solid surface. It is typically quantified by using contact angle measurements at mineral surfaces in the case of the flotation of mineral particles. Contact angle research to describe wettability has been investigated at different scales by sessile drop measurements, molecular dynamic simulation, and atomic force microscopy. In this study, the density functional theory (DFT) was employed for predicting the surface free energy and contact angles of a well-known hydrophobic phyllosilicate mineral talc and a well-known hydrophilic phyllosilicate mineral muscovite based on the calculated interfacial energy and surface charge. The results revealed that the predicted contact angle at the atomic scale was larger than the experimental value, and identified two interactions: electrostatic interaction and hydrogen bonding, between the hydrophilic muscovite surface and the water layer, while a water-exclusion zone of 3.346 Å was found between the hydrophobic talc surface and the first water layer. This investigation gives a new perspective for wettability determination at the atomic scale.

Geochemistry & Geophysics↗

Microstructure Scale Lithium-Ion Battery Modeling, Part IV: The Representativity of Microstructure Parameters and Electrochemical Response

Lithium-ion battery electrochemical models require an accurate description of the electrodes microstructures to be predictive, that can be achieved through nanoscale imaging. Such observations are however limited by their field of view (FOV), as they provide only a subset of the whole electrode volume that does not necessarily represent the whole electrode microstructure heterogeneity, and therefore can bias the analysis. A representativity analysis has been performed on the microstructure parameters and, in a novel way, on the full cell electrochemical response to evaluate the predictions representativeness, and thus relevance, of a microstructure scale electrochemical model. The microstructure parameter deviation propagations to the electrochemical response have been quantified for different charge rates. This defines a threshold for the microstructure parameters FOV for a desired maximum deviation of the electrochemical response. Electrochemical model shows cell representative section areas are increasing with C-rate, due to higher in-plane heterogeneities, indicating larger FOVs are required specifically for fast charge modeling. Representativity analysis determines a cell FOV of 144.4 × 154.4 μm 2 is large enough to establish a convergence on the representative section areas for low-intermediate C-rate (≤2.5 C), therefore positively concludes on the model representativeness for these rates, but is not large enough to conclude for higher rates.

25 ENERGY STORAGE↗

Chiral population analysis: a real space visualization of X-ray circular dichroism

The microscopic understanding of probing and controlling molecular chirality is of considerable interest. Numerous spectroscopic techniques are capable of monitoring molecular asymmetry and its consequences, ranging from the infrared to the X-ray regime. Resonant X-rays have long been used to investigate local atomic sites within molecules thanks to the localized nature of core electronic transitions. These techniques can be used to determine the extent to which chirality is a local versus a delocalized property. However, how to systematically partition dichroic contributions from the point of view of electronic structure simulations remains an open question. Here, we introduce the concept of chiral population analysis that connects chirality to the atomic orbital picture. In analogy with Mulliken population analysis, which assigns charges to atomic orbitals, chiral populations allow the dichroic response to be distributed among the participating atomic orbitals. This decomposition can be further visualized in real space by representing it in terms of isosurface plots, providing an intuitive way to connect the dichroic response to its origins. Thus chiral population analysis can be particularly useful to assess the extent to which a given electronic transition is sensitive to chirality as a local or global feature of the molecular geometry.

36 MATERIALS SCIENCE↗

Techno-Economic Analysis and Market Potential of Geological Thermal Energy Storage (GeoTES) Charged With Solar Thermal and Heat Pumps

In this project, we developed a techno-economic analysis (TEA) model that can be used to evaluate the viability of a proposed Geological Thermal Energy Storage (GeoTES) design. This MATLAB-based model integrates distinct subsystem models for the reservoir, wells, power cycle, and solar field to capture their distinct characteristics. It applies this approach in simulating GeoTES storage and dispatch operations for durations ranging from hourly to seasonal. Using cases studies based on GeoTES designs provided by industry partners - Premier Resource Management (PRM) and EarthBridge Energy - we validated the TEA model estimations of system performance and costs (such as thermal and electrical power/energy inflow and outflow, capital costs, and levelized costs of energy and storage) for both concentrating solar thermal (CST) and Carnot Battery (CB) pairings with GeoTES (CST-GeoTES and CB-GeoTES). For the CST-GeoTES case, the model was validated against the proposed system designed by PRM. It showed good agreement with PRM's estimations when well and pump costs derived from PRM's estimations were used. When GETEM-based costs were used, there was a slight overprediction due to GETEM's project/site agnostic assumption of these costs. From a sensitivity analysis perspective, the levelized cost of electricity (LCOE) of the CST-GeoTES case was most sensitive to well flow rate and the charging temperature. An optimal design scenario resulted in an LCOE of 0.11 $\$$/kWhe. CST-GeoTES can also provide a source of heat to meet seasonal demands. With 12-hour and 24-hour levelized cost of heat (LCOH) of 0.018 $\$$/kWhth and 0.022 $\$$/kWhth, respectively, CST-GeoTES could be competitive in the California market with an average industrial price of natural gas in California between 0.041-0.047 $\$$/kWhth. The levelized cost of storage (LCOS) for CST-GeoTES depends on the energy storage duration. Although the LCOS is relatively higher for shorter durations (e.g., ~0.50 $\$$/kWhe for 1 hour of storage), it is an order of magnitude lower (0.06 $\$$/kWhe) for longer storage durations and competitive with lithium-ion batteries (beyond 12 hours of storage) and molten-salt thermal energy storage (beyond 32 hours). Energy. Three options were explored and applied to the EarthBridge case study: (1) A Carnot Battery design using R125 working fluid with both hot and cold storage; (2) A Carnot Battery design using R125 working fluid with only hot storage; (3) A Carnot Battery using a commercially available heat pump with carbon dioxide (CO2) working fluid and hot storage only. The CB-GeoTES with cold storage only had a slight (round-trip) efficiency advantage over the system without (43.4% vs. 42.8%). This is because the cold storage is limited by the freezing point of water, so the cold storage is not much colder than the environment. The system using commercially available technologies was the least efficient - partly because different cycles were used in the heat pump (CO2) and heat engine (binary cycle) which leads to some inefficiencies. Using the commercially available design, the levelized cost of energy (LCOS) from the model (0.10 $\$$/kWhe) was higher than that estimated by EarthBridge (0.068 $\$$/kWhe). This is because of the low round-trip (38.7%) efficiency of the commercially available design. Sensitivity analysis reveals that the model is most sensitive to electricity price. Including electricity price in the TEA for CB-GeoTES leads to an increase in LCOS from the base value to 0.25 $\$$/kWhe. To determine storage sites suitable for GeoTES, we gathered and analyzed geological, petrophysical, and geophysical data of oil and gas reservoir and aquifers in California and Texas. We down-selected possible sites based on cut-off values for site characteristics (e.g., reservoir temperature, formation thickness, permeability, porosity, depth, and brine salinity) and preliminary costs. Using this approach, the Carrizo-Wilcox, Yegua-Jackson, and Dockum brackish aquifers in Texas were identified as having the highest suitability. Similarly, in the central California region, the White Wolf, Belridge South Tulare, and Belridge South Reef Ridge were the most suitable. Going further, we assessed the storage potential in the selected sites. To do this we developed distributions of reservoir characteristic data and applied a Monte Carlo-based analysis to account for intrinsic uncertainty in the acquired data. The analysis revealed that the Carrizo-Wilcox aquifer had the highest storage potential with a mean capacity of 554 TWhth (i.e., 63 TWhe). The estimated capacity serves as an upper limit of storage potential given that not all fields in the basin will be developed. We participated in multiple outreach activities including conference presentations, panel session discussions, and the facilitation of a GeoTES workshop at the NREL Golden campus.

15 GEOTHERMAL ENERGY↗

One-loop analysis of β decays in SMEFT

We perform a loop-level analysis of charged-current (CC) processes involving light leptons and quarks within the Standard Model Effective Field Theory (SMEFT). This work is motivated by the high precision reached in experiment and Standard Model calculations for CC decays of mesons, neutron, and nuclei, and by a lingering tension in the Cabibbo universality test. We identify the SMEFT operators that induce the largest loop-level contributions to CC processes. These include four-quark and four-fermion semileptonic operators involving two third-generation quarks. We discuss the available constraints on the relevant effective couplings and along the way we derive new loop-level bounds from K → πν$\overline{v}$ on four-quark operators involving two top quarks. We find that low-energy CC processes are quite competitive with other probes, set constraints that do not depend on flavor-symmetry assumptions, and probe operators involving third-generation quarks up to effective scales of Λ ≃ 8 TeV. Finally, we briefly discuss single-field ultraviolet completions that could induce the relevant operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE↗

Search for Quasiparticle Scattering in the Quark-Gluon Plasma with Jet Splittings in 𝑝⁢𝑝 and Pb-Pb Collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The ALICE Collaboration reports measurements of the large relative transverse momentum (𝑘 𝑇 ) component of jet substructure in 𝑝⁢𝑝 and Pb-Pb collisions at center-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV. Enhancement in the yield of such large-𝑘 T emissions in head-on Pb-Pb collisions is predicted to arise from partonic scattering with quasiparticles of the quark-gluon plasma. The analysis utilizes charged-particle jets reconstructed by the anti-𝑘 T algorithm with resolution parameter 𝑅 = 0.2 in the transverse-momentum interval 60 < 𝑝 T,ch,jet < 80 GeV/𝑐. The soft drop and dynamical grooming algorithms are used to identify high transverse momentum splittings in the jet shower. Comparison of measurements in Pb-Pb and 𝑝⁢𝑝 collisions shows medium-induced narrowing, corresponding to yield suppression of high-𝑘 𝑇 splittings, in contrast to the expectation of yield enhancement due to quasiparticle scattering. The measurements are compared to theoretical model calculations incorporating jet modification due to jet-medium interactions (“jet quenching”), both with and without quasiparticle scattering effects. These measurements provide new insight into the underlying mechanisms and theoretical modeling of jet quenching.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗

Exploring the Effects of Varying Pre-Chamber Geometry in a Heavy-Duty Natural Gas Optical Engine under Dilution Conditions

Pre-chamber combustion is an advanced ignition strategy that has been shown to enhance spark ignition (SI) combustion stability in natural gas (NG) engines by providing distributed ignition sites from turbulent jets and enhancing main-chamber turbulence. Pre-chamber combustion has been proven especially advantageous compared to SI in ultra-lean and dilute operating conditions. This work involves experimental investigation of the effects of varying passive pre-chamber nozzle configuration on pre-chamber and main chamber combustion under simulated exhaust gas recirculation (EGR) dilution (0 and 20%) conditions in a heavy-duty, single-cylinder, optically accessible NG engine at stoichiometric fuel-air ratio. Pre-chamber nozzle configurations include four pre-chambers with constant nozzle area to pre-chamber volume ratio (A/V) with different nozzle sizes and orientations and one configuration with larger nozzles. The optical engine is operated in a skip-fire sequence consisting of 18 motored cycles followed by two consecutive fired cycles to elucidate the effect of combustion residuals (internal EGR) on combustion evolution. Pressure-based diagnostics are used to monitor pre-chamber and subsequent main chamber combustion, and optical diagnostics include high-speed OH* chemiluminescence to visualize the development of pre-chamber jets and the resulting ignition of the main chamber charge. Heat release analysis of the in-cylinder pressure data indicates that a faster pre-chamber pressure rise does not always translate into faster main-chamber combustion. The pre-chamber with the smallest nozzle diameter produced the highest pre-chamber pressure rise and fastest combustion under non-diluted conditions. However, dilution delays the main chamber ignition for smaller nozzles despite a comparable rise in pre-chamber pressure compared to configurations with larger nozzles. This effect is more pronounced for cycles with in-cylinder combustion residuals in addition to external dilution. Additionally, it was observed that pre-chambers with swirling nozzles have a faster pressure rise in the pre-chamber and main chamber under dilute conditions. Optical diagnostics suggest that the main reason for the delay between the pre-chamber pressure rise and main-chamber combustion lies in jet quenching and delayed re-ignition, which can even lead to misfire if jets emitted from small nozzles combined with dilution fail to re-ignite.

Dhotre, Akash [University Of Minnesota-Twin Cities↗

Electric Vehicle Supply Equipment Cybersecurity Through Emulation

As the grid evolves, it is paramount to understand the risks that cyberattacks pose before assets are deployed. Leveraging the ARIES Cyber Range, NREL has created a platform to conduct analysis of EV charging protocol cybersecurity to understand the risks and impacts that cyberattacks may pose to critical infrastructure.

bug bounty prize↗

Electric Vehicle Supply Equipment Security Solutions

The EVs@Scale cybersecurity pillar deep dive focuses on the three main tasks that NREL has in the consortium: (1) Analysis of EV charging mobile applications, (2) Cybersecurity risk assessments of EVSE through DER-CF, and (3) PKI for EVSE.

ADVANCED PROPULSION SYSTEMS,MATHEMATICS AND COMPUT↗

Theoretical model of dielectrophoretic water vapor condensation during multiphase flow of electrically charged water droplets in air

Water droplet growth in air, influenced by vapor diffusion and electrical charging, has significant implications for indoor air quality, cloud formation, and aerosol dynamics. Here, this study presents a comprehensive theoretical model that describes the growth dynamics of electrically charged water droplets under external electric fields, incorporating diffusion, electrostatic, and dielectrophoretic forces to enhance understanding and practical applications in air dehumidification and environmental control. A theoretical model was established to analyze the growth dynamics of the charged water droplets. The approach extended the thermodynamic framework of the Modified Kelvin-Thomson equation by integrating electrostatic and dielectrophoretic effects, enabling a comprehensive evaluation of vapor concentration near the surface of charged droplets. A sensitivity analysis was performed to assess the roles of droplet size, electric field intensity, surface charge, and vapor concentration on condensation rates. The study progressed beyond the classical diffusion-limited growth of neutral droplets, addressing the behavior of charged droplets in electrohydrodynamic environments. The model highlighted the significance of dielectrophoretic drift flux for dipolar vapor molecules, especially in micron-sized droplets ranging from 10 to 100 µm. Findings demonstrated that, for smaller droplets subjected to elevated electrical charge, average growth rates varied from 0.3 to 1.6 µm/s if the air relative humidity increased from 50 to 100 %. Analytical results are corroborated by numerical simulations, which exhibited strong concordance, thereby validating the proposed theoretical approach. This research introduced an innovative and scalable analytical methodology for electric field-enhanced droplet growth, offering valuable insights into the advancement of energy-efficient air dehumidification, water desalination, and atmospheric water harvesting technologies.

42 ENGINEERING↗

Single-Bunch Instabilities Driven by Space Charge During Low-Energy Cooling at Injection in the EIC Hadron Storage Ring

This paper presents a simulation-based study of single-bunch dynamics at the injection energy of 23.8 GeV for protons in the EIC hadron storage ring, focusing on the impact of space-charge–driven instabilities. The analysis demonstrates that at this energy, the proton bunch experiences significant transverse space-charge forces, which can reduce the stability margin in the presence of the geometric and resistive wall impedance. Various collective effects were considered, with particular attention to the nonlinear nature of the transverse space charge. To stabilize the beam, high chromaticity and octupoles were introduced and their effects analyzed using the ELEGANT code. The results provide a quantitative assessment of the stability thresholds and offer guidance for the machine design and operational strategy at injection.

43 PARTICLE ACCELERATORS↗

MAS-CCD: New technique for measuring low-level charge content based on the multiple amplifier architecture

Low-noise detectors are a key technology for the next generation of astronomical instruments aimed at spectroscopy of faint objects and the search for exoplanets. In this context, the multiple-amplifier sensing charge-coupled device (MAS-CCD) emerges as a promising technology for future scientific instruments. A critical parameter affecting the performance of these devices is spurious charge, produced by the clocking of the gates. Its measurement is typically challenging with existing methods. In practice, the optimization of this parameter often relies on empirical procedures that require significant time and careful consideration of the trade-off with full-well capacity. In this work, we present a new technique to estimate spurious charge based on covariance analysis of the output amplifiers of the MAS-CCD, which measures the same charge packet in different amplifiers at different times. The method enables fast and precise measurements of spurious charge under operating conditions where conventional approaches are difficult to apply. We develop the theoretical framework of the method and validate the model through simulations. The results demonstrate the feasibility of this approach and suggest that it could serve as a basis for reliable large-scale characterization of sensor performance.

Gamero, Miqueas E. [Buenos Aires, CONICET; Fermila↗

EVI-EnSitePy (Electric Vehicle Infrastructure – Energy Estimation and Site Optimization Tool in Python) [EVI-X Modeling Suite] [SWR-25-07]

EVI-EnSitePy is a comprehensive agent-based tool designed for the analysis and design of high-power charging sites, encompassing a wide array of site agents including Electric Vehicles (EVs), chargers, energy storage units (ESS), renewable energy resources (DER), and loads. This versatile tool offers diverse functionalities and a modular modeling approach, allowing detailed configuration of agents based on power ratings, port numbers, energy capacities, demand requirements, charger interfaces, and flexibility to customize the tool for project specific goals. By simulating agent interactions and employing various metrics, EVI-EnSitePy enables the assessment of site performance, exploration of energy management systems (EMS), and implementation of innovative EV charging policies. Utilizing EV charge schedules and arrival states, the tool performs thorough charging site simulations, with outputs consisting of agent and site-level power profiles and statistical metrics. Employing a tree graph structure, EVI-EnSitePy supports nested site structures and power distribution modeling. The tool's ability to generate charging schedules deterministically or via stochastic analysis further enhances its versatility. Through its features and capabilities, EVI-EnSitePy offers a powerful platform for informed decision-making in the realm of high-power charging site design and operation.

Jackson, Derek [National Renewable Energy Laborato↗

Indications for Freeze-Out of Charge Fluctuations in the Quark-Gluon Plasma at the LHC

The D-measure of net-charge fluctuations quantifies the variance of net charge in strongly interacting matter. It was introduced over 20 years ago as a potential signal of quark-gluon plasma (QGP) in heavy-ion collisions, where it is expected to be suppressed due to the fractional electric charges of quarks. Measurements have been performed at RHIC and LHC, but the conclusion has been elusive in the absence of quantitative calculations for both scenarios. We address this issue by employing a recently developed formalism of density correlations and incorporate resonance decays, local charge conservation, and experimental kinematic cuts. We find that the hadron gas scenario is in fair agreement with the ALICE data for $\sqrt{{𝑠}_{\textrm{NN}}}$ = 2.76 TeV Pb–Pb collisions only when a very short rapidity range of local charge conservation is enforced, while the QGP scenario is in excellent agreement with experimental data and largely insensitive to the range of local charge conservation. A Bayesian analysis of the data utilizing different priors yields moderate evidence for the freeze-out of charge fluctuations in the QGP phase relative to hadron gas. The upcoming high-fidelity measurements from LHC Run 2 will serve as a precision test of the two scenarios.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scalable Truck Charging Demand Simulation for Cost-Optimized Infrastructure Planning

This project developed a scalable, high-resolution model to simulate medium- and heavy-duty (MHD) electric truck charging demand and assess its impact on grid infrastructure. Using generative modeling, simulation, and cost optimization, the project delivered an end-to-end software pipeline and a library of 96 real-world scenarios for the Dallas–Houston megaregion. We demonstrated a modular architecture for transportation and grid modeling, implemented cost-optimized infrastructure planning methods, and quantified grid capital, operational, and environmental costs across a wide range of truck electrification scenarios. The results have been adopted by major utility stakeholders and contributed to regional planning efforts.

24 POWER TRANSMISSION AND DISTRIBUTION↗