Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ILs”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Thermal Hydraulic Experimental Test Article: Second Year of Testing with Secondary Sodium System (Fiscal Year 2025 Final Report)

The Thermal Hydraulic Experimental Test Article (THETA) is currently installed in the Mechanisms Engineering Test Loop (METL) 28” test vessel #4. Both the primary and secondary sodium systems remain online to facilitate continued testing. This fiscal year, work was performed using a COMSOL Multiphysics magnetohydrodynamic model to characterize flow more accurately in the secondary electromagnetic flowmeters. Experimental campaigns were then performed to study the thermal hydraulic differences between sodium and water as a surrogate fluid in the THETA geometry as well as a study to better characterize and understand temperature oscillations that exist at the outlet of the core to the hot pool. A peer-reviewed article was published in the ASME Journal of Nuclear Engineering and Radiation Science detailing the THETA facility and providing an overview of a test that was performed with the primary and secondary system online [1]. Work continues to develop a database to house experimental THETA data to better facilitate collaboration with industry and laboratory partners for their use of the data for code benchmarking/validation. THETA remains fully operational and is positioned for continued testing in fiscal year 2026.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Multi-trait multi-environment genomic prediction strategies for Miscanthus sacchariflorus

Genomic selection holds the potential to serve as a strategic tool to enhance the genetic gain of complex traits in Miscanthus breeding programs. The development of improved cultivars requires their assessment for various traits across diverse environments to ensure suitable overall performance. Hence, the multi-trait multi-environment (MTME) genomic prediction (GP) models offer an opportunity to improve selection accuracy. This study aims to evaluate the potential of five GP models: (1) three MTME models including genotype-by-trait-by-environment interaction (G×E×T) and (2) two single-trait multi-environment (STME) models (with and without G×E interaction). A Miscanthus sacchariflorus population comprising 336 genotypes evaluated in three environments and scored for four traits (biomass yield YDY, total culm number TCM, average internode length AIL, and culm node number CNN) was analyzed. The predictive ability of the models was evaluated considering three cross-validation schemes resembling realistic scenarios (CV1: predicting new genotypes, CVP: predicting missing traits in a given environment, and CV2: predicting partially observed genotypes). On average, in all cross-validation schemes compared to the STME the predictive ability of the MTME models was 10% to 70% higher for TCM and AIL. On the other hand, for YDY and CNN, both STME models performed similarly or slightly better (between 5 to 64%) than the MTME models in most environments. While the MTME models were not successful for all traits when compared to their STME counterparts, MTME models improved the prediction of the performance of genotypes that were untested across environments or lacked trait information in a specific environment. Overall, our study suggests that MTME GP models can be implemented in Miscanthus breeding programs to improve the predictive ability of the complex traits, shorten breeding cycles, and accelerate selection decisions.

genomic prediction (GP)↗
DOE OSTI2025

Interphase Engineering Enabled by Using a Separator with Electrochemically Active Carbazole Polymers for Lithium-Ion Batteries

Separators are generally considered inert components in lithium-ion batteries. In the past, some electroactive polymers have been successfully applied in separator modifications for overcharge protection or as acid scavengers. This study highlights the first use of two “electroactive” carbazole polymers (copolymer 9-phenyl-9H-carbazole-phenyl [PCP] and poly(9-vinylcarbazole) [PVC]), which were each applied separately as coatings on the cathode-facing side of commercial Celgard 2325 separators, respectively, to enhance the cycling performance of 0.3Li 2 MnO 3 ·0.7LiMn 0.5 Ni 0.5 O 2 //graphite (LMR-NM//Gr) full cells through interphase engineering. The team observed an irreversible polymer oxidation process of the carbazole-functionalized polymers—occurring only during the first charge—for the modified separator cells, and the results were confirmed by dQ/dV analysis, cyclic voltammetry measurements, and nuclear magnetic resonance characterizations. During this oxidation, carbazole polymers participate in the process of interphase formation, contributing to the improved cycling performance of LMR-NM//Gr batteries. Particularly, oxidation takes place at voltages of ~4.0 and ~3.5 V when PCP and PVC are used as separator coatings, which is highly irreversible. Further postmortem examinations suggest that the improvements using these modified separators arise from the formation of higher-quality and more inorganic SEI, as well as the beneficial CEI enriched in LixPOyFz. These interphases effectively inhibit the crosstalk effect by reducing TM dissolution.

carbazole↗
DOE OSTI2025

Over three decades, and counting, of near-surface turbulent flux measurements from the Atmospheric Radiation Measurement (ARM) user facility

Processes mediating the coupling of terrestrial, aquatic, biospheric, and atmospheric systems influence weather, climate, and ecosystem dynamics via transfer of energy, momentum, water, and carbon (or other species). These exchange processes are quantified by measurements of near-surface turbulent fluxes. Understanding processes at these interfaces provides insight toward understanding and predicting current and future states within the Earth system. The Atmospheric Radiation Measurement (ARM) user facility has been conducting measurements of near-surface turbulent fluxes since the early 1990s at long-term fixed locations and shorter-term mobile deployments across the Earth. ARM has utilized two established methods for conducting these measurements: energy balance Bowen ratio (EBBR) and eddy covariance (EC). Primary measurements from the former include sensible and latent heat flux, while the latter also measures fluxes of momentum and carbon (primarily carbon dioxide, with methane fluxes measured at two locations to date). The EBBR systems have been deployed at 22 locations, and, to date, the EC systems have been deployed at over 50 sites, with plans for additional novel site locations in the future. Herein, the history, evolution, and key aspects of these instrument systems are documented, along with information on data quality assurance and post-processing, as well as best use practices. Additionally, three data validation experiments were recently conducted, and their key findings are summarized. Finally, ancillary datasets acquired by ARM, which can contextualize and aid interpretation of the near-surface turbulent flux measurements, are discussed. The datasets described herein include the eddy correlation flux measurement system: 30ECOR (https://doi.org/10.5439/1879993, Sullivan et al., 1997), 30QCECOR (https://doi.org/10.5439/1097546, Gaustad, 2003), ECORSF (https://doi.org/10.5439/1494128, Sullivan et al., 2019a), and associated AmeriFlux and Methane Value-Added Product, AMCMETHANE (https://doi.org/10.5439/1508268, Billesbach, 2011); the energy balance Bowen ratio system: 30EBBR (https://doi.org/10.5439/1023895, Sullivan et al., 1993) and 30BAEBBR (https://doi.org/10.5439/1027268, Gaustad and Xie, 1993); and the carbon dioxide flux measurement system: CO2FLX (https://doi.org/10.5439/1287574, https://doi.org/10.5439/1287575, https://doi.org/10.5439/1287576, Koontz et al., 2015a, b, c; https://doi.org/10.5439/1989774, https://doi.org/10.5439/1989776, https://doi.org/10.5439/1992202, Biraud and Chan, 2002a, b, c). These data can be found by searching the above data stream names at https://adc.arm.gov/discovery/#/results/ (last access: 8 September 2025).

Sullivan, Ryan C. [Argonne National Laboratory (AN↗
DOE OSTI2025

Vapor Phase Installation of CpCo(CO) 2 in MOF‐808

Metal-organic frameworks (MOFs), including MOF-808 present an opportunity for vapor phase installation of reactive centers that may serve as uniform sites for precision catalysis. The mechanism of cobalt installation in MOF-808 through atomic layer deposition (ALD) upon CpCo(CO )2 exposure is investigated through in situ FTIR spectroscopy complemented by density functional theory modeling. The role of subsequent H 2 O exposure in the hydrolysis of carbonyl and cyclopentadienyl ligands to complete the installation is also investigated. In situ FTIR study reveals that upon exposure of CpCo(CO) 2 to MOF-808 at 115 °C, a long-lived, stable intermediate is formed that presents carbonyl stretching vibrations similar to the undissociated precursor that are attributed to multiple chemisorbed carbonyl complexes. In MOF-808, cobalt species are initially absorbed near the BTC linker rather than directly on the nodes, revealing the importance of non-covalent interactions in the installation process. DFT suggests that H 2 O exposure promotes carbonyl elimination but leaves a cyclopentadienyl-capped cobalt that is more stable than the complete hydrolysis product. The simulations are consistent with cobalt installation in MOF-808 using CpCo(CO) 2 and H 2 O at 115 °C and further consistent with the lack of ALD thin film deposition of cobalt oxide thin films on planar supports.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tunable 3D Aerosol Jet Printing of Low‐Power Redox‐Gated Transistors with Multicomponent Inks

Printed hybrid electronics (PHE) offer a promising alternative for microelectronics fabrication, addressing some limitations of traditional subtractive manufacturing. Despite the versatility of PHE, particularly in the customization of printing inks, these devices have not yet matched the performance of silicon-based electronics due to challenges in gating mechanisms and operational stability. However, the potential of low-voltage redox-gating to achieve significant carrier modulations in correlated metal oxides remains unexplored in PHE. This study systematically investigates vanadium dioxide (VO 2 ) nanoparticles and redox inks, linking their organization in solution to their morphology, phase state, and properties in solid films and multilayered structures. Using an aerosol jet printer (AJP), a solid-state VO 2 transistor is fabricated, operating at just 0.4 V gating voltage. The printed VO 2 films demonstrate redox-modulated conductivity and consistent transistor behavior. The solid-state redox gating materials also provide long-term stability, with the device maintaining performance over 6000 cycles without degradation. These results highlight the potential of redox gating to enhance the application of functional nanoparticles in printed hybrid microelectronics, especially for flexible, low-voltage, and energy-efficient devices.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Particle‐Size‐Dependent Lithium‐Ion Transport in PEO/LLZO Composite Electrolytes

Lithium-metal batteries with solid electrolytes can deliver higher energy density and improved safety than conventional Li-ion batteries. Among solid electrolyte candidates, polymer/ceramic composite electrolytes are attractive because they combine polymer flexibility with the high ionic conductivity of ceramics. However, whether ceramic fillers synergistically reduce polarization losses in the polymer matrix remains unclear. A central unknown is the critical polymer/ceramic interfacial resistance (Rint,crit), below which adding ceramics lowers electrolyte overpotential. Here, we present the first macroscale model framework to quantify R int,crit for composite electrolytes based on polyethylene oxide (PEO) and Ta-doped Li 7 La 3 Zr 2 O 12 (LLZO). A 1D model for DC-polarization of tri-layer cells (PEO-LiTFSI/LLZO/PEO-LiTFSI) shows that LLZO surface functionalization reduces the PEO/LLZO interfacial resistance, consistent with electrochemical impedance measurements. Extending to a 2D composite model, we show notably that Rint,crit scales linearly with LLZO particle diameter and shifts toward experimentally accessible values (e.g., 28.8 Ωcm 2 ) as particle size increases. At fixed ceramic volume fraction, larger LLZO particles reduce the number of interfacial crossings, driving more current through the ceramic phase and lowering concentration polarization. In contrast, R int,crit is largely independent of ceramic volume fraction. These results demonstrate that ceramic filler-size engineering can enable synergistic, energy-efficient transport in polymer/ceramic composite electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2026

Nonequilibrium defect-phase nanostructures stabilized by irradiation in undersaturated Ni-Si nanocrystalline alloy

Nanocrystalline thin films of the undersaturated alloy Ni-8.5 at% Si were subjected to 2 MeV Ti irradiation at temperatures ranging from 450˚C to 550˚C. Correlative microscopy combining transmission electron microscopy (TEM), scanning-TEM and atom probe tomography (APT revealed that large dose irradiation at 450˚C of samples with initial grain sizes below 100 nm stabilized a novel nanostructure which surprisingly contained three co-existing phases, the γ face-centered-cubic (FCC) matrix, γ' L12 ordered precipitates on intragranular dislocation loops and Ni 31 Si 12 precipitates at triple junctions (TJs). In contrast, irradiation at 550˚C and irradiation of larger grain-size samples at 450˚C only produced a γ-γ' two-phase coexistence. Analysis of the three-phase nanostructure and phase field simulations indicates that radiation-induced segregation is most pronounced at TJs, thus triggering the formation of Ni 31 Si 12 precipitates. These incoherent precipitates, in turn, are expected to stabilize the grain size under irradiation. The results are generalized using the concept of driven defect-phases. It is suggested that the stabilization of driven defect-phases may impart radiation resilience by providing localized relaxation modes to the microstructure evolution during and after temporary perturbations in irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Scale-up study of enhancing algal growth through bioaugmentation with the indole-3-acetic acid producing bacteria Azospirillum brasilense

Algae-based wastewater treatment technologies can recover nutrients while generating algal biomass with diverse application potentials. While numerous bacteria have demonstrated the ability to enhance algal growth through symbiotic relationships in laboratory settings, scaling this benefit to larger applications remains challenging due to complex wastewater conditions and microbe-algae dynamics. This study explores the effectiveness of bioaugmentation with Azospirillum brasilense, a plant growth-promoting bacterium known for producing the phytohormone indole-3-acetic acid (IAA), in boosting algal productivity within an Algaewheel wastewater treatment system. Through a series of experiments conducted at lab-, pilot-, and full-scale levels, we evaluated the impact of bioaugmentation on algal growth and developed a bioaugmentation strategy for an Algaewheel system. The results indicate that a biweekly dosing, achieving a final bioaugmenting cell density of 0.2–4 × 10 9 cells/L, effectively doubles the biomass productivity in a residential subdivision wastewater treatment plant in Northern Illinois. The bacterial production of IAA was the primary mechanism driving this enhancement, as demonstrated by the similar growth and yield improvements observed with chemically supplemented IAA and bioaugmentation across all experimental scales. High-throughput sequencing of 18S rRNA genes revealed that the algal community in the bioaugmented tank exhibited more stable biodiversity than the control. Additionally, bioaugmentation improved nutrient removal efficiency during winter, highlighting its potential to enhance the overall performance and sustainability of algae-based wastewater systems. The developed bioaugmentation and chemical treatment methods offer further operational solutions for managing yield, optimizing biochemical profiles, and enhancing biofuel production potential.

Algae-bacterial symbiotic relationship↗
DOE OSTI2025

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗
DOE OSTI2025

Site-specific surface reactivity on MgO for atomic layer deposition via selective hydration

Atomic layer deposition (ALD) is a powerful technique for thin film synthesis, offering atomic-scale precision and conformality. While ALD of MgO has been widely studied for applications in energy storage and microelectronics, its potential as surface on which deposition may be selective and defects repaired remains underexplored. Here, we present a combined theoretical and experimental investigation of MgO surface hydration and its implications for targeted ALD growth using water and dimethyl aluminum isopropoxide (DMAI) as reactants. We perform density functional theory (DFT) calculations to examine molecular and dissociative H 2 O adsorption on MgO (100) terraces and step-edge sites, including pristine surfaces and those with Mg/O vacancies. Reaction Gibbs free energies are calculated under various conditions to quantify surface reactivity. Our findings reveal facet- and defect-dependent hydration behaviors that align with experimental ALD growth trends on MgO (100). This study provides a molecular-level understanding of MgO surface chemistry critical for optimizing ALD processes for thin film growth and defect repair.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Data-driven equation-free dynamics applied to many-protein complexes: The microtubule tip relaxation

Microtubules (MTs) constitute the largest components of the eukaryotic cytoskeleton and play crucial roles in various cellular processes, including mitosis and intracellular transport. The property allowing MTs to cater to such diverse roles is attributed to dynamic instability, which is coupled to the hydrolysis of GTP (guanosine-5'-triphosphate) to GDP (guanosine-5'-diphosphate) within the β-tubulin monomers. Understanding the equilibrium dynamics and the structural features of both GDP- and GTP-complexed MT tips, especially at an all-atom level, remains challenging for both experimental and computational methods because of their dynamic nature and the prohibitive computational demands of simulating large, many-protein systems. This study employs the “equation-free” multiscale computational method to accelerate the relaxation of all-atom simulations of MT tips toward their putative equilibrium conformation. Using large MT lattice systems (14 protofilaments × 8 heterodimers) comprising ~21-38 million atoms, we applied this multiscale approach to leapfrog through time and nearly double the computational efficiency in realizing relaxed all-atom conformations of GDP- and GTP-complexed MT tips. Commencing from an initial 4 μs unbiased all-atom simulation, we interleave coarse projective “equation-free” jumps with short bursts of all-atom molecular dynamics simulation to realize an additional effective simulation time of 1.875 μs. Our 5.875 μs of effective simulation trajectories for each system expose the subtle yet essential differences in the structures of MT tips as a function of whether β-tubulin monomer is complexed with GDP or GTP, as well as the lateral interactions within the MT tip, offering a refined understanding of features underlying MT dynamic instability. Furthermore, the approach presents a robust and generalizable framework for future explorations of large biomolecular systems at atomic resolution.

Wu, Jiangbo [University of Chicago, IL (United Sta↗
DOE OSTI2024

Salt matters: How ionic strength and electrolytes impact redox polymer reactivity and dynamics for energy storage

As the global demand for sustainable energy grows, redox-active polymers (RAPs) have emerged as promising materials for batteries due to their advantages in stability, ease of preparation, and low-cost processability. Despite factors traditionally known to impact polymer dynamics (e.g., temperature, viscosity, and structure), we posit that investigating the effect of ionic strength and/or supporting electrolyte types on the electrochemical performance of RAP systems is crucial, both in aqueous and nonaqueous systems. Here, we first highlight recent findings on RAP-electrolyte interactions, elucidating how their polyelectrolyte nature determines their redox activity. Then, we focus on strategies to enhance RAP performance for energy storage through ionic strength optimization and tailored electrolyte composition. These insights into the modulation of RAP reactivity provide a foundation for improving battery performance in both flow and stationary configurations, thus facilitating progress toward next-generation energy storage solutions.

25 ENERGY STORAGE↗
DOE OSTI2026

Comparison of ammonia with methanol, liquefied natural gas and conventional marine transportation fuels through life cycle cost and emissions analysis

This work evaluates ammonia as a potential marine fuel for a SUEZMAX tanker and compares it with methanol, liquefied natural gas, and conventional fuel oils. The motivation arises from the need to identify low-emission, cost-competitive fuel options that can reduce greenhouse gas emissions from international shipping. The central hypothesis is that ammonia produced from renewable energy sources can achieve lower well-to-wake greenhouse gas emissions with varying life cycle costs based on the region. Life cycle assessment and techno-economic analysis were performed for a thirty-year vessel lifetime on two representative trade routes: from Saudi Arabia to Japan and from Saudi Arabia to the Netherlands. Four ammonia production pathways were assessed: natural gas, natural gas with carbon capture, natural gas pyrolysis, and renewable electricity–based synthesis. Results show that wind-based ammonia produced in Saudi Arabia achieved the lowest life cycle well-to-wake greenhouse gas emissions, between 0.58 and 0.64 million metric tons, among all fuels when using regional grid process electricity. With renewable process electricity, ammonia produced from natural gas pyrolysis in Saudi Arabia showed comparable emissions of 0.37 to 0.44 million metric tons with wind-based ammonia of 0.37 to 0.43 million metric tons. Liquefied natural gas exhibited the lowest life cycle cost, between 402 and 412 million United States dollars, and the only negative carbon abatement cost, ranging from −277 to –322 United States dollars per metric ton of greenhouse gas, compared with high sulfur fuel oil. The findings indicate that renewable ammonia offers a promising long-term pathway for reducing shipping emissions, while liquefied natural gas remains the most cost-effective option in the near term.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Impact of non-irrigation on 1G and 2G bioethanol potential of oilcane feedstock: A field to fuel pipeline study

This study evaluates the bioethanol potential in response to irrigation (IR) and non-irrigation (NIR) of oilcane (OC) during a seasonal drought prior harvest. The juice was extracted through mechanical pressing of stems and fermented by Ethanol Red® yeast to produce first-generation bioethanol. Hydrothermal pretreatment followed by enzymatic hydrolysis of bagasse was performed to produce monomeric sugars from structural carbohydrates. The hydrolysates were fermented with engineered yeast for second-generation bioethanol production. The irrigated oilcane juice (276.3 ± 8.9 g/L) constitutes higher sugar concentrations than non-irrigated oilcane juice (236.5 ± 2.2 g/L). The enzymatic hydrolysis of IR-OC and NIR-OC pretreated bagasse yielded similar concentrations of 247.5 ± 2.22 and 249.7 ± 4.98 g/L fermentable sugars. Industry-relevant bioethanol titers of ≥99 g/L and ≥75 g/L were achieved from juice and hydrolysates, respectively. Therefore, the non-irrigation regime did not impact the 1G and 2G bioethanol titers. However, the overall bioethanol yield can be lower due to the reduction of stem yield (8 %) per hectare.

09 BIOMASS FUELS↗
DOE OSTI2025

Bat activity at ecovoltaic solar energy developments in the Midwestern United States

As global photovoltaic (PV) solar electricity generation continues to increase, some PV sites are co-prioritizing electricity generation and ecosystem function (“ecovoltaics”) to align renewable energy development with biodiversity conservation. Thus far, positive responses of plant and insect communities to native habitats at ecovoltaic sites have been observed, but there has been little research on bat responses to ecovoltaic designs in the U.S. We conducted passive ultrasonic monitoring in 2023 and 2024 at 12 solar sites and paired reference sites (agricultural fields) to investigate bat activity responses to ecovoltaic facilities in the Midwestern U.S. We found that average weekly overall bat activity was approximately 50?% higher within ecovoltaic sites than reference sites in the first half of the monitoring season. We also found species-specific differences in bat responses to ecovoltaic sites, with Hoary Bats showing higher activity on ecovoltaic sites throughout most of the monitoring season, Big Brown Bats showing higher activity on ecovoltaic sites during the first one-third of the monitoring season, and Silver-haired Bats showing no difference in activity between ecovoltaic sites and reference sites. There were no weeks in which bat activity was statistically greater on reference sites, suggesting that bats in the Midwestern U.S. do not avoid ecovoltaic solar sites. Rather, our results suggest that ecovoltaic sites in this region may provide early season (May-June) habitat at a time of year when resources may be limited in the surrounding landscape. These findings support a growing body of evidence on the positive ecological outcomes of ecovoltaic designs. Further investigations of the types of bat calls being recorded at PV sites and relationships with insect prey abundance are needed to understand the underlying drivers of species-specific responses to PV developments.

14 SOLAR ENERGY↗
DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗