SEARCH · Search NASA
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.
Thermionic conversion system performance analysis based on space charge theory and isothermal diode sheath theory
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.
The NASA-Goddard procedure for destructive physical analysis (DPA) of nickel-cadmium cells contains a method for analysis of residual charged nickel as NiOOH in the positive plates at complete cell discharge, also known as nickel precharge. In the method, the Ni(III) is treated with an excess of an Fe(II) reducing agent and then back titrated with permanganate. The Ni(III) content is the difference between Fe(II) equivalents and permanganate equivalents. Problems have arisen in analysis at NAVSURFWARCENDIV, Crane because for many types of cells, particularly AA-size and some 'space-qualified' cells, zero or negative Ni(III) contents are recorded for which the manufacturer claims 3-5 percent precharge. Our approach to this problem was to reexamine the procedure for the source of error, and correct it or develop an alternative method.
Single nucleon removal in relativistic and intermediate energy nucleus-nucleus collisions is studied using a generalization of Weizsacker-Williams theory that treats each electromagnetic multipole separately. Calculations are presented for electric dipole and quadrupole excitations and incorporate a realistic minimum impact parameter, Coulomb recoil corrections, and the uncertainties in the input photonuclear data. Discrepancies are discussed. The maximum quadrupole effect to be observed in future experiments is estimated and also an analysis of the charge dependence of the electromagnetic cross sections down to energies as low as 100 MeV/nucleon is made.
Single-nucleon removal in relativistic and intermediate energy nucleus-nucleus collisions is studied using a generalization of Weizsacker-Williams theory that treats each electromagnetic multipole separately. Calculations are presented for electric dipole and quadrupole excitations and incorporate a realistic minimum impact parameter, Coulomb recoil corrections, and the uncertainties in the input photonuclear data. Discrepancies are discussed. The maximum quadrupole effect to be observed in future experiments is estimated and also an analysis of the charge dependence of the electromagnetic cross sections down to energies as low as 100 MeV/nucleon is made.
Participation of U.S. scientists on the COPERNIC (COmplete Positive ions, Electrons and Ram Negative Ion measurements near Comet Halley) plasma experiment on the Giotto mission is described. The experiment consisted of two detectors: the EESA (electron electrostatic analyzer) which provided three-dimensional measurements of the distribution of electrons from 10 eV to 30 keV, and the PICCA (positive ion cluster composition analyzer) which provided mass analysis of positively charged cold cometary ions from mass 10 to 210 amu. In addition, a small 3 deg wide sector of the EESA looking in the ram direction was devoted to the detection of negatively charged cold cometary ions. Both detectors operated perfectly up to near closest approach (approx. 600 km) to Halley, but impacts of dust particles and neutral gas on the spacecraft contaminated parts of the data during the last few minutes. Although no flight hardware was fabricated in the U.S., The U.S. made very significant contributions to the hardware design, ground support equipment (GSE) design and fabrication, and flight and data reduction software required for the experiment, and also participated fully in the data reduction and analysis, and theoretical modeling and interpretation. Cometary data analysis is presented.
The onset of electrohydrodynamic motion associated with the imposition of an electric field across a thin layer of liquid has been investigated for the case in which the electrical conductivity varies linearly over the depth of the layer. The variation of the conductivity is due to concentration gradients in the charge-carrying solutes and its spatio-temporal evolution is represented by a convective-diffusion equation. When the viscous relaxation time is long compared to the time for charge relaxation, the analysis reveals that the neutral stability curves for the layer can be characterized by three dimensionless parameters: Ra(sub e) is equivalent to d(epsilon)E(sup 2, sub 0) Delta sigma/mu K(sub eff) sigma(sub 0), an electrical Rayleigh number; delta sigma/sigma(sub 0), the relative conductivity increment; and alpha, the transverse wave number of the disturbance. Here d is the thickness, epsilon is the dielectric constant, and mu is the viscosity of layer, E(sub 0) is the applied field strength at the lower conductivity boundary, and K(sub eff) is an effective diffusivity associated with the Brownian motion of the charge-carrying solutes. With viscous-stress-free boundaries, at which the electrical conductivity and the normal component of the electric field are prescribed, the critical Ra(sub e) is 1.504 x 10(exp 4) at a critical transverse wave number of 1.97 when Delta sigma/sigma(sub 0) is 10. As Delta sigma/sigma(sub 0) increases, the critical Ra(sub e) increases and shifts to shorter wavelength disturbances; the critical imposed field strength, however, passes through a minimum because the lower-conductivity boundary exerts a considerable stabilizing influence in the presence of steep conductivity gradients. Similar trends were obtained for liquid layers with rigid boundaries.
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.
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.
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.
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.
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.
Solar flare charged particle composition analysis with nuclear emulsion detectors in Nike-Apache sounding rockets
The radial variation of the interplanetary propagation of solar cosmic rays is investigated on a grand scale by performing a statistical analysis of charged-particle data obtained by Pioneers 10 and 11 from launch to Jovian encounter (i.e., at distances of 1 to 5 AU). A numerical propagation model is adopted which includes diffusion, convection, and adiabatic deceleration, as well as a variable power-law diffusion coefficient. The study of solar particle events is carried out by analyzing individually each of five key parameters that are uniquely defined in each solar particle event: the time to maximum flux starting from the initial release at the sun, the anisotropy at the time to maximum flux, the flux at that time, the time width of the event at half the flux at the time to maximum, and the decay time. Combined results for the proton energy ranges from 3.4 to 5.2 MeV and 24 to 30 MeV are presented, the average radial diffusion coefficients between 1 and 5 AU are determined for each energy range, and it is found that the combination of a near-impusive injection with a temporally constant spectral index produces a reasonably good fit to the high-energy data but not as good a fit to the low-energy data. Implications of the results obtained are discussed for cosmic-ray propagation theory, solar modulation studies, and interplanetary acceleration of solar particles.
The advent of the concept of holes and careful determination of mobilities have made it possible to realize that the potential difference between the contacts of an irradiated semiconductor in the absence of current flow is due to the difference in mobility of electrons and holes. In this paper, a preliminary analysis shows that charge neutrality and zero-current condition are compatible. This situation is clarified and it is found that the equation for the electric field E in the absence of current is not absolutely correct even under low-level injection conditions. In particular, it is shown that fairly large currents can be generated by suitable light sources, currents which are entirely due to the disparity between electron and hole mobilities.
A worst case assessment was performed on the Command Data System (CDS) multilayer printed circuit board and an output power transformer module in the power subsystem. An estimate of the Jovian environment during the 35 hour orbit insertion was supplied by JPL and used as an input to calculate the electron transport into the Galileo components. A radiation shielding analysis computer code, CHARGE, calculated the electron transport deposition trapped in the anticipated sensitive areas of the multilayer board and transformer module. Based on these trapped charge calculations electric fields were calculated between the identified isolated areas and the spacecraft ground. The results of the assessment of electrostatic discharge (DSD) in the CDS multilayer printed circuit board indicate that the probability of ESD in the FR4 is low. The probability of ESD in the components attached to the multilayer board, however, is uncertain based on a lack of prior experimental data.
Magnetic field models for Saturn based on data obtained by the Pioneer 11, Voyager 1, and Voyager 2 spacecraft are discussed. Saturn is distinguished form earth and Jupiter by the symmetry of its internal field about its axis of rotation; however, the strong periodic modulation of Saturn's radio emission is strong evidence of some departure from axial symmetry. Analysis of energetic charged particle absorption signatures observed by Pioneer and Voyager spacecraft are consistent with axisymmetric models of Saturn's magnetic field and demonstrate that any equatorial displacement of the dipole is limited to less than 0.01 R(S). Saturn's magnetosphere appears to be intermediate in configuration to those of earth and Jupiter. An equatorial ring current of about 10 million A, confined to a 5 R(S) thick annulus with inner and outer radii of 8 and 16 R(S), has a major effect on the geometry of the outer magnetosphere. Field lines are moderately but measurably stretched out in the equatorial plane.
Theoretically, charge exchange between solar wind alpha-particles and geocoronal hydrogen at the magnetosheath results in the production of helium ions in an excited state, i.e, He(2+) + H (right arrow) He(+*) + H(+) resulting in the emission of 30.4-nm radiation upon relaxation of the helium ions. It has been estimated that the magnetosphere was compressed to inside the orbit of IMAGE during the Halloween storm of 2003. This configuration should have placed the subsolar point of the magnetosheath well beyond the limb of the solar disk and allow EUV the potential of imaging this solar wind charge exchange interaction Analysis of EUV full-frame data acquired during 31 October 2003 (Day 304) indicates the detection of 30.4-nm photons in the likely vicinity of the magnetosheath either as a direct consequence of the emission proposed above or, provided the magnetosheath He(+) density was sufficient to exceed the EUV sensitivity threshold, via resonant scattering of solar UV. We present initial analysis of EUV imaging during the Halloween 2003 storm from the context of solar wind parameters, including composition, and from the perspective of the generation of simulated EUV imaging of the event both as an emitting and as a scattering process.